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CHAPTER  1  BACKGROUND 


1 . 1  INTRODUCTION 

Adaptive  pi i.  sed  array  radar  (APAR)  represents  one  of  the  most  demanding 
signal  processing  problems  that  merits  consideration  of  the  use  of  advanced 
signal  processing  techniques  such  as  optical  processing.  In  this  report,  we 
summarize  the  results  of  .  one  v  ar  advanced  study  of  the  use  of  optical  pro¬ 
cessing  for  APAR.  The  basic  concepts  and  the  three  major  optical  signal  pro¬ 
cessing  (OSP)  approaches  used  are  reviewed  in  Sects.  1.3-1. 5.  In  Sect.  1.2, 
we  provide  an  overview  and  summary  of  the  relevant  APAR  issues.  In  Chapters 
2-4,  each  of  the  new  OSP  techniques  is  considered  in  more  detail.  Our  summary 
and  conclusions  are  advanced  in  Chapter  5. 

In  Sect.  1.2,  the  salient  issues  of  APAR  as  they  apply  to  our  OSP  solutions 
are  reviewid.  For  a  more  global  description  of  APAR  research,  see  |1].  In  Sect. 
1.3,  we  outline  our  coherent  optical  correlator  (COC)  algorithm  and  approach. 

In  this  system,  an  optical  processor  is  used  to  compute  the  far  field  noise  distri¬ 
bution  in  angle  (or  space)  and  time  (or  frequency)  for  a  phased  array.  OSP  systems 
using  acousto  optic  (AO)  transducers  are  used  to  achieve  the  necessary  system  real¬ 
ization  in  the  desired  performance.  A  post  processor  can  then  compute  the  set  of 
adaptive  weights  to  apply  to  the  elements  of  the  antenna  to  null  this  far  field 
noise  pattern.  The  details  of  this  system  and  our  recent  research  on  it  are  in¬ 
cluded  in  Chapter  2. 

In  Sect.  1.4,  we  outline  and  highlight  our  second  OSP  technique  for  APAR  pro¬ 
cessing.  This  system  is  basically  an  optical  matrix-vector  multiplier  with  an 


2 


electronic  feedback  loop.  The  matrix-vector  multiplier  is  realized  by  a  linear 
array  of  LEDs,  fiber  optic  interconnections,  a  mask  and  a  linear  photo  detector 
array  with  parallel  output.  The  addition  of  an  electronic  feedback  system  results 
in  an  optical  system  capable  of  solving  a  general  matrix-vec tor  equation.  We 
refer  to  this  as  an  iterative  optical  processor  (IOP).  For  APAR  applications, 
this  system  is  used  to  compute  the  set  of  adaptive  weights  W  given  the  covariance 
matrix  M  in  the  steering  vector  S.  This  system  has  been  fabricated  and  evaluated. 

It  is  described  in  Chapter  3. 

The  third  and  final  OSP  system  for  APAR  uses  multiple  input  wavelengths  as  an 
adjunct  to  the  IOP  system  of  Chapter  3.  The  resultant  wavelength  diversity  processor 
(WDP)  is  highlighted  in  Sect.  1.5  and  described  in  detail  in  Chapter  4.  It  greatly 
enhances  the  capability  of  the  IOP  system  and  allows  the  use  of  many  new  algorithms 
and  optical  data  processing  operations  for  diverse  APAR  and  signal  processing  appli¬ 
cations. 

1.2  APAR  OVERVIEW  AND  REVIEW 

A  simplified  block  diagram  of  an  adaptive  phased  arra>  processor  is  shown  in 
Fig.  1.1.  The  basic  concept  is  to  multiply  each  of  the  received  signals  v^  by  an 
appropriate  weight  w^  and  to  then  sum  these  products  to  produce  the  output 

N 

E  =  "  w  v  .  (1.1) 

n  n  n 

If  the  same  weight  is  applied  to  all  elements,  the  beam  formed  is  normal  to  the 
array  and  described  by 

sin  (n  Mxsin  °) 

/'  d  sin  f'  (  ’ 

sin  ( - x - ) 


E  (ft)  =  K 


(1.2) 


1 


where  d  is  the  element-element  spacing,  "  is  the  steering  vector  of  the  array, 
and  N  is  the  number  or  elements.  By  varying  the  phases  of  the  weights  w^,  the 
beam  can  be  steered  to  different  directions.  This  is  achieved  by  selecting  the 
phases  such  that  a  signal  incident  from  the  desired  scan  angle  adds  in  phase 
across  the  array. 

By  adjusting  the  amplitudes  and  phases  of  the  w^,  the  sidelobe  levels  can  be 

decreased  and  the  effects  of  other  noise  sources  in  the  antenna's  field  of  view 

can  be  reduced.  In  such  APAR  systems,  nulls  are  placed  at  angles  and  frequencies 

in  the  antenna  pattern  corresponding  to  different  noise  sources.  The  adaptive 

control  is  achieved  by  separate  adaptive  loops  on  each  of  the  antenna  elements  as 

shown  in  Fig.  1.2.  The  steering  signals  s*  indicate  the  direction  in  which  the 

n 

arrav  is  steered  and  hence  the  location  of  the  antenna's  main  beam.  The  z  values 

n 

are  the  correlation  of  the  received  arrav  signals  v  and  the  output  g  =  '  w  v  . 

n  n  n 

When  w  t  ©  g  is  large,  u^  and  hence  change  rapidly.  The  purpose  of  an  adaptive 
array  is  to  reduce  the  noise  in  g.  When  u^  in  a  given  channel  is  large,  the  cor¬ 
responding  weight  w^  has  a  larger  effect  on  reducing  the  residue  noise  in  g. 

For  receiver  noise  only,  the  same  for  all  channels,  the  weights  w^  all  approach 

the  same  value.  If  one  w  is  larger,  the  corresponding  w  v  term  in  g  is  larger 

n  n  n 

and  this  will  cause  larger  z  and  u  values  for  that  channel,  which  will  decrease 

n  n 

w^.  Thus,  an  adaptive  array  achieves  uniform  illumination  only  for  receiver  noise 
alone.  When  the  received  signal  energy  is  less  than  the  interference  and  noise 
energy,  the  adaptive  loops  attempt  to  minimize  the  input  power  (subject  to  the  steer¬ 
ing  vector  constraint).  If  the  noise  in  g  is  approximately  zero,  the  z  values  will 

n 

be  small  and  the  w  will  be  relatively  constant, 
n 

The  weights  can  be  described  bv 


(1.0 


i ocps 


mui 


Figure  1.1 


F iguro  1  .  2 


ARRAY 

it  i’MKTS 


Simplified  block  diagram  of  an  adaptive  array 


Diagram  of  one  type  of  adaptive  control  loop. 
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The  low  pass  filter  with  time  constant  r  can  be  described  by 


i  u  +  u  =  z  , 
n  n  n 


(1.4) 


and 


z  =  v*  g  =  v*  w  v  , 
n  n  n  n  n 


(1.5) 


These  equations  describe  the  response  of  each  of  the  N  adaptive  loops.  The  vari¬ 


ation  of  w^  depends  on  the  inputs  from  all  N  channels  and  on  the  weights  of  the 


other  N-l  channe’s.  In  matrix  form,  we  describe  the  loops  by 


W  =  G  (S*  -  U) 

:  l’  +  l!  =  Z 
Z  =  V*  (V  -  W), 


(1.6) 


where  V  is  a  row  vector  and  W  etc.  are  column  vectors.  In  terms  of  average 


values 


Z  =  v"  V1  W  =  MW, 


(1.7) 


where  M  is  the  covariant  matrix 


M. .  =  Avg 
il 


(v;  v.).  (i.8) 

If  S  is  a  constant  steering  vector,  l'  =  -W/G  and  the  loop  is  described  by 


i  W/G  +  (M  +  J/G)  W  =  S  , 


(1.9) 


where  I  is  the  identity  matrix.  Equation  (1.9)  describes  a  set  of  N  different 
equations  that  describe  the  average  response  of  the  array  weights  W.  The  values 
W  depend  upon  the  external  noise  (through  M)  and  the  control  loop  parameters  (S  , 
G  and  ).  The  conventional  approach  is  to  solve  (1.9)  for  W  and  then  apply  these 


W  to  the  received  array  signals. 


The  covariance  matrix  M  is  fundamental  to  all  AFAR  processing  theory  since 
it  describes  the  noise  environment.  The  diagonal  terms  in  M  are  a  measure  of  the 
power  in  each  channel,  whereas  the  off-diagonal  elements  of  M  describe  the  direc¬ 
tion  of  arrival  of  the  noise.  The  rate  of  convergence  of  (1.9)  depends  upon  the 
noise  environment.  The  steady  state  solution  to  (1.9)  is  found  with  W  =  0 
(assuming  G  s--  1)  to  be 

WM  =  S*.  (1.10) 

The  desired  weights  W  can  then  be  found  from  S  and  M  ^  as 

W  =  M_1  S*.  (1.11) 

Equation  (1.11)  can  be  obtained  by  inverting  the  matrix  M,  and  forming  the  indi¬ 
cated  vector /mat rix  product.  Computation  of  M  requires  2  N  samples  at  the  signal 
frequency  to  yield  acceptable  statistics. 

Various  correlation  loop  processors  exist.  The  maximum  SNR  circuit  is  the 
most  popular.  It  uses  the  Widrow  least  mean  square  circuit.  It  requires  an 
initial  estir.iate  for  W,  computation  of  the  gradient  at  different  surface  points, 
etc.  until  the  minimum  concave  surface  is  obtained.  The  system  in  Fig.  1.2  uses 
a  typical  Howe  1 ls-Appe Ibaum  loop  in  which  the  residue  is  fed  back  and  correlated 
with  each  of  the  received  signals.  The  filtered  outputs  are  proportional  to  the 
gradient.  When  subtracted  from  the  steering  signal,  they  yield  the  adaptive  weights 
A  modified  random  search  technique  with  many  variations  also  exist. 

1 • '■*  COHERENT  OPTICAL  CORRELATOR  SYSTEM 

in  the  COC  system,  we  compute  the  angular  and  temporal  noise  distribution  of 
the  far  field  antenna  pattern  using  acousto  optic  correlators.  A  dedicated  hard¬ 
ware  digital  post  processor  can  then  compute  the  set  of  adaptive  weights  from  the 
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optically  computed  N 


or  N 
m  m 


(■-) 


i  m) 


noise  power  N  of  the  various  noise  sources  at  ancles 
m 


,  f  (noise  power  versus  angle  and  frequency  distribution)  patterns. 


The  details  of  this  system  and  our  recent  research  on  it  are  included  in  Chapter 
2.  in  this;  section,  we  highlight  the  svstem  and  algorithm  concepts  involved  for 
background  purposes. 


Since  the  basic  concept  bv  which  we  determine  the  power  N  and  the  location 

m 

angle  ■  of  each  noise  source  is  quite  different  fr"m  the  conventional  approach, 
we  briefly  review  the  system  philosophy  here. 


We  assume  a  phased  array  radar  with  N  receiving  elements  with  parameters 

given  in  Table  1.1.  The  signal  s  received  at  antenna  element  n  is  correlated 

n 

with  a  reference  signal  (the  central  element  n  =  M/2  of  the  array  is  used  as  the 

reference  r  for  simplicity).  We  denote  the  total  received  signal  at  element  n 

due  onlv  to  the  noise  source  m  at  angie  '  bv  f  .  The  total  received  signal  at 

m  •  nm 

element  n  due  to  all  \  sources  is  then 


M 

f  =  :  f  . 
n  m  nm 


(1.12) 


With  respect  to  the  reference  element  r,  the  signal  received  at  n  can  be 
described  by 

f  (t)  =  f  ft--  )  (1.1 

n  m  nm  ^  nm  / 

where 


■  =  k(r-n)d  cos  (1.14) 

nm  m 

where  k  =  2"/'.  For  a  single  noise  source  m  at  '•  the  s.  ’1  received  at  elements 

m 

N/2  =  r  and  n  are  the  same;  they  are  simplv  delaved  in  time  bv  ■  .  With  M  noise 

nm 

sources  present,  we  simply  sum  over  m  to  obtain  the  total  signal  f  received  at 


arrav  element  n. 
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The  time  dependence  of  the  received  signals  was  explicitly  included  in 

(1.13)  to  allow  the  delay  in  arrival  time  of  the  signal  at  two  elements  of  the 

array  to  be  more  conveniently  written  as  :  .  Of  utmost  importance  is  the  fact 

nm 

that  the  delay  between  elements  r  and  n  is  also  a  function  of  the  angle  of  inci¬ 
dence  0  or  the  variable  m  of  the  noise  source.  This  delav  i  is  known  in  ad- 
n  '  nm 

vance  (for  a  given  array  element  n  and  direction  “  ),  thus  we  know  where  to  look 

for  a  given  channel  n  and  angle  r'  . 

m 

The  reference  signal  r  sees  the  sum  of  M  noise  sources 

f  (t)  =  f  (t)  +  f  ,(t)  +  .  .  .  +  f  M(t).  (1.15) 

r  rl  r2  rM 

At  array  element  n,  we  find 

f  (t)  =  f  ./  t-r  \  +  f  Jt-T  ,)+  .  .  .  +  f  M  ft-  ).  (1.16) 

n  rll  nil  r2\  n2  /  rM  \  nm/ 

liquations  (1.15)  and  (1.16)  express  (1.13)  and  (1.14)  in  detail.  The  delay  depends 

on  the  element  n  as  well  as  the  angle  a.  The  separate  delays  !.  and  ,  are  re- 

tm  km 

lated  bv  : .  =  K:,  .  For  different  array  elements  n,  the  delav  between  the  signal 

tm  km  - 

received  at  successive  array  elements  (due  to  a  noise  source  at  a  fixed  •)  equals 
the  linear  relationship  shown. 


The  operations  required  on  the  reference  array  element  f  are  to  correlate  it 

with  all  N  other  received  signals  f  .  This  is  the  operation  we  propose  to  optically 

n 

perform.  We  denote  the  correlation  of  the  reference  signal  f  and  all  other  N  sig¬ 
nals  bv  C  where 
rn 

Cr„  '  Crl(-  nl)  +  c,2(‘- „2)  +  '  '  '  +  O'17’ 

Hach  of  these  correlations  f  (*)  f,,  f  (*)  f.,,  etc.  contains  manv  terms. 

r  1  r  2 

All  cross  terms  are  zero  if  the  separate  noise  sources  are  assumed  to  he  inde¬ 
pendent  or  not  coherent.  The  result  of  the  correlation  of  a  11  pairs  of  corrcla- 


id 
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Cions)  for  all  m  thus  yield  the  autocorrelation  of  all  M  noise  sources  with  each 

correlation  located  at  its  respective  t  value.  Since  t  is  related  to  f-  by  (1.14), 

m 

the  correlation  of  f^_  and  all  f  yields  the  desired  information,  the  power  and 
angular  distribution  of  all  noise  sources  in  the  antenna's  far  field  pattern,  i.e. 

M  (  f  V 

ml  m  J 

In  Fig.  1.3,  we  shown  the  general  form  for  such  a  multi-channel  correlation 
output.  Slit  integration  (with  a  specially-shaped  detector  or  by  post  processing) 
can  provide  the  desired  output  information.  In  Fig.  1.4a,  we  show  the  output  pattern 
from  such  a  system  in  the  form  of  Fig.  1.3.  After  slit  integration,  the  pattern 
of  Fig.  1.4b  results  with  two  peaks  at  the  correct  spatial  locations  corresponding 
to  the  angular  locations  of  the  two  noise  sources  and  with  the  amplitude  of  each 
proportional  to  the  power  of  the  corresponding  noise  source. 


correlator  for  AFAR. 


m'  ANGLE  of  noise  sourc 
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"lie  details  of  this  approach  to  AFAR  processing  are  described  elsewhere 
|  2 — 4 ] •  For  review,  we  highlight  how  this  correlator  approach  to  AFAR  processing 
relates  to  the  more  conventional  techniques.  In  Fig.  1.3,  the  horizontal  axis 
is  and  the  vertical  axis  is  array  element  number  n.  The  horizontal  row  of 

M 

spots  at  n  =  X  in  Fig.  1.3  represent  the  correlation  f  (*)  f  .  Since  f  =  ’  f 

•  r  w  N  r  m  rm 

and  f  =  1  f  are  each  composed  of  the  sum  of  the  signals  received  from  the  M 
X  m  Nm 

noise  sources  and  since  these  M  noise  sources  are  independent,  the  location  of  the  M 

spots  along  the  horizontal  line  at  n  =•  X  occur  at  the  '  values  given  bv  (1.14), 

rm 

These  locations  are  thus  proportional  to  the  M  angles  ■  ■  of  the  M  noise  sources. 

m 

The  intensities  of  these  M  spots  are  proportional  to  the  energies  present  in  each 

of  the  M  noise  sources.  These  M  spots  of  light  are  thus  the  M  autocorrelations  of 

the  M  noise  sources,  each  located  at  t  =  :  given  bv  (1.14),  and  thus  are  propor- 

nm 

tional  to  the  noise  source  angles, 
m 


We  now  show  the  equivalence  of  these  optical  and  analytical  f ormula t i ons .  Re¬ 


call  that  the  signals  received  at  array  elements  1  to  N'  are  f  and  the  portion  of 

n 

f  due  to  a  noise  source  at  is  given  bv  (1.12).  For  a  fixed  ,  each  element 
n  m  m 

n  of  the  array  sees  the  same  signal  f  delaved  in  time  bv  a  different  amount 

m 


at  each  detector,  i.e.. 


r  (t)  =  r  ft-  )  (i.i8) 

nm  m \  nm  I 

where  is  seen  to  depend  on  the  arrav  element  n  and  the  angle  "  of  the  noise 

m 

source.  From  (1.12)  and  (1.18),  we  can  write  the  received  signal  at  element  n 

as  in  (1.13)  as  the  supi,  over  the  M  source  angles,  of  the  signals  f  received  at 

nm 

arrav  element  m. 


t  he 


F  1  e pit'll l  (  i  ,  j ) 
product  of  thi¬ 


ol  the  covariance 
signals  received 


mat  r  i  x  >1 
it  arrav 


is  f  .  !  . 
i  ! 


e 1  emeu  t  s 


or  the  average  over  time  of 
i  anti  j.  We  write  this  as 


M 
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M .  .  =  f  .  f  .  =  /  f.(t)f.*(t)dt.  (1.19) 

i.l  i.l  i  I 


Substituting  (1.13)  into  (1.19)  we  find  (with 

M.  .  (■-  .  \f 

m  m\  im  j  m 


jm 


the 


( L )dt . 


zero 


reference) 


(1.20) 


Without  the  : .  zero  reference,  we  have 
jm 


M.  . 
1J 


(1.21) 


where  -  • .  .  Thus  element  (i,  j)  of  M  is  the  sum  over  M  of  the  correla- 

ij  im  j  m 

tions  of  the  received  signals  evaluated  at  ... 


For  a  fixed  target  angle  ‘‘  eacli  element  of  the  covariance  matrix  M  is  thus 

m 

the  autocorrelation  of  the  signal  (due  to  that  fixed  u)  at  1 .  -  .  Different 

im  j  m 

elements  in  M  correspond  to  different  array  elements  and  thus  different  shifts. 
For  each  source  m,  each  row  of  M  is  the  autocorrelation  of  that  noise  source  for 
all  ’  and  a  given  time.  Thus  M  =  M  (;,t).  As  shown  in  (1.21)  and  noted  above, 
the  covariance  matrix  M  is  equivalent  to  a  correlation  as  used  in  (1.17)  and  Figs. 
1.3-1. 4.  if  we  correlate  the  reeeived  signal  at  one  element  of  the  array  with  the 
received  signals  at  the  other  elements  of  the  array  we  thus  obtain  (1.21)  and  thus 
M.  In  the  optical  system,  is  continuous  of  course. 


The  above  analysis  has  shown  that,  when  written  in  our  correlation  terms,  the 
conventional  method  of  AFAR  processing  is  equivalent  to  ours.  In  Chapter  2.  we 
deseribe  our  reeent  resear  h  approaches,  algorithms  and  tecnniques  to  realize  such 
an  AFAR  processor  in  teal  time  with  presently  available  AO  transducers. 


1 r) 


1.4  ITERATIVE  OPTICAL  PROCESSOR  (  I  OP  )  CONCEPT 

The  basic  1  OP  system  is  a  vector /mat  r  i  x  multiplier  (lip,.  1.4)  in  which  the 
input  vector  is  an  estimate  of  tile  adaptive  weights  and  is  realized  as  the  out¬ 
put  from  a  linear  arrav  of  I.EIXs.  The  matrix  (1  -  M)  is  a  fixed  J— I )  mask.  Associ¬ 
ated  intermediate  optics  correctly  image  W  onto  (I  -  M>  and  I  -  M  onto  a  linear 

* 

output  detector  to  whose  outputs  S  are  added.  The  resultant  output 


W.  ,  ,  =  W.  ( X  — M )  +  S*  (1.22) 

i+l  i 

is  a  new  estimate  of  W.  it  is  fed  back  to  the  LED  input  and  the  vector /matrix 
multiplication  is  repeated  until 


W. 

t 


W 


i  +  1 


within  some  tolerance  >  .  When  this  occurs,  the  resultant  output 

W  =  M_1  S 


(1.23) 


(1.24) 


is  the  desired  solution  of  (1.23)  for  W^.  This  system  and  the  associated  iterative 
algorithm  used  are  discussed  in  more  detail  in  [5,6]  and  in  Chapter  3. 

1 . 5  WAVELENGTH  DIVERSITY  PROCESSOR  (WPP)  CONCEPT 

An  important  adjunct  to  the  10P  system  is  the  use  of  wavelength  as  an  extra 
input  variable  dimension.  The  resultant  system  is  shown  schematically  in  Fig.  1.6 
and  described  in  detail  in  Chapter  4.  We  refer  to  this  system  as  a  wavelength 
diversity  processor  (WDP) .  As  the  input  to  such  a  system,  we  use  three  (or  more) 
linear  arrays  of  laser  diodes  (LDs).  The  output  from  each  laser  diode  uniformly 
illuminates  t lit*  corresponding  row  of  the  mask  as  in  the  10P.  The  multi-color  light 
leaving  the  mask  is  then  split  into  three  or  more  original  input  colors  by  a  grating. 
The  corresponding  three  vector-matrix  proem ts  in  the  different  wavelengths  of  light 
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are  thus  formed  and  collected  under  different  corresponding  output  detectors. 
Various  alternate  versions  of  this  system  have  been  developed  and  demonstrated 
during  the  past  year.  These  and  the  new  algorithms  and  operations  possible  on 
such  a  U'DP  system  are  described  in  Chapter  4. 


Schematic  diagram  of  the  iterative  optical  processor  (IOP). 


MASK 


Schematic  diagram  of  the  wavelength  diversity  processor  (WDP) 
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CHAPTER  2  COHERENT  OPTICAL  CORRELATOR  (COC) 

2. 1  INTRODUCTION 

In  this  chapter  we  describe  our  recent  work  on  the  new  coherent  optica] 

correlator  (COC)  concept  for  APAR  signal  processing.  As  noted  in  Chapter  1, 

the  purpose  of  this  COC  system  is  to  obtain  an  estimate  of  the  far  field  noise 

distribution  N  ]  as  a  function  of  the  intensity  N  and  the  angular  location 
m\  ml  m 

'■  of  each  source.  In  our  new  research,  we  have  altered  the  original  space  inte¬ 
grating  (SI)  correlator  system  described  in  Chapter  1  to  a  time  integrating  (TI) 
system  ''Sect.  2.2).  This  new  system  permits  longer  integration  and  correlation 
times  a.d  “bus  better  noise  field  statistical  estimates.  We  have  also  changed 
the  real  time  spatial  light  modulators  (SLMs)  used  from  2-D  transducers  to  1-D 
acousto  optic  (AO)  cells.  This  was  done  because  AO  cells  are  more  easily  avail¬ 
able  than  2-D  SLMs  and  because  of  their  higher  bandwidth  and  superior  performance, 
compared  to  other  SLMs.  The  AO  cells  we  fabricated  during  this  past  year  are 
described  in  Sect.  2.1  together  with  the  initial  tests  we  were  able  to  perform  on 
these  devices. 

A  new  adjunct  antenna  concept  for  APAR  processing  on  the  COC  system  was  formu¬ 
lated  (Sect.  2.4)  and  a  digital  simulator  for  this  new  COC  system  was  written 
(Sect.  2.5).  A  new  hybrid  time  and  space  integrating  (TS1)  AO  system  architecture 
was  developed  (Sect.  2.6)  that  extends  the  COC  concept  to  wideband  receivers  and 
wideband  noise  sources.  This  new  TST  system  can  compute  both  the  angular  (spatial) 
and  temporal  (frequency)  distribution  of  the  antenna's  noise  field  pattern  in  parallel. 
Use  of  this  system  can  thus  enable  one  to  perform  adaptivity  in  both  time  anu  space. 

The  two  electronic  support  systems  (a  computer-driven  one  and  a  hardware  system) 
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that  we  are  assembling  for  the  AO  COC  systems  are  described  in  Sect.  2.7.  The 
issue  of  quantization  of  the  input  data  and  performing  complex  correlations  are 
key  features  in  the  use  of  both  systems.  In  Sect.  2.8,  we  present  two  new  post¬ 
processing  algorithms  for  the  narrow-band  and  wide-band  noise  cases.  The  post¬ 
processor  is  used  to  compute  the  adaptive  weights  from  the  spatial  and  temporal 
field  strength  distribution  of  the  noise  field  to  be  nulled.  In  Sect.  2.9,  exper¬ 
iments  on  the  AO  TI  correlator  for  residue  arithmetic  operations  are  presented 
(residue  arithmetic  was  one  of  the  candidate  techniques  for  APAR  that  we  studied 
in  our  earlier  report  [2]).  Initial  experimental  simulations  of  the  hybrid  TSI 
system  were  then  performed  (Sect.  2.10)  with  encouraging  results  that  demonstrated 
the  basic  concept.  The  status  of  the  COC  system  is  then  summarized  in  Sect.  2.11. 

2 . 2  TIME  INTEGRATING  (TI)  CORRELATOR  CONCEPT 

The  prior  COC  systems  we  considered  were  space  integrating  correlators  in 
which  the  correlation  was  performed  by  multiplication  in  the  frequency  domain 
and  the  resultant  correlation  was  displayed  in  space.  For  these  systems,  the 
correlation  of  the  two  signats  and  was  realized  as 

s,  ©  s2  =j.  S*""!  =  R]2(x),  (2.1) 

where  capital  letters  represent  the  Fourier  transforms  of  the  corresponding  space 
functions.  In  Fig.  2.1,  we  show  the  schematic  diagram  of  a  time  integrating  corre¬ 
lator  [7|.  This  is  the  basic  system  architecture  we  propose  for  the  COC  processing 
of  adaptive  array  data.  This  system  realizes  the  correlation  of  s^  and  s0  by  inte¬ 
grating  in  the  time  domain,  i.e. 

Sj  ©  s?  =  Sj  (t)  s*  (t-x)  dt  =  R j  2 (x) .  (2.2) 


20  - 


The  general  description  of  the  simplified  system  shown  is  quite  direct.  An 
input  source  (such  as  a  LED  or  LD)  is  time  sequentially  modulated  with  the  received 
signal  s^  from  one  antenna  element.  Thus  the  output  from  the  l.ED  as  a  function  of 
time  is  s^  (t).  The  output  from  this  point  source  is  expanded  to  uniformly  illu¬ 
minate  an  acousto  optic  (AO)  cell  fed  with  an  input  signal  s^  (t).  The  transmittance 
of  the  AO  cell  is  a  function  of  time  (t)  and  spatial  location  (x)  is  described  by 
(t-x/v).  The  light  distribution  leaving  the  AO  cell  is  thus  s^  (t)  s0  (t-x), 
where  x  =  x/v  and  where  v  is  the  velocity  of  the  acoustic  wave  in  the  cell.  For 
simplicity,  we  set  v  =  1  to  simplify  the  associated  mathematics.  The  AO  cell  is  then 
imaged  onto  the  output  plane  where  a  time  integrating  detector  is  placed.  This  forms 
the  integral  of  the  product  of  these  two  signals  in  time.  The  resultant  integration 
is  the  desired  correlation  in  (2.2).  It  is  displayed  in  space  x  at  the  output  plane. 
The  immediate  attractiveness  of  this  optical  signal  processor  for  the  adaptive  radar 
problem  is  that  long  integration  times  or  equivalently  long  duration  signals  of  large 
time  bandwidth  product  can  be  correlated.  This  system  has  a  low  range  delay  search 
window  (equal  to  the  aperture  time,  or  transit  time,  T  of  the  signal  across  the  acousto 

optic  cell).  However,  this  is  the  exact  case  that  occurs  in  adaptive  radar  processing 

(low  range  windows  and  long  integration  times).  Thus  this  basic  correlator  archi¬ 
tecture  appears  most  appropriate  for  the  adaptive  algorithm  required. 

The  system  of  Fig.  2.1  can  easily  be  extended  to  a  multi-channel  time  inte¬ 
grating  (TI)  correlator  as  shown  in  Fig.  2.2  by  replacing  the  input  light  source 

with  a  linear  array  of  hF.Ds  or  LDs.  The  resultant  output  plane  pattern  now  con¬ 

tains  N  1-D  correlations  of  the  reference  signal  s  with  all  N  other  received 
antenna  signals.  This  output  pattern  is  the  dual  of  the  one  obtained  in  our 
original  COC  system  described  in  Chapter  1.  However,  no  2-D  SLMs  are  required  in 
this  system  and  it  provides  longer  integration  times  and  hence  better  noise  estima¬ 
tion  than  the  space  integrating  system. 


2.3  ACOUSTO-OPTIC  (AO)  CF.LT.S 


Considerable  time  was  required  to  complete  fabrication  ol  the  AO  cells  to 
be  used.  Much  time  was  also  necessary  to  adequately  test  the  cells.  This  fea¬ 
ture  is  often  ignored  in  assembly  of  a  system.  Since  the  cell  s  specifications 
directly  affect  the  performance  of  the  system  and  provide  necessary  data  for  new 


improve 


d  cell  fabrication,  such  a  test  program  should  be  completed  in  the  next 


phase  of  this  program.  All  AO  cell  fabrication,  polishing  and  transducer  bondin'-' 
was  performed  by  V.'est inghouse  Corporation  (Pittsburgh). 

The  structure  of  the  AO  cells  fabricated  is  shown  in  Pig*  2.3.  Two  cells 

were  produced:  30  mm  long,  33  mm  high  and  6  mm  thick  (see  Fig.  2.4).  The  AO 

cells  consisted  of  a  LXB  transducer  bonded  to  the  TeO?  AO  material. 


Two  LNB  transducers  about  6  x  14  mm"  were  bonded  to  the  right  face  and  con- 
nected  in  series  as  shown.  An  8.3  x  10  mm  thick  x-cut  LNB  crystal  is  used. 

The  crystal  thickness  was  chosen  to  yield  the  desired  center  frequency  f^  [8]  of 
operation  from  t  =  =  6.6  x  l(2/2(4  x  102).  The  light  wavelength  used 

effects  the  center  frequency  and  bandwidth  of  the  TeO^  AO  material  used  as  shown 
in  Fig.  2.5.  We  selected  a  light  source  of  633  nm  (corresponding  to  the  HeNe 
laser  line)  and  thus  a  center  frequency  of  about  40  MHz  and  a  bandwidth  of  about 


20  Mllz  or  more  (before  matching)  , 


TeC^  was  selected  because  of  its  good  optical  quality,  ready  availability 
in  large  size  and  its  low  optical  and  acoustic  attenuation  [8].  TeO^  has  a  high 
M?  figure  of  merit  relating  diffraction  efficiency  ri,  acoustic  power  P  ,  height 
of  the  acoustic  beam  H  and  interaction  length  L  by  [8] 


5  M.,  >  2  H  -1  BP 
2  0  a 


For  TeO^ ,  M 7  =  795  compared  to  =  1  for  quartz.  Bandwidth  is  another  parameter 

2 

of  concern  in  AO  cells.  A  second  figure  of  merit  M1  =  M0  r<V  describes  the  band- 
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width  BW  performance  of  the  cell,  since  BW  r*  nV~.  TeO^  offers  an  M  =  13.1  versus 
Mj  =  1.0  for  quartz.  Other  materials  such  as  As^Se^  have  higher  Mj  =  204  but  are 
not  available  commercially.  Recent  work  may  change  this  condition,  but  moreso 
As ,;Sc is  of  use  only  in  the  IR  region  (0.9-11  pm). 

A  slow  shear  TeO?  crystal  was  thus  selected.  Its  acoustic  attenuation  in  this 
9 

mode  is  16  dB/usec  Hz“.  This  is  larger  than  the  value  obtainable  with  other  materi- 

2 

als,  but  it  is  much  lower  (better)  than  for  As^Se.^  (27.5  dB/usec-Hz  ).  For  mo¬ 
rel  1,  this  parameter  is  thus  negligible  since 

-2  2  2 

(16)  (50  usee)  (4  x  10  )  GHz  =  1.28  dB 

or  less  than  a  10%  loss  in  acoustic  power  across  the  length  of  the  line.  Few 

other  materials  (besides  KRS-5)  allow  operation  in  visible  light  with  good  optical 
quality  (KRS-5  is  quite  soft  and  thus  exhibits  poor  optical  quality).  Present 
technologi . al  capability  thus  dictates  the  use  of  TeO^ 


'V**. 


giirv 


Photograph  of*  orio  of  the  AO  rolls 
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The  design  and  selection  of  the  cells  were  completed  early  in  month  2.  The 
TeO^  material,  obtained  from  Crystal  Technology,  Incorporated,  arrived  in  month  6. 
The  optical  faces  of  the  cells  were  then  polished  at  Westinghouse  Corporation. 

2 

Interferograms  were  made  and  polishing  was  repeated  until  the  two  large  30  x  30  mm 
faces  were  flat  to  A/2  over  the  central  70%  of  the  area.  The  transducer  surface  face 
was  then  polished  to  A/2.  Gold  chrome  electrodes  were  deposited  on  the  TeO^  side 
of  the  LNB  and  the  LNB  was  then  cut  to  the  desired  size  and  bonded  to  the  crystal 
with  epoxy  and  then  lapped  from  its  original  0.010"  thickness  to  its  final  thickness 
of  3.3  mils. 

O 

The  epoxy  is  first  mixed  to  low  velocity  to  allow  a  1000  A  thick  layer  to  be 

deposited.  Small  droplets  of  epoxy  are  then  deposited  over  the  surface  (in  a  clean 

2 

room)  to  produce  750  gm/cm  deposition  and  baking  at  50°C  for  24  hours  and  cooling 
at  room  temperature  for  eight  hours  were  used.  The  top  Gold  electrodes  were  then 
deposited  and  wires  were  bonded  to  the  top  edges  of  the  LNB  to  obtain  the  series 
connection  shown  in  Fig.  2.3.  The  bottom  electrode  (Fig.  2.3)  is  deposited  to  the 
Te09  before  the  epoxy  bond  is  applied.  These  steps  were  completed  in  months  7  and 
8  at  which  time  the  cells  were  now  ready  for  initial  tests. 


Initial  AO  cell  tests  were  performed  in  month  9.  These  included  initial  verifi¬ 
cation  that  diffraction  occurred,  a  recheck  of  the  cell's  optical  quality,  and  point 
laser  tests  of  the  integrity  and  uniformity  of  the  LNB  bond.  These  latter  fre¬ 
quency  response  plots  are  necessary  to  design  the  impedance  matching  network.  These 
tests  consisted  of  illuminating  the  device  quite  close  to  the  LNB  with  a  point  laser 
beam.  For  different  positions  of  the  laser  beam,  the  frequency  of  the  input  signal 
was  varied  and  the  first  order  diffracted  light  intensity  versus  input  frequency  was 


measured.  (A  single  photo  diode  was  used  and  moved  to  different  positions  correspond- 


ing  to  the  different  diffracted  angles  for  the  different  input  frequencies). 

After  correctly  matching  cables  and  cable  lengths,  the  test  was  repeated  with 
a  swept  input  frequency  (20-60  MHz)  and  the  first  order  diffracted  pattern  was 
detected  on  a  reticon  linear  detector  array.  The  input  angle  oi  illumination 
of  the  laser  light  was  chosen  to  yield  the  f ] at  test  spec t rum  for  the  transducer¬ 
cell  combination  design. 

The  first  cell  had  a  nice  response  but  with  low  bandwidth  (Fig.  2.6a). 

The  second  cell  had  an  anomily  in  the  response  above  50  KHz  and  a  large  but  non-uniform 
bandwidth.  By  varying  the  laser  i 1 luminat ion  angle  and  by  proper  design  of  the 
impedance  matching  network,  we  can  reduce  the  central  peak  response  and  increase 
the  side  lobe  levels.  The  non-uniform  response  peaks  above  50  MHz  do  not  concern 
us  since  we  will  operate  the  device  at  35-40  MHz  with  a  10-20  MHz  bandwidth. 


We  now  consider  the  design  of  the  impedance  matching  network.  This  absorbed 
most  of  our  month  10  AO  cell  research.  The  objective  was  to  design  a  Rl.C  model 
of  tiu>  AO  device  (transducer  and  cell)  and  to  include  the  appropriate  Rl.C  impedance 
matching  network  to  match  the  frequency  response  impedance  data  on  the  system.  A 
parallel  model  yielded  poor  results,  with  a  fairly  constant  C  70-80  pF  value  but 
with  a  large  R  range  (  500-50  ohms)  from  25-45  MHz.  Thin  large  R  range 

would  make  it  difficult  to  couple  input  power  to  the  cell  and  would  result  in  re¬ 
flected  or  standing  waves.  A  series  model  was  used  and  gave  better  results  with: 

I!  =  15.6  olims,  C  =  71  pF  at  25  MHz;  R  =  18.7,  C  =  94  pF  at  40  MHz;  and'  R  = 

21  ohms,  C  =  126  pi"  at  45  MHz.  A  smaller  R  range  but  an  unacceptably  large 

0  range  occurtid.  It  was  thus  determined  to  cancel  the  capacitive  reactance  of  the 
system  (55  (dins)  at  14  MHz  with  an  0.25  ill  inductor  in  series  with  the  cell 

and  thus  increase  the  response  ol  tin1  cell  at  this  frequency  and  at  its  sidel.obes. 


In  Fig'..  2.6a  and  2.6b,  wt  show  the  frequency  response  of  the  first  AO  cell 
before  ,.nd  after  impedance  matching.  In  Fig.  2.6a,  we  see  a  resonance  peak  at 
f  =  47  MHz,  with  a  quite  low  2  MHz  bandwidth.  In  Fig.  2.6h,  the  impedance 
matching  network  is  seen  to  give  a  much  better  response  with  a  center  frequency 
of  31  MHz  and  a  3  dB  bandwidth  of  18  MHz  extending  from  22-40  MHz.  Similar  per¬ 
formance  was  obtained  for  the  second  cell. 


In  month  11,  further  initial  tests  were  performed  on  the  device  including  the 
measurement  of  near  field  point  laser  beam  probing  at  three  vertical  locations 
(center,  top  and  bottom  of  the  33  mm  cell  dimension  and  at  2  mm,  12  mm  and  24  mm 
from  the  LN'B  transducer).  This  data  provides  necessary  AO  field  strength  informa¬ 
tion  from  which  the  usable  vertical  cell  size  and  the  specific  vertical  portion  of 
the  cell  to  use  (from  uniformity  and/or  diffraction  efficiency  n  considerations)  emerges. 
The  strongest  acoustic  field  (and  the  largest  n  by  507)  occurred  along  the  center 
of  the  cell.  This  uneven  acoustic  field  restricts  the  vertical  aperture  of  use  to 
approximately  1  cm  or  less  in  the  center  of  the  cell. 

In  month  12,  with  the  integrity  of  the  optical  polish  ana  transducer  bond 
verified,  the  AO  cells  were  returned  to  Kestinghouse  to  h  vo  an  absorption  wedge 
ground  onto  the  left  edge  of  the  cell  as  shown  in  rig.  2.4.  The  angle  at  the  far 
end  of  the  cell  should  be  chosen  to  cause  reflected  waves  to  diffract  back  at 
different  angles  and  in  different  directions  than  the  original  acoustic  wave.  An 
angle  of  22° degrees  was  chosen  (as  a  compromise  between  not  decreasing  the  cell's 
usable  width,  not  damaging  the  crystal,  and  vet  producing  large  angular  and  direc¬ 
tional  differences  in  the  reflected  waves). 

The  completed  cells  arrived  just  as  the  present  contract  period  ended  and  thus 
only  several  preliminary  final  tests  were  possible.  These  were  conducted  after  the 
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final  contract  date  and  are  included  for  completeness.  The  inter f erogram  for  cell 
1  was  repeated  and  its  optical  quality  was  found  to  be  >/2  over  the  central  $07  of 
its  area.  Scatter  level  measurements  of  the  bulk  properties  of  the  unit  showed 
about  a  30  dB  scatter  level  measured  at  first  order.  N'ew  fre¬ 
quency  response  data  obtained  on  cell  1  verified  an  f  =  31  MHz,  a  3  dB  bandwidth 
of  18  MHz  from  22  MHz  to  40  MHz.  Unless  reflections  were  very  large,  we  would  not 
expect  the  f  and  the  bandwidth  of  the  cells  to  change.  This  was  verified  bv  the 
new  data  above.  A  reasonable  n  =  2 07  was  obtained  with  0.5  watt  of  input  electrical 
power  at  31  MHz.  Complete  tests  of  n  versus  f  for  different  input  electrical 
power,  input  light  angle  and  vertical  beam  position  will  be  done  in  a  latter 
phase  of  this  work.  Care  must  be  taken  not  to  burn  out  the  cell  with  large 
input  power.  0.1  watt  appears  to  allow  the  present  cells  to  perform  adequately. 

A  second  AO  cell  was  similarly  tested  with  an  0.25  uH  inductor  in  series. 

Its  center  frequency  (27  MHz)  and  3  dB  bandwidth  (14  MHz,  from  20  MHz  to  34  MHz) 
were  acceptable  and  compatible  with  the  intended  svstem. 

Point  probe  testing  of  both  cells  (or  Schliercv.  wtaging  of  the  acoustic 
wave)  were  repeated  and  n  versus  f  was  measured  as  a  function  of  cell  position 
(vertical  and  horizontal).  A  1 5  mm  cell  region  in  the  center  was  chosen  for 
use  because  it  demonstrated  the  maximum  uniformity  and  maximum  r,.  By  measuring 
the  frequency  response  of  the  cell  for  different  input  powers,  the  linear 
operating  region  of  cells  1  and  2  were  found  to  be  0.02  watts  to  0.2  watts. 

We  will  thus  operate  cells  1  and  2  over  this  input  power  range. 

Because  of  the  strong  near  field  effect  observed  (variation  in  the  intensity 
of  the  acoustic  interaction  versus  cell  position),  no  cell  phase  response  measure¬ 
ments  of  the  AO  transducer  could  be  made  at  present.  The  final  operating  specifi¬ 
cations  for  the  AO  cells  are  summarized  in  Table  2.1.  They  appear  quite  adequate 
for  the  intended  COG  system  for  AFAR. 
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TABLE  2.1 

AO  CELL  SPECIFICATIONS 

Parameter 

Cell  1 

Cell  2 

Time  Aperture  (I’sed) 

30  nsec 

30  nsec 

Center  Frequency 

31  MHz 

27  MHz 

Bandwi  dtli 

18  MHz 

14  MHz 

Time  Bandwidth  Product 

340 

420 

Optical  (Quality 

‘  /  2 

A/2 

Scatter  Level 

30  dB 

30  dB 

I)i  f  I  met  ion  Ef  f  ic  i  enc  v 

20" 

20% 

Input  Power 

0.1  Watt 

0.1  Watt 
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In  the  operation  of  the  T1  correlator  of  Fig.  2.1  and  its  modified  forms 
(Sect.  2.6),  care  must  be  taken  to  match  the  center  frequencies  f  and  01 
the  input  source  and  the  Art  cell.  In  our  initial  experiments,  we  will  use  a 
point  source  AO  modulator  in  place  of  the  LEI).  If  ^  f  ,  the  TI  correlator 
output  has  a  difference  term 

Sj(t)  s>2  ( t—  r )  cos  (2t'Aft)dt 

and  a  sum  term 

s  (  ( t )  s2  ( t- r  )  cos  2tt  f  2)  dt  * 

If  the  'f  or  fj  +  is  greater  than  the  bandwidth  of  s^,  then  the  integral 
appri'aclus  zero  since  integration  of  an  even  function  for  a  time  much  greater 
than  its  period  is  zero.  The  sum  term  above  will  always  be  zero  and  we  must  ad¬ 
just  f  to  equal  f0  and  thus  decrease  ,‘.f  to  less  than  0.1  of  the  bandwidth  of  s^  s^ 
or  else  electronic  detection  and  external  integration  are  necessary.  We  will  thus 
operate  the  Art  cells  slightly  off  from  their  central  resonance  frequencies  to  ensure 
■f  =  0  and  that  the  cosine  term  in  the  difference  term  above  approaches  1. 
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2 . 4  ADJUNCT  ANTENNA  CONCEPT 

The  initial  ( OC  concept  advanced  in  [2|  and  highlighted  in  Chapter  1  used 
correlations  to  obtain  the  angular  location  of  the  target.  This  operation  re¬ 
quires  correlation  measurements  to  determine  the  relative  time  delays  line-to- 
line  between  received  elements  of  the  phased  array.  A  numerical  analysis  of 
the  AO  cell  requirements  for  such  a  system  were  conducted.  We  assumed  a  con¬ 
servative  radar  center  frequency  of  1  GHz  (corresponding  to  a  )  =  0.3  m)  and 

R 

a  linearly  spaced  phased  array  with  spacing  d  =  1/2  =  0.15  m  between 

elements.  For  a  target  at  an  angle  0  =  60^*,  the  relative  time  delay  between 
adjacent  elements  is  thus 


,  =  (d/c)  sin  =  0.433  nsec.  (2.4) 

d 

If  an  AO  cell  with  1  .isec  aperture  time  (and  1-2  GHz  bandwidth)  were  used,  re¬ 
solving  the  time  delay  in  (2.4)  would  require  an  AO  cell  TBWP  of  1  Usee /0. 4 33 
nsec  2300.  Such  AO  cell  spec i f i cat  ions  [9]  are  possible  but  require  quite 
extensive  and  sophisticated  AO  cell  fabrication  aid  more  extensive  support  elec¬ 
tronics  (with  1-2  GHz  bandwidth).  Such  efforts  were  beyond  the  scope  of  our 
present  research  and  the  multi-channel  svstem  required  to  process  such  data  is 
more  complex  than  our  simpler  2  channel  svstem. 

For  these  reasons,  we  chose  to  consider  use  of  AO  cells  with  larger  time 
apertures  (e.g.  40-50  .'.sec),  lower  bandwidths  (10-20  MHz),  but  with  comparable 
TBWP  1000).  To  use  such  AO  devices  for  the  given  AFAR  problem,  we  must  in¬ 
crease  the  time  delay  to  be  resolved  e 1 ement-to-el ement  in  the  phased  array  by 
t he  AO  correlator.  This  resulted  in  the  new  system  concept  of  an  adjunct  antenna 
(with  one  element  widely  spaced  by  manv  wavelengths  from  the  main  phased  array). 
Fig.  2.7  shows  the  associated  adjunct  antenna  system  concept. 
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ADJUNCT 
ANTENNA 
(5>1  ELEMENTS) 


Figure  2.7  Adjunct  antenna  concept. 


Assuming  a  20  MHz  bandwidth  AO  cell  with  a  50  usee  aperture  time  (approximately 
equivalent  to  the  performance  of  our  AO  cells)  with  a  TBWP  of  1000,  the  minimum  time 
resolution  possible  is  50  usec/1000  =  50  nsec.  To  realize  such  a  time  delay,  we  must 
increase  the  spacing  d  between  adjacent  antenna  elements  (i.e.  in  this  case  the  spac¬ 
ing  between  the  phased  array  and  the  adjunct  omni  antenna  in  Fig.  2.7).  The  time  de¬ 
lay  :  in  (2.4)  can  be  rewritten  in  units  of  radians  as 

t  =  (2’>d/>)  sin  0.  (2.5) 

The  antenna  spacing  d  required  to  yield  a  '  =  50  nsec  for  a  target  at  0  =  60^ 

(sin  =  0.866)  is 

d  =  -c/sin  <•  =  17.3  m  =  57.6  A  at  f  =  1  GHz.  (2.6) 

If  the  input  signal  is  heterodyned  to  a  center  frequency  f  =40  MHz  or  1  =  7.5  m 
(corresponding  to  the  center  frequency  of  our  AO  cell,  a  normal  phased  array  antenna 
requires  d  =  ■/ 2  =  3.75  m  spacing  at  an  associated  TBWP  of  50,000/10.8  =  4629.  Since 
this  is  beyond  the  scope  of  presently  available  AO  cells,  we  consider  the  use  of  a 
adjunct  antenna  with  spacing 

d  =  nijn  c/sin  ■■  =  (so  x  10"9)  (t  x  l08\/0.866  =  2.3  A  =  17.3  m.  (2.7) 
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This  antenna  element  spacing  is  quite  realistic  and  is  compatible  with  the 
TBWP  =  1000  obtainable  directly  with  current  AO  technology. 

This  adjunct  antenna  system  thus  appears  capable  of  resolving  the  necessary 
50  nsec  minimum  time  delays  (between  the  signals  received  at  the  phased  array  and 
the  adjunct  antenna  element)  within  the  TBWP  =  1000  and  the  40  MHz  heterodyne  cen¬ 
ter  frequency  of  the  available  AO  cells.  As  we  shall  see  later,  processing  of  the 
data  from  such  a  system  requires  only  a  2  element  or  2  channel  processor  rather  than 
a  n  channel  system  (Sect.  2.6). 

2.5  COG  SIMULATOR 

Analysis  of  the  COC  processor  (Sect.  2.2)  and  the  adjunct  antenna  concept 

(Sect.  2.4)  were  best  facilitated  by  simulation.  The  signals  received  at  one 

element  of  a  phased  array  is  the  sum  of  N  different  noise  sources  at  angles  f-  , 

m 

range  delays  :  ,  frequencies  f  and  bandwidths  B  .  Generation  of  such  composite 
m  m  m 

signals  for  one  or  more  receiving  antenna  elements  is  quite  complex  considering 
the  multitude  of  possible  noise  source  parameters. 

COC  simulator  for  APAR  was  developed  with  all  the  necessary  features  to 
permit  study  of  advanced  phased  array  systems  with  space  and  time  diversity.  The 
currently  operational  simulator  can  handle  20  adaptive  elements  (determined  by 
specifying  the  variable  RN) .  The  number  of  noise  sources,  their  locations  and 
their  frequencies  are  controlled  by  the  simulation  parameters  RNO,  XSOl'R,  and  XFRTO. 
Similarly,  the  antenna  element  spacing  is  t^ken  to  be  1/2,  but  can  easily 
be  altered  within  the  program  (in  practice  this  is  done  most  simply  by  varying  the 
frequency  and  lienee  the  wavelength  of  the  input  signal).  The  power  levels  of  the 
different  noise  sources  are  taken  to  be  unity  unless  otherwise  specified.  We  can 
also  produce  random  uncorrelated  received  signals  as  are  necessary  in  practice  bv 
using  different  phases  for  different  received  sinewaves.  This  can  be  used  to  repro- 
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sent  uncertainty  in  the  synchronization  of  different  noise  sources  as  well  as 
producing  noise  sources  with  different  degrees  of  coherence. 

Let  us  first  consider  the  case  of  sine  wave  signals,  since  these  represent 
the  simplest  type  of  received  signal.  In  this  case,  the  time  delay  t  between 
adjacent  antenna  elements  is  easily  expressed  by  including  a  simple  phase  differ¬ 
ence 

4>  =  2u  ft  (2.8) 

in  the  received  signals.  For  each  frequency,  the  necessary  received  signals  at 
different  array  elements  can  thus  be  described  by  introducing  0  in  (2.8).  In  our 
Fortran  program,  we  vary  the  sampling  rate  of  the  input  signal  to  represent  differ¬ 
ent  receivers.  To  verify  the  fidelity  of  such  monof ri quency  signals,  we  take  their 
FFT  and  observe  that  it  has  a  non  zero  value  at  only  one  frequency. 

To  model  uneorrelated  noise  sources,  we  introduce  a  small  frequency  deviation 
.'.f  from  f  such  that 

T.'.f  =  +n  (2.9) 

where  n  is  an  integer  and  T  is  the  total  signal  duration.  The  associated  Fortran 
program  is  listed  in  Table  2.2. 

In  Fig.  2.8,  we  shown  a  pseudo  isometric  display  of  the  received  signals  at 
the  N  elements  of  the  phased  array  with  time  and  antenna  element  as  independent 
variables  (horizontally  and  vertically).  In  this  example,  11  independent  receiving 
elements  are  assumed  with  a  single  noise  source  of  amplitude  0.7  at  30^*  with  inte¬ 
gration  over  500  signal  points  employed.  As  seen,  different  sine 

wave  noise  sources  align  along  lines  at  angles  corresponding  to  0  of  the  source. 

m 

The  integration  along  these  lines  gives  fine  peak  location  accuracy  as  noted  in 


Chapter  1 . 
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To  describe  wide-band  noise  sources  or  noise  sources  with  n  given  bandwidth, 
our  simulator  employs  LFM  signals,  since  their  center  frequency  and  bandwidth  can 
easily  be  controlled  and  since  their  correlation  is  well  predictable.  Such  signals 
are  most  necessary  in  determining  the  performance  of  a  APAR  system  with  wide-band 
signals  and  wide-L.and  receivers  with  known  deterministic  outputs.  The  Fortran 
program  for  this  is  listed  in  Table  2.3.  It  is  capable  of  generating  all  necessary 
received  signals  from  multiple  sources,  each  characterized  by  a  different  LFM 
sequence  and  each  with  a  different  center  frequency  and  bandwidth.  The  initial 
phase  and  hence  the  coherence  of  the  different  noise  sources  can  be  controlled  by 
randomizing  the  phase  of  the  different  signals.  The  program  parameters  XFREQ, 

XRATE,  AMPL  and  N’Osamples  determine  the  parameters  of  the  LFM  signal.  The 

LFM  waveform  produced  with  XFREQ  =  1.0,  XRATE  =  0.1,  AMPL  =  1.0  and  NO  samples  =  10.23 

is  shown  in  Fig.  2.9  and  is  seen  to  be  of  the  desired  form. 

For  such  a  wide-band  signal,  time  delays  between  corresponding  received  ele¬ 
ments  cannot  be  represented  as  simple  phase  shifts  since  a  given  time  delay  can 
correspond  to  different  phase  delays  at  different  frequencies.  Thus,  for  a  given 
array  with  element  spacing  d,  is  fixed  and  in  our  routine  we  simulate  delayed 
signals  hv  constructing  the  corresponding  signal  and  shifting  it  in  the  time  domain. 

To  more  correctly  model  true  APAR  noise  sources,  random  signal  sequences  with 
controlled  statistical  properties  rather  than  the  deterministic  sine  wave  and  LFM 
signals  considered  thusfar  are  necessary.  To  describe  such  noise  sources,  we 
specify  their  amplitude  probability  density  functions  and  the  average  signal  power 
(variance)  and  we  characterize  their  frequency  response  by  their  power  spectrum.  We 
select  the  Gaussian  PDF  as  an  appropriate  choice  because  it  is  easily  mathematically 
modelled  and  because  it  represents  manv  possible  and  realistic  noise  scenarios. 

Central  limit  theorem  [10)  cons  ider.it i ons  oas  i  1  v  verify  the  appropriateness  of  such 
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a  model.  Various  noise  source  spectra  such  as  white  noise  can  be  generated 
by  various  digital  filters  (separate  routines  exist  for  this  [2))  to  produce  noise 
sources  of  any  desired  bandwidth. 

The  Fortran  program  in  Table  2.4  produces  the  APAR  required  signals  for 
two  Gaussian  noise  sources.  It  uses  the  library  function  RAN  from  the  DEC 
operating  system.  This  routine  generates  random  numbers  uniformly  distributed 
between  0.0  and  1.0  that  are  then  used  to  generate  a  Gaussian  random  variable 
sequence  as  in  [11].  The  program  parameters:  RN  (number  of  receiving  elements), 

RNO  (number  of  noise  sources),  and  XSOUR  (angular  location  of  the  noise  sources) 
are  used  to  control  the  signal  produced.  The  probability  density  of  the  data 
produced  by  this  routine  was  verified  by  generating  10,000  samples  of  Gaussian 
noise  and  inspecting  the  histogram  of  the  samples  to  verify  its  Gaussian  shape. 

With  687.  of  the  samples  found  to  lie  between  +  0.40  out  of  the  full  +  1.0  range 
the  data  also  has  the  desired  standard  deviation  '  =  0.40.  The  fidelity  of  this 
random  data  is  also  verified  by  averaging  the  FFT  of  many  sampled  realizations  to 
decrease  the  estimation  variance.  A  smooth  and  flat  spectrum  resulted  from  this 
experiment.  However,  for  white  random  noise,  observing  that  its  autocorrelation 
is  a  delta  function  is  a  simpLer  test.  To  vary  the  bandwidth  of  such  white  random 
noise,  digital  filters  using  Butterworth  polynomials  and  bilinear  transformations 
arc  used  as  explained  in  |2|. 

Tiles.-  new  routines,  combined  with  those  in  |2)  now  enable  us  to  produce  mono, 

I.FM  and  random  signals  whose  frequency,  bandwidth,  number,  duration,  delay  and 
mutual  coherence  can  he  independently  varied.  They  now  include  multiple  signals 
that  are  adaptive  in  spate  and  time.  These  signals  tvpil'v  the  type  of  data  to  be 
expected  at  the  receiving  elements  of  a  phased  array  and  adjunct  antenna  and  will 
he  tin-  tines  used  in  the  ('<)('  anil  TI  AH  systems  both  tor  simulation  ind  in  the  computer 
driven  elf  t runic  hardware  support  svstem. 
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2 . 6  HYBRID  TIME  AND  SPACE  INTEGRATING  (TST )  APAR  PROCESSOR 

The  basic  AO  TI  correlator  was  shown  in  Fig.  2.1  aid  a  uu 1 t i-channel  system 

in  Fig.  2.2.  They  are  described  in  Sect.  2.2.  In  this  section,  we  consider  a 

new  AO  processor  that  can  provide  both  the  angular  and  temporal  noise  distribution 

X  (•'  ,  f  I  of  the  antenna '  s far-f ield  noise  pattern.  This  is  necessary  when  wide- 
m  \  m  m  / 

band  noise  sources  and/or  wide-band  receivers  are  employed  and  is  necessary  to 
produce  a  APAR  system  with  adaptivity  in  both  angle  and  time  (i.e.  space  and  fre¬ 
quency)  . 

We  concentrate  on  the  adjunct  antenna  system  (Sect.  2.4)  for  the  reasons 
described  earlier,  i.e.  it  allows  the  necessary  target  angles  to  be  resolved  by 
correlation.  Only  a  2  channel  processor  (e.g.  Fig.  2.1)  is  thus  necessary  with 
one  signal  being  the  received  signal  from  the  adjunct  antenna  and  the  second 
signal  being  the  received  signal  at  the  central  reference  element  of  the  main  X 
element  phased  array.  If  a  multi-channel  system  were  used  (e.g.  Fig.  2.2)  in 
which  the  N  phased  array  received  signals  were  correlated  with  the  reference 
adjunct  antenna  element,  it  would  not  be  possible  to  resolve  the  locations  of  t  lie 
correlation  peaks  on  the  different  output  channels  as  noted  in  Sect.  2.4.  Since 
no  added  resolvable  information  is  obtained  from  the  different  channels,  only  two 
are  used. 

We  now  consider  how  to  use  the  TI  system  of  Fig.  2.1  in  the  adjunct  antenna 
phased  array  scenario  and  how  to  modify  it  to  produce  both  angle  and  frequency 
information  on  t he  ant cnna ' s  far-field  noise  pattern.  The  case  of  a  single  wide-hand 
noise  source  simultaneously  generating  multiple  discrete  frequencies  is  considered 
first  because  of  the  notutional  and  conceptual  simplicity  it  provides.  The  APAR 


problem  is  thus  effectively  reduced  to  a  2  antenna  element  case  (with  a  large 
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d  '/-  spacing  between  antenna  elements).  The  two  received  signals  are  described 

by 

s  (t)  =  s(t-r/2)  (2.10a) 

a 

s,(t)  =  s(t+r/2),  (2.10b) 

b 

where  the  time  delay  (in  seconds)  between  the  two  received  signals  is 

!  =  (d/c)  sin  f‘.  (2.11) 

Note  that  c  is  the  velocity  of  propagation  of  the  radiation  and  that  t  is  indepen¬ 
dent  of  the  frequency  of  the  radiation  and  hence  the  frequency  distribution  of  the 
noise  sources  and  depends  only  on  the  angle  ••  at  which  the  noise  sources  are  located. 

The  basic  concept  used  in  the  TSI  processor  is  to  produce  frequency  filtered 
versions  of  the  received  signal  s(t)  with  different  bandpass  center  frequencies 
and  to  correlate  these  with  the  original  signal  s(t).  Since  the  correlation  peak 
value  is  proportional  to  the  energy  of  the  signal,  we  can  thus  determine  the  amount 
of  signal  energy  in  different  frequency  bands.  [Although  similar  information  can 
be  obtained  from  a  simpLe  Fourier  transform  of  s(t),  this  will  not  work  when  s(t) 
contains  multiple  signals  at  different  angles,  each  with  a  different  bandwidth  and 
I requenev  distribution].  This  concept  is  shown  in  block  diagram  form  in  Fig.  2.10 
where  N  +  1  correlations  are  produced:  N  correlations  R  -R  of  the  narrow-band 
filtered  signals  s^-s  (with  band  pass  filters  h^-h^)  with  the  wide-band  reference 
signal  s(t-'/2)  and  the  correlation  R^,  of  the  wide-band  signals  and  s^ . 
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Figure  2.10  Simplified  block  diagram  of  a  two-channel  adjunct  antenna  TSI 
processor . 

The  peak  values  of  R^-R  denote  the  energy  contained  in  the  N  frequency  bands 
chosen  of  the  wide-band  signal  s.  If  we  increase  the  number  of  filters  N,  we  in¬ 
crease  the  frequency  resolution.  However,  increasing  N,  decreases  the  width  of 
the  h^  filters,  thus  producing  weaker  and  .much  wider  output  correlations,  hence 
making  correlation  plane  detection  more  difficult.  A  compromise  in  N  (large  N  for 
resolution  and  small  N  for  accurate  detection)  is  necessary.  The  choice  of  N  deperi  s 
on  the  signal  bandwidth,  the  detector  sensitivity,  the  integration  time  and  the 
number  of  adaptive  weights  available.  Since  there  is  no  need  to  estimate  the  fre¬ 
quency  and  angular  resolution  of  the  far-field  noise  to  a  resolution  better  than 
that  for  which  the  radar  can  cancel  and  adapt;  the  number  of  adaptive  weights  avail¬ 
able  also  determine  the  choice  of  N. 


The  output  R^  from  the  top  channel  is  the  wide-band  correlation  of  the  wide¬ 
band  signals  s^  and  s^ .  This  output  will  have  a  much  narrower  correlation  peak 
width  than  the  other  N  correlations  and  is  thus  used  to  determine  the  target's 
angle  whereas  the  other  N  correlations  provide  information  on  the  frequency 


distribution  of  the  noise  at  different  angles. 
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To  realize  Che  processing  depicted  in  Fig.  2.  10, an  AO  TI  correlator  is  used. 

The  block  diagram  of  such  a  system  is  shown  in  Fig.  2.11.  In  this  system,  shifted 

versions  of  s,  are  produced  bv  a  AO  cell.  Different  band  pass  filtered  versions 
b 

of  s,  are  produced  bv  frequency  plane  filtering  (BP  filts)  and  correlation  is 
b  " 

achieved  by  multiplying  (X)  the  shifted  signals  and  time  integrating  (Tl)  the  prod¬ 
ucts  on  an  output  detector. 

A  schematic  diagram  of  the  AO  hybrid  TSI  correlator  is  shown  in  Fig,  2.12, 

This  AO  architecture  combines  the  best  features  of  the  SI  and  the  TI  AO  correlators 
into  a  new  hybrid  architecture.  This  system  uses  only  one  real-time  1-D  trans¬ 
ducer  (AO  cell),  a  time  sequentially  modulated  LED  light  source  (or  a  point  AO 
modulator)  and  a  fixed  mask.  The  received  signal  s^lt)  is  used  to  time  sequentially 
modulate  the  output  from  a  LED  light  source,  whose  output  is  expanded  to  uniformly 

illuminate  an  AO  cell  at  P„.  The  AO  cell  is  fed  with  the  signal  s,  (t).  With  s 

Z  D  3 

and  s,  as  in  (2.10),  the  light  distribution  leaving  P  is 
b  Z 

u„(x,y,t)  =  s(t-' /2)s(t+7 /2-x) ' (y) ,  (2.12) 

where  the  AO  cell  produces  the  continuous  shift  x  and  the  ' (y )  function  notes  that 
this  is  the  pattern  at  v  -  0. 
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At  P .j,  we  see  that  each  vertical  line  passes  a  different  spatial  frequency 
part  of  u,^  due  to  the  apertures  located  at  different  horizontal  positions  on 
each  line  of  the  mask.  We  describe  the  light  leaving  the  n-th  line  at  P  by 

Ul(  x,n’ t  )  =  u3RECT(-x-n)*  (2-14) 

where  the  R1XT  function  describes  the  band  pass  filter  at  P„  and  .  is  the  center 

3  n 

frequency  of  the  n-th  channel  of  this  hand  pass  filter  mask. 


I.enses  Ln  image  P_  onto  P  in  the  vertical  direction  and  produce  the  Fourier 
2  3  s 

transform  of  u^  in  the  horizontal  direction  along  The  light  distribution 

incident  on  P.  is 

■4 

U4  ('  x  ’ "  ’  t)  =  s(t“  /2)sn(t+  /2-’x)’  U-,5) 

where  s  denotes  the  n-th  narrow-band  filtered  version  of  s(t)  and  '  is  the  spatial 
n  x 

shift  variable.  The  detector  at  P,  time  integrates  the  light  distribution  in  (2.15) 

4 

for  T  seconds  and  produces  the  correlation  output 

U,T)  0  *<t-  /2).n(t+i/2-Tx)dt 

.  s<o  ©»„<■>- »„('«)• 

The  top  row  of  the  filter  mask  at  P  is  transparent  and  thus  along  the  top 

line  in  P,  we  find  the  wide- band  correlation  of  the  wide-hand  unfiltered  signals 

s  and  s, 
a  b 

R 

Thus,  as  before,  R  provides  narrow  correlation  peak  widths  and  thus  good  time 

deluv  =  or  angle  resolution  ■  of  the  target,  whereas  the  R  correlations 
x  m  m  n 

contain  the  necessary  frequency  distribution  information  on  the  target.  We  discuss 
these  systems  later  in  Sect.  2.10  and  for  now  note  that  the  narrow-hand  and  wide¬ 
band  sional  correlations  must  be  considered  in  an  APAR  adaptive  in  both  space  and 


(1/T)  .  s(t-r/2) 

0 


s|t+T/2--ixjdt. 


(2.17) 


t  i  nu* . 
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2.7  ELECTRONIC  SUPPORT  SYSTEM 

In  this  section,  we  describe  the  two  electronic  support  systems  being  assembled 
to  provide  the  necessary  drive  and  synchronization  for  the  COC  APAR  correlators.  We 
refer  to  the  two  systems  as  computer-driven  and  hardware  electronic-support  systems 
respectively.  Since  the  received  signal  at  one  element  of  a  phased  array  is  the  sum 
of  the  signals  from  N  different  noise  sources  each  of  which  can  be  at  a  different 
angle,  range,  strength,  center  frequency  and  bandwidth,  a  single  received  signal  is 
quite  complex.  Since  each  of  its  constituent  component  subsignals  (from  the  differ¬ 
ent  individual  noise  sources)  changes  with  the  noise  source  scenario  and  with  the 
antenna  spacings,  producing  the  different  received  signals  for  each  case  is  in  it¬ 
self  a  major  effort.  We  intend  to  use  a  combination  of  simulation  and  experimentation 
(with  the  two  electronic  support  systems)  to  conduct  our  APAR  study.  In  this  section, 
we  consider  the  different  electronic  support  svstems,  the  advantages  and  disadvantages 
of  each,  the  design  and  performance  of  each.  An  overview  of  both  techniques  is 
given  to  convey  the  key  points  of  the  diverse  signal  generation  concepts. 

In  the  computer-driven  system,  the  two  received  signals  are  digitally  computed 
using  the  simulator  (Sect.  2.5)  and  their  time  histories  are  produced,  stored  in  a 
digital  file  and  recorded  on  tape.  This  magnetic  tape  is  then  brought  from  the  main 
computer  facility  to  our  own  dedicated  PDP-11  where  it  is  loaded  into  our  tape  file 
and  then  placed  in  the  large  2  M-bit  on-line  memory  peripheral  available  on  our  sys¬ 
tem.  The  signals  are  then  read  from  this  memorv  as  base  band  data,  D/A  converted, 
up-converted  to  the  center  frequency  of  the  AO  cells  by  heterodyning  with  an  oscilla¬ 
tor  using  two  mixers.  Bias  is  then  added  through  T  splitters,  the  signals  arc  then 
attenuated  (ATT)  to  bring  them  to  the  proper  range  for  the  final  power  amplifier 
(AMP)  to  drive  the  point  source  (LED  or  AO  modulator)  and  the  AO  line  itself.  This 
svstpm  is  shown  in  block  diagram  form  in  Fig.  2.15. 
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Figure  2.13  Block  diagram  of  the  computer  driven  electronic  support  system. 

In  the  hardware  system  (Fig.  2.14),  separate  noise  generators  (XO)  are  used 
(one  for  each  noise  source).  The  bandwidth  and  amplitude  of  each  can  be  controlled 
directly  on  die  generators.  Delay  lines  are  used  to  delay  and  band  pass  filter  each 
signal.  Separate  oscillators  enable  the  center  frequency  of  each  noise  source  to 
be  independently  controlled.  The  signals  from  the  different  noise  sources  with  the 
proper  delays  and  bandwidths  and  center  frequencies  are  then  combined  and  fed  to 
the  standard  T,  ATT,  AMP  circuits  to  provide  input  signals  of  the  necessary  power 
and  linearity  for  the  point  modulator  on  the  AO  lines. 

From  tlie  above  two  paragraphs,  the  basic  advantages  and  disadvantages  of  the 
two  techniques  should  be  apparent.  The  computer-driven  system  is  capable  of  pro¬ 
ducing  essentially  any  type  of  signal  and  is  thus  much  more  flexible  than  the  hard¬ 
ware  system.  It  also  offers  a  larger  range  of  delays  than  one  can  achieve  with 


hardware  delay  lines,  probably  better  synchronization,  and  moreso  the  ability  to 
accurately  repeat  long  (32  msec)  duration  signals.  This  system  is  more  complex 
in  the  sense  that  it  requires  operation  of  the  entire  PDP-11  system  and  its 
peripherals,  but  moreso  the  dynamic  range  and  bandwidth  of  the  data  it  can  pro¬ 
duce  are  severely  limited  by  the  speed  and  size  of  the  on-line  memory.  Conversely, 
the  hardware  support  system  is  severely  limited  in  the  flexibility  of  the  signals 
and  delays  it  can  produce,  but  the  svstem  is  more  directly  controlable  and  can 
produce  larger  dynamic  range  data.  The  bandwidth  and  duration  of  the  signals 
can  be  longer  in  the  hardware  system,  but  random  signals  cannot  be  accurately 
repeated  on  this  system.  Inclusion  of  LFM,  PRN  and  other  deterministic  signal 
sources  in  the  system  is  possible  ant!  would  improve  its  performance  at  the  cost 
of  addition  expense  and  complexity.  Both  signal  generation  electronic  support 
systems  are  thus  seen  to  be  mutually  compatible  and  of  use  for  different  cases. 

Let  us  now  first  quantify  the  performance  of  the  computer-driven  system.  The 

on-line  memory  is  the  major  limitation  in  this  system.  It  is  a  2-M-bit  memory, 

normally  configured  as  512  x  512  x  8-bits.  Its  entire  contents  can  be  randomly 

accessed  in  30  msec.  Serial  mode  read  out  is  normally  used  with  one  8-bit  word 

accessed  about  every  100  nsec.  Thus,  the  minimum  delay  increment  possible  is 

',t  .  =  100  nsec  and  delnv  increments  in  multiples  of  this  are  possible  (there  is 

min 

no  point  in  delays  greater  than  0.5  times  the  duration  of  the  AO  cell,  or  in  our 
case  25  i.sec) .  Similarly,  signal  durations  up  to  30  msec  (and  approximately  250,000 
samples)  are  possible  and  they  can  be  accurately  repeated  by  recycling  the  memory. 

The  bandwidth  and  dynamic  range  of  the  data  from  this  system  are  its  major 
limitations.  With  one  8-bit  output  word  every  100  nsec,  the  signal  bandwidth 
possible  is  10  MHz  (base  band)  or  20  MHz  USB.  Two  signals  are  necessary  and  s  . 
If  these  are  real,  then  A-bit  words  (16  gray  levels)  can  be  used  for  o;u  h ,  with 
four  of  the  eight  ouLput  bits  at  each  sample  corresponding  to  one  A-bit  sample  ol 
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oath  signal.  This  is  the  situation  shown  by  the  two  memory  ports  in  Fig.  2.14. 

If  complex-valued  data  is  to  be  processed,  then  two  signals  (real  and  imaginary 
are  necessary  for  each  waveform  and  four  output  signals  must  be  obtained  from 
each  8-bit  output  word  from  the  memory  (i.e.  2-bit  or  four  level  accuracy  is 
possible).  If  less  bandwidth  is  acceptable,  256  level  or  8-bit  sampies  can  be 
obtained  at  5  MHz  by  assigning  each  alternate  output  word  to  be  a  sample  of  a 
different  signal. 

These  specifications  are  quite  compatible  with  the  available  AO  cell  and  the 
intended  COC  TS I  APAR  processor .  Moreso,  a  separate  analysis  has  been  performed  by 
us  that  indicates  that  2-4  bit  quantization  of  random  signal  data  does  not  apprecia 
blv  effect  the  SNR  of  the  resultant  correlation,  even  when  reasonable  input  noise 
is  present  [  1 2  J .  Thus,  this  computer-driven  system  is  quite  attractive  for  the 
intended  application.  It  is  clearly  the  most  flexible  system  capable  of  providing 
the  widest  range  of  waveforms  and  delays. 

It  is  possible  to  extend  botli  electronic  support  systems  to  provide  the  neces- 
sarv  input  signals  for  complex-correlations.  To  achieve  this,  four  signals  are 
necessary  from  the  memory,  each  pair  is  quadrature  modulated  and  fed  to  the  point 
modulator  and  the  AO  cell  respectively.  Although  complex  correlations  yield  3  dB 
more  SNR  than  do  real  correlations,  in  many  adaptive  radars,  the  phase  of  the  input 
signal  might  not  be  preserved  with  heterodyning  and  thus  complex  correlation  would 
not  be  of  use.  In  practice,  one  can  achieve  complex  correlation  by  inputing  differ 
ent  signals  and  cycling  the  system  twice  and  adding  the  separate  outputs  electroni¬ 
cal  lvwitli  quadrature  detection.  Thus,  complex  correlations  can  be  performed  on  the 
existing  system  with  added  complexity.  Since  the  central  point  and  purpose 
of  the  D)C  correlator  can  be  proven  without  complex  correlations,  we  have  elected 
to  perform  real  rather  than  complex  operations  on  this  system. 
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The  computer-driven  electronic  support  system  is  nearly  operational.  Addi¬ 
tional  funds  are  necessary  for  a  new  operating  disc  system  and  for  redesign  of  the 
sync  board  to  suppress  horizontal  and  vertical  blanking  pulses  during  fly  back. 

These  are  expected  within  several  months.  Increasing  the  bandwidth  of  the  memory 
to  20  MHz  is  possible  but  would  require  extensive  rewiring  (twisted  pairs)  and 
additional  high  bandwidth  clock,  driver  and  receiver  circuits  plus  higher  frequency 
n/A  converters.  These  improvements  are  possible  but  again  are  not  necessary  to 
prove  the  feasibility  of  the  COC  concept.  All  of  the  associated  mixers,  T  connectors 
oscillators  and  amplifiers  have  been  unified  and  are  common  to  both  electronic 
support  systems. 

For  the  above  reasons,  for  near-term  experiments  on  the  COC  APAR  radar  pro¬ 
cessors  we  will  use  the  hard  wired  electronic  support  system  of  Fig.  2.14.  The 
use  of  real  not  complex  correlations  is  quite  attractive  for  the  hardware  system 
since  absolute  synchronization  of  the  real  and  imaginary  parts  of  the  signal  are 


not  needed  then.  The  range  of  the  delays  possible  in  the  hardware  system  is  less 
than  in  the  computer-driven  one  (400  nsec  minimum),  but  this  is  still  quite  accept¬ 
able  for  APAk  applications.  By  varying  the  delay  line  taps,  incremental  delays 
in  units  of  400  nsec  can  be  realized.  At  400  nsec,  the  clock  to  the  delay  line 
permits  a  maximum  signal  bandwidth  of  2.5  MHz.  It  is  possible  to  continuously 
vary  the  delay  but  as  the  delay  is  linearly  increased,  the  bandwidth  is  linearly 
decreased.  This  clocking  feature  of  the  delay  lines  will  be  used  to  vary  the 
bandwidth  of  the  different  noise  sources  in  our  electronic  system. 

In  the  system  of  Fig.  2.14,  the  frequency  of  the  different  noise  sources  (maximum 
of  three)  is  set  by  separate  oscillators.  The  bandwidth  of  each  source  will  he  con¬ 
trolled  by  varying  the  clock  frequency  to  the  delay  lines.  The  delays  will  be  ad¬ 
justed  with  different  tap  weights  on  the  delay  lines.  A  low  pass  filter  follows 
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each  delay  line  to  remove  aliasing  and  to  reconstruct  the  continuous  bandwidth 
signal  from  the  discrete  one.  This  filter  is  chosen  to  cover  the  necessary 
source  bandwidth  range  anticipated.  The  strength  of  each  noise  source  is  adjusted 
directly  on  the  noise  generator.  A  combiner  (microwave  splitter)  is  necessary 
rather  than  an  operational  amplifier  to  combine  the  signals  because  of  the  high 
frequencies  involved. 

The  dynamic  range  of  the  electronic  system  is  excellent,  60  dB,  limited  by  the 
delay  lines,  compared  to  23.5  dB  for  the  two-channel  4-bit  computer  system.  The 
oscillator  used  (Techtronics  5  G  503)  operates  from  250  KHz-250  MHz.  The  mixer 
is  the  Mini-Circuits  model  ZAY-3  with  frequency  response  from  70  KHz-200  MHz  and 
linear  operation  up  to  15  dBm  signal  input.  The  Mini-Circuits  model  ZAY-3  atten¬ 
uator  has  a  frequency  range  from  1-200  MHz  and  continuously  controlable  attenua¬ 
tion.  The  broad  band  amplifiers  are  the  EMI  model  300  L  with  an  operating  range 
of  250  KHz-110  MHz.  Their  maximum  RF  output  (3  W)  is  much  larger  than  the  0.1  W 
needed  for  linear  operation  of  the  AO  cells,  thus  their  linearity  at  the  operating 
power  levels  used  is  much  better  than  the  2. 52  obtained  at  3  W. 

The  hardware  electronic  support  system  is  fully  operational.  The  computer- 
driven  system  represents  a  new  feature  emphasizing  flexibility  and  multi-purpose 
operation  of  optical  signal  processing  systems.  Both  systems  are  necessary  to  fully 
pursue  a  COC  approach  to  APAR  radar.  Use  of  these  and  other  support  hardware  in 
simulation  systems  are  included  in  Sects.  2.9  and  2.10. 

2 . 8  POST  PROCESSING 

i 

i 

In  this  section,  we  discuss  various  post-processing  techniques  by  which  the 
J  adaptive  weights  for  the  APAR  can  be  obtained  from  data  in  other  forms.  i'he  I0P 

t 


(Chapter  3)  and  the  WDP  (Chapter  4)  techniques  operate  on  the  covariance  matrix  M 
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of  the  antenna  and  thus  can  also  be  thought  of  as  optical  post-processing,  whereas 
calculation  of  M  (by  digital  or  other  techniques)  is  the  primary  preprocessing  in  this 
case.  We  refer  to  this  method  as  sampled  matrix  inversion.  Obtaining  an  accurate 
estimate  of  M  must  be  given  some  attention.  In  (13],  K  >  2  N  samples  (for  an  N 
element  array)  were  found  to  give  an  average  loss  ratio  better  than  3  dB  in  the 
sampled  estimate  of  M.  However,  more  recent  work  [14]  has  pointed  out  that  this 
implies  that  one  half  of  the  time  the  loss  will  be  above  3  dB  and  half  of  the  time 
it  will  be  below  3  dB.  Thus,  recent  work  suggests  that  K  ->  4  N  be  used  and  shows 
that  this  insures  that  the  probability  of  the  loss  being  above  3  dB  is  then  only 
0.0032.  More  attention  to  this  issue  of  accurate  estimation  of  the  sampled  co- 
variance  matrix  is  necessary  and  should  be  pursued  in  a  later  phase  of  this  work. 

Our  major  concern  in  post-processing  arises  in  conjunction  with  the  COC  system. 

In  this  system,  the  location  (in  angle  or  time  delay  and  frequency)  of  the  noise 
sources  is  calculated  in  the  optical  system  and  the  adaptive  weights  W  must  still 
he  obtained  from  these  data  in  a  digital  post-processor.  The  weights  are  chosen 
to  produce  nulls  in  the  antenna  pattern  at  noise  locations  while  maintaining  response 
in  the  signal  direction,  with  no  stipulations  on  the  antenna  response  at  other 
locat ions . 


From  the  preliminary  discussion  in  Sect.  2.6,  we  saw  that  the  cases  of  narrow- 

band  and  wide-band  noise  and  noise  at  the  same  angle  and  different  frequencies  and 

noise  sources  at  the  same  frequency  but  different  angles  must  be  treated  separately 

with  different  post-processing  necessary  (as  described  below)  and  with  different 

slit  filter  widths  h  in  Fig.  10  chosen  for  each  case.  In  Sect.  2.10,  we  will  unify 
n 

these  different  cases.  For  now  we  first  consider  the  case  of  a  narrow-band  noise 


source  at  an  angle  A  and  at  a  frequency  f  .  A  simple  post-processor  can  be  used 
m  '  m 

in  this  case  as  we  now  describe. 
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The  time  delay  i  in  units  of  seconds  (the  quantity  measured  in  the  COC  processor) 
is 

;  =  (d/c)  sin  n  (2.18a) 

m 

and  is  thus  independent  of  the  frequency  f  of  the  noise  source.  However,  the  time 
delay  in  radians 

’  =  1 2  ■  f  \  (d/c)  sin  P  (2.18b) 

l  ml  m 

depends  upon  the  frequency  f  of  the  noise  source.  Moreso,  when  the  adaptive  weights 

m 

w^  are  computed,  then  the  resultant  antenna  response  E(n)  has  a  simple  Fourier  trans¬ 
form  relationship 

E(0)  =  1  w  exp/j2iinf  sin  P  ]  .  (2.19) 

n  n  V  m  m' 

Thus,  a  simple  digital  FFT  post-processor  can  be  used.  However,  note  that 

the  frequency  f  of  the  noise  source  must  be  known  in  (2.19).  Moreso,  if  the  noise 
m 

source.,  exist  at  different  frequencies,  then  a  different  FFT  is  necessary  to  compute 
the  weights  w^  to  place  nulls  at  each  f‘mj  parameter  pair.  If  this  is  done  sep¬ 

arately  for  each  noise  source,  no  convenient  way  to  weight  the  sum  of  the  individual 

weights  w  to  null  each  noise  source  has  vet  been  found.  Since  one  set  of  weights 
n 

may  place  a  null  in  the  desired  location,  whereas  another  set  of  weights  may  place 
a  peak  there,  the  resultant  response  in  the  desired  location  will  be  non-zero. 

Thus,  we  first  consider  the  problem  of  the  pos t-process  ing  required  when  the 
noise  sources  lie  at  the  same  frequency  and  are  narrow-band.  We  then  consider  the 
case  of  wide-band  noise  sources  separately  later.  Each  case  will  be  shown  to  re- 
ati ; re  a  different  post-processing  algorithm.  Solutions  for  each  of  the  different 
no :  ;e  pattern  cases  will  be  shown  to  exist  (here  and  in  summary  in  Sect.  2.10). 

Thus,  all  AFAR  cases  can  be  solved.  From  the  form  of  the  output  data  at 


of  Fig. 
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2.12,  Che  post-processing  necessary  can  be  determined  from  the  noise  field  present 
at  a  given  time.  Thus,  the  resultant  APAR  processor  appears  to  he  applicable  to 
a  variety  of  different  cases. 

Let  us  first  consider  the  narrow-band  noise  case,  when  all  noise  sources  are 
at  the  same  frequency.  From  (2.18b),  we  see  that  noise  sources  at  different  fre¬ 
quencies  f  and  angles  0  can  result  in  the  same  time  delav  '  in  units  oi  radians, 
m  m  m 

This  implies  that  if  we  chose  an  arbitrary  frequency  f  at  which  we  assume  the  noise 
source  to  exist,  the  angle  ''  computed  for  the  noise  source  will  satisfy  the  equation 

f.  sin  A  =  f  sin  <•  ,  (2.20) 

0  mm 

where  f  and  0  represent  the  correct  frequencv  and  angle  of  the  noise  source.  This 
mm 

result  in  (2.20)  is  important  because  it  implies  that  we  can  choose  an  arbitrary 
frequency  f^  at  which  to  assume  the  noise  source  exists,  calculate  the  weights 
to  produce  a  null  at  the  angle  ('■  and  that  the  resultant  set  of  weights  will  auto¬ 
matically  produce  a  null  at  the  correct  frequency  f  and  angle  f*  as  in  (2.20). 

m  in 

We  refer  to  this  adaptive  array  post-processing  algorithm  as  the  self-compensating 

or  self-correcting  algorithm.  It  is  only  of  use  with  narrow-band  sources  and  sources 

at  the  same  frequency  or  at  only  a  few  different  frequencies.  To  demonstrate  the  use 

of  this  algorithm,  we  simulated  the  multi-channel  COC  system  for  the  case  of  a  noise 

source  at  angle  >:  and  frequencv  f  .  The  multi-channel  output  from  the  COC  system 
m  '  m 

was  obtained  and  slit  integration  of  the  pattern  was  performed  as  described  in  [2], 

The  resultant  output  shown  in  Fig.  2.15a  exhibits  a  peak  at  0,  corresponding  to  a 
time  delay  ;  calculated  from  (2.18a)  that  is  independent  of  the  frequency  of  the 


noise  source. 


The  simulator  (section  2.5)  was  then  used  to  compute  w  for  this  N 

n  m 

pattern  and  to  then  calculate  E(()  from  (2.19)  using  an  FFT  routine.  For  an 

assumed  frequency  f^,  the  resultant  E(i)  pattern  of  Fig.  2.15b  was  obtained. 

It  exhibits  a  null  at  c  =  >*  as  is  desired.  We  then  repeated  the  simulations 

for  a  noise  source  frequency  f  and  obtained  the  same  set  of  weights  but  the 

m 

different  E(i')  pattern  shown  in  Fig.  2.15c..  As  seen  this  pattern  exhibits  a 

null  at  the  desired  P-  location  rather  than  at  •  . 
m 

From  this  discussion  and  simulation  example,  we  see  that  the  self-correcting 
post-processing  technique  works.  We  also  see  other  features  and  uses  of  the 
simulator  for  calculating  the  adaptive  weights  and  the  resultant  far-field  noise 
distribut ion. 

We  now  consider  a  new  algorithm  by  which  to  calculate  the  set  of  adaptive 
weights  for  the  case  of  many  noise  sources  at  the  same  frequency  but  different 
angular  locations.  This  technique  uses  a  new  version  of  a  projection  mapping  (PM) 
technique  [lb]  that  has  been  used  thusfar  only  for  image  restoration  [17].  We 
apply  this  technique  to  the  problem  of  calculating  the  adaptive  weights  for  a 
AFAR  and  we  suggest  a  modified  algorithm  that  results  in  much  faster  convergence 
of  this  routine.  Such  a  technique  is  necessary  to  make  the  speed  of  the  post-pro¬ 
cessor  compatible  with  the  high  speed  of  the  COC  system  that  estimates  the  noise 
field.  A  digital  post-processor  solution  is  envisioned,  although  modifications 
to  the  [OP  (Chapter  3)  or  the  WDP  (Chapter  4)  system  may  make  an  optical  version 
of  this  algorithm  possible. 

To  describe  the  new  PM  algorithm,  we  reformulate  the  APAR  problem  as  below. 

We  consider  an  array  of  N  receiving  elements  at  vector  locations  in  3-11 

space.  We  denote  the  unit  vector  in  the  direction  of  the  signal  by  ij?.  We  de¬ 
scribe  the  noise  field  by  M  unit  vectors  u.  (where  i  =  1  .  .  .  M)  in  the  direction 

—l 
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of  the  M  noise  sources.  Our  problem  is  to  compute  the  adaptive  weights  W  such 
that  we  maintain  full  response  in  the  direction  of  the  signal  u  ,  i.e. 

C1  W  =  1  ,  (2.21) 

and  null  the  response  in  the  ii  directions,  i.e. 

C.‘  »  =  0,  (2.22) 

T  t 

where  i  =  1 .  .  .  M  and  where  the  vectors  C  and  C.  describe  the  sienal  and  noise 

-  _ x  ^ 

response  as 

CJ  =  exp  j2r|fo/c]|p1.usQ  .  .  .  exp [j  2  ( f Q / c)  ^ J  (2.23) 


■r  r  -  / 

\  /  \~ 

i  ( ~  i  \ 

C.  =  /exp  j 2 ’  p 

ro/cH£i*Hi)_ 

...  exp  j  2t ( fg/c  1  j 

h-i).: 

(2.24) 

Although  the  present  formulation  assumes  all  narrow-band  noise  signals  to  be  at  the 
same  center  frequency  f^,  the  proposed  solution  can  later  be  extended  to  the  case 
of  non-equal  noise  source  frequencies  using  techniques  in  [19-20]. 


(2.27) 
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-en  M  >  N  -  1 ,  no  solution  exists.  In  most  cases  of  interest,  the  number 
of  noise  sources  M  will  be  less  than  the  number  of  adaptive  elements  N  or  M  N  -  1 
(or  in  practice  only  the  M  largest  noise  source  will  be  used).  In  this  case,  the 
pseudo-inverse  method  [15]  can  be  used  and  the  solution  is 

W  =  AT  A  _1  AT  U  (2.28) 

where  A  is  now  an  (M  +  1)  x  N  matrix. 

Computing  the  weights  can  thus  be  reduced  to  the  common  problem  of  solving 
linear  equations  for  which  many  solutions  exist.  However,  most  are  computa t iona 1 Iv 
complex,  especially  for  complex  array  geometries. 

In  the  modified  PM  technique  we  propose,  each  of  the  M  +  1  constraints  in 

(2.22)  is  treated  as  a  hyperplane  in  an  N  dimensional  space  whose  axes  are  the  weights 

w^  to  be  found.  The  problem  is  to  find  the  intersection  of  all  of  these  hvper- 

planes.  The  projection  method  begins  with  an  initial  guess  W^  for  the  intersection 

point.  We  then  project  this  point  onto  the  first  hvperplane  and  obtain  a  new  point 

which  we  then  project  onto  the  second  hyperplane  to  obtain  W^,  etc.  until  the 

vector  point  ^  is  obtained  after  projection  onto  the  last  M  +  1  hvperplane. 

This  completes  the  first  iteration  cvcle.  W_.  ,  ,  is  then  orthogonally  projected 

- —  — ■ — —  —  ■  — M  +1  1 

back  onto  the  first  hvperplane,  giving  +  ^  and  the  process  is  continued  until 
W_„  .  _  is  obtained  at  the  end  of  the  second  iteration. 

The  procedure  is  repeated  until  convergence  occurs.  The  speed  of  convergence 
is  improved  if  the  hyperplanes  are  orthogonal  and  the  larger  the  angle  between  ad¬ 
jacent  hyperplanes  (closer  to  90^)  the  faster  the  convergence  [18].  To  speed  up 

convergence,  we  suggest  a  pair  wise  orthogonal izat ion  technique  that  makes  adjacent 
hvperplanes  orthogonal.  The  method  proposed  is  similar  to  one  used  in  image  restora¬ 
tion  | 17).  We  wish  to  implement  both  the  original  and  the  fast  I'M  techniques  and  to 
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npplv  them  to  the  AI'AR  optimal  weight  problem.  The  modified  PM  technique  should 
produce  deeper  nulls  at  the  correct  locations  in  less  iterations  than  did  the 
original  technique. 


The  conventional  PM  technique  requires  MN  initial  complex  operations  and 

2 

M“  complex  operations  per  iteration.  Conversely,  the  modified  PM  technique  re¬ 
quires  only  2  MM  initial  operations  and  MM  complex  operations  per  iteration.  Since 
M  ■  N  in  general,  a  large  computatioral  savings  is  obtained  by  the  new  second  tech¬ 
nique.  We  thus  advocate  such  a  digital  post-processing  method  for  computation  of 
the  adaptive  weights  in  tue  COC  system.  When  the  frequencies  of  the  N  sources  differ 
and  when  wide-hand  jammers  arc-  present,  the  technique  in  [19-20]  should  be  used  and 
incorporated  into  the  modified  PM  technique.  In  these  advanced  methods,  narrow- 
band  noise  sources  can  be  assumed  and  the  algorithm  modified  to  produce  wider  nulls 
at  the  central  frequencies  of  the  noise  sources.  We  hope  to  pursue  such  techniques 
in  later  phases  of  this  program.  For  now  we  note  that  the  modified  PM  technique 
can  be  used  for  multiple  narrow-band  sources  and  the  advanced  PM  technique  for  wide¬ 
band  noise  sources.  We  also  note  that  both  techniques  operate  directly  on  the  data 
obtainable  from  the  C.OC  output  system  in  the  form  of  the  TSI  processor  of  Fig.  2.12. 

2 . 9  T 1  COC  EXPERIMENTS 

To  obtain  initial  experimental  results  on  the  AO  system  of  Fig.  2.1  in  a  new 
configuration  and  application,  we  considered  the  problem  of  the  accuracy  of  an  opti¬ 
cal  processing  system  for  APAR.  Oik  approach  to  this  problem  is  the  use  of  residue 
arithmetic  [2,  21,  22],  In  such  a  system,  the  input  data  to  he  operated  upon  is 
presented  as  its  residue  values  in  a  chosen  basis  set.  Operation  on  such  data  is 
attractive  for  optical  systems  because  the  dynamic  range  of  the  bit  information  is 
restricted  to  the  corresponding  residue  value  and  moreso  because  the  associated  data 


is  in  the  format  in  which  additions,  multiplications,  anil  subtractions  are  possible 
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without  carries.  The  major  problem  such  a  data  representation  is  how  to  convert 
input  information  into  a  residue  number  system  representation  and  the  associated 
formatting  required. 

To  achieve  such  operation  at  high  data  rates,  we  have  devised  a  version  of  the 

TI  AO  system  of  Fig.  2.1  that  produces  the  associated  conversion  needed  in  a  time 

pulse  position  coding  scheme  rather  than  the  conventional  spatial  pulse  position 

coding  technique  [21].  Tn  residue  arithmetic,  we  represent  an  integer  .J  by  the  n- 

tuple  set  of  remainder  residues  R  .  witli  respect  to  the  N  integer  moduli  m  . 

mi  i 

The  system  we  used  to  demonstrate  the  above  concept  is  shown  in  Fig.  2.16  in 
schematic  form.  In  this  system,  the  temporal  source  modulation  is 

s0(t)  =  -(t-JAt),  (2.29) 

where  we  represent  the  value  of  the  data  by  the  time  of  occurrence  of  the  pulse. 

The  AO  cell  at  Pj  is  uniformly  illuminated  with  this  light  distribution.  The  AO 
cell  has  a  transit  time  m."  t  and  is  fed  with  a  signal  consisting  of  pulses  of  width 
At  and  period  n't  where  m  is  the  residual  modulus. 

We  describe  the  pattern  at  I*  hy 


s,(x,t)  =  '  ‘  t  x-nm.'.t  +  t )  . 

1  n 


(2.  SO) 


I.eav  ing, 
the  form 


f  i  nd  s 


fl 


integrated  at  P  to  produce 


f 


n 


i 

(d-nm) Ax  . 


(2.11) 


From  (2.11),  we  see  that  the  spatial  location  of  the  output  peak  at  P  is  propor¬ 
tional  to  the  desired  residue  K  of  .1  modulo  m. 

m 


To  provide  parallel  conversion  of  .1  into  R  for  N  different  bases  m .  ,  we 

mi  l 

modify  s^(t)  to  be  a  sum  of  signals,  each  consisting  of  pulses  of  width  t  and 

different  periods  T.  =  m.'.t  and  different  associated  carrier  frequencies  f.  within 
11  i 

each  't  for  each  rrc  .  This  causes  each  modulus  to  appear  spatially  separated  at  P9, 
where  we  block  the  dc  and  one  first  order  term. 

At  P  ,  we  place  N  vertical  grating.  These  cause  the  correlations  s  (*)  s 
1  U  1  i 

of  s^  with  different  modulo  m .  signals  to  appear  at  different  spatial  frequency 
locations  in  P.^.  The  P  pattern  is  thus  the  N  correlations  (*)  Sjj  on  N'  differ¬ 
ent  horizontal  lines.  The  horizontal  position  of  each  correlation  corresponds  to 

R  .  with  the  m.  encoded  in  the  vertical  direction.  Thus,  the  sv-tem  of  Fig.  2.1b 
mi 

provides  parallel  conversion  of  the  input  integer  data  .1. 

The  system  of  Fig.  2.1b  was  assembled  using  an  AO  cell  with  s^  =  '38  MHz,  band¬ 
width  of  +5  MHz,  and  T  =  32.25  ..sec.  Moduli  m.  -  7,  9,  10  were  chosen  with  t  =  606 

—  c  1 

nsec  and  associated  periods  T.  corresponding  to  pulse  frequencies  f ^  =  34.7,  38,  and 
41.3  MHz.  These  correspond  to  spatial  frequencies  f =  55.9,  61.29,  and  66.61. 

These  result  in  a  reasonable  separation  of  2.62  mm  between  gratings  at  P Appropri¬ 
ate  grating  spatial  frequencies  were  chosen  for  use  at  P.,  (15,  35,  55  cy/mm) . 

The  corresponding  outputs  are  shown  in  Fig.  2.17.  The  position  of  the  peak  on 
each  line  corresponds  to  the  residue  R  .  for  the  corresponding  moduli  m.  when  evalu- 

TTU  1 

a  ted . 

An  older  AO  cell  of  lower  optical  quality  was  used  in  these  experiments.  The 
final  AO  cells  arrived  too  late  in  the  program  to  permit  the  Tl  and  APAR  experiments 
to  be  performed  on  then.  Similarly,  equipment  delays  and  lack  of  funds  prohibited 
full  use  of  the  electronic  support  system  until  after  completion  of  the  contract 
period.  These  experiments  still  demonstrate  the  use  of  the  AO  cell  and  much  of  the 
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support  electronics  as  well  as  the  T1  AO  correlator  in  a  new  application  (residue 
arithmetic)  addressing  accurate  optical  data  processing  for  APAR  use.  In  the  next 
section,  the  COC  simulator  is  used  with  the  TS1  architecture  to  show  computation  of 
a  far-field  noise  pattern  adaptive  in  angle  and  frequency. 

2.10  TSI  EXPERIMENTS 

In  this  section,  we  provide  simulated  verification  of  the  TSI  system  of  Fig.  2,12 
and  discuss  the  associated  post-processing  required  to  calculate  the  adaptive  weights 
from  its  output  data  for  diverse  noise  cases. 

To  demonstrate  use  of  the  TSI  system  of  Fig.  2.12,  two  random  uncorrelated 

2 

noise  sources  with  Gaussian  statistics,  zero-mean  values  and  with  variances  o  =  1.0 

2  0  0 
and  '  =  1,4,  located  at  =  +23  and  =  -22  with  different  bandwidths  of  0.352" 

and  0.234"  (2/3  of  the  first  source)  were  produced  using  the  simulator  of  Sect.  2.5. 

We  let  ,  0,  +-•  in  frequency  or  radian  bandwidth  correspond  to  samples  0,  256  and 

512  respectively.  The  spectrum  of  the  two  noise  sequences  are  shown  in  Fig.  2.18. 

For  the  first  noise  source,  its  bandwidth  was  90  points  out  of  256  centered  at  the 

point  256  +  45  =  301;  whereas  for  the  second  noise  source  its  60  point  bandwidth  was 

centered  at  286. 

The  angular  locations  of  the  noise  sources  were  verified  by  performing  the  wide¬ 
band  correlation  corresponding  to  the  top  line  in  the  output  pattern  of  Fig.  2.12. 

The  results  are  shown  in  Fig.  2.19.  The  relationship  between  the  delay  ;  in  radians 
and  the  angle  of  the  noise  source  is 

■  =  (2"d/>)  sin  (2.32) 

In  the  simulations,  we  assumed  d  =  100  i/2"  or 


povr 
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t  =  100  sin  0  (2.33) 

We  employed  a  detector  integration  time  T  =  512  samples.  The  two  noise  source 
angles  0  =  +23®  and  -22®  correspond,  using  (2.33),  to  delays  of  +39  and  -37 
points  respectively.  The  data  in  Fig.  2.19  confirm  these  noise  source  angle 
locations.  The  amplitude  of  the  two  noise  peaks  also  agrees  with  the  intended 
signal  strengths  with  the  second  noise  source  (with  the  larger  variance)  pro¬ 
ducing  the  larger  peak. 

The  N  =  9  multi-channel  outputs  in  Fig.  2.12  are  shown  in  Fig.  2.20. 

The  bandwidth  of  each  band  pass  filter  h  was  chosen  to  be  0.1171'  (30  points). 

This  resulted  in  nine  separate  frequency  filters.  From  Fig.  2.18  we  see  that 
the  first  source  lies  at  23®  and  has  zero  response  beyond  0.352";  whereas  the 
second  source  lies  at  -22®  and  has  zero  response  beyond  0.234tt.  These  values 
are  in  exact  agreement  with  the  expected  results.  The  strength  of  the  response 
in  the  different  frequency  bands  in  Fig.  2.20  are  approximately  equal  correspond¬ 
ing  to  the  approximately  flat  spectrum  (Fig.  2.18)  of  the  noise  sources  over  their 
bandwidth.  Note  the  wideband  (WB)  output  on  the  top  line. 

We  now  briefly  consider  the  use  of  the  TSI  system  (Sect.  2.6)  and  the  vari¬ 
ous  post-processing  algorthms  (Sect.  2.8)  for  different  noise  source  cases  in 
APAR. 


For  the  case  of  a  single  source  at  one  frequency,  low  correlation  SNR  is 
expected  because  of  the  zero  bandwidth  signal.  We  expect  some  signal  bandwidth 
in  practice  and  thus  the  TT  correlator  will  yield  adequate  estimation  of  the 
target  angle  .  because  of  the  self-compensating  algorithm  feature  (Sect.  2.8) 
the  noise  frequency  estimate  is  not  of  concern  in  this  case.  Multiple  sources 
of  the  same  frequency  perform  similarly  and  if  their  bandwidth  is  low,  their 
correlation  peaks  will  be  wide  and  accurate  discrimination  of  the  different  angles 
of  the  noise  sources  may  be  difficult.  For  such  cases,  beam  forming  techniques  may 


be  preferable. 


POWER 


For  the  case  of  a  single  noise  source  with  multiple  frequencies,  we  have  a 
wide-band  jammer.  In  this  case,  the  self-compensating  feature  is  not  appropriate. 
Instead,  we  use  the  TSI  system  and  the  modified  projection  method  with  additive 
features  to  produce  a  wider  null  at  the  center  frequency  (if  the  noise  source 
frequencies  are  in  one  band)  or  the  advanced  projection  technique  in  other  cases. 

Choices  of  the  slit  widths  (or  the  effective  bandwidth  of  the  different  filters) 
in  the  TST  mask  and  other  issues  such  as  the  combined  use  of  the  wide-band  corre¬ 
lation  output  and  the  multi-channel  narrow-band  correlations  (the  first  to  accurately 
determine  the  angular  target  locat ions  and  the  combination  to  deconvolve  the  fre¬ 
quency  response  of  each  target)  were  described  in  Sect.  2.6. 

2.11  SUMMARY  AND  CONCLUSION 

Our  new  CO C  work  on  APAR  processing  has  resulted  in  many  new  algorithms  and 
system  architectures  as  well  as  in  the  fabrication  of  many  components  for  the  COC 
processor.  In  the  past  one  year  we  have: 

(1)  Developed  a  new  COC  concept  using  1-D  AO  cells  rather  than  2-D  SLMs 
(because  they  are  more  easily  fabricated  and  more  readily  available) 
and  TI  rather  than  SI  correlators  (because  their  longer  integration 
times  allow  better  noise  statistical  estimates)  (Sect.  2.2). 

(2)  Fabricated  and  performed  initial  testing  of  two  AO  cells  (Sect.  2.3). 

(3)  Developed  a  new  simulator  to  handle  multiple  phased  array  signals 

with  adaptivity  in  space  and  frequency,  to  compute  weights  and  calculate 
corrected  antenna  patterns,  plus  conventional  Fourier  transform  and 
correlation  routines  (Sect.  2.5). 

(4)  Devised  a  new  adjunct  antenna  concept  that  allows  full  space  and  fre¬ 
quency  AFAR  data  to  be  obtained  using  only  a  two-channel  processor 
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(5)  Developed  a  new  hybrid  time  and  space  integrating  AO  correlator 
architecture  that  combines  the  best  features  of  the  TI  and  the 

SI  systems,  and  provides  a  2-1)  display  of  the  angle  and  frequency 
location  of  the  far-field  antenna  noise  pattern  from  the  two- 
element  adjunct  antenna  data  (Sect.  2.6). 

(6)  Designed  and  fabricated  a  hardware  electronic  support  system  and 
designed  and  nearly  completed  fabrication  of  a  multi-purpose 
computer-driven  electronic  support  system  for  complex  noise 
source  scenarios  (Sect.  2.7). 

(7)  Developed  a  self-correcting  post-processing  system,  a  new  fast 
iterative  modified  projection  method  technique  and  an  advanced 
projection  concept  to  compute  the  adaptive  weights  for  narrow- 
band  and  wide-band  jammers  that  differ  in  both  angle  and  fre¬ 
quency  (Sect .  2.8) . 

(8)  Performed  an  experimental  verification  of  a  TI  correlator  and 
used  it  to  demonstrate  residue  arithmetic  computations  in  a 

new  time  position  coding  scheme.  Residue  arithmetic  is  an  advanced 
technique  whereby  numerical  computations  can  be  performed  in 
parallel  with  no  carries  and  with  reduced  dynamic  range  require¬ 
ments.  This  is  necessary  to  realize  an  accurate  optical  computer 
(Sect.  2.9). 

(9)  Demonstrated  the  use  of  the  TSI  system  for  computation  of  the  2-D 
space  and  frequency  output  antenna  pattern  from  an  adaptive  array 
using  an  adjunct  antenna  (Sect.  2.10). 

Thus,  in  conclusion,  we  have  devised  a  new  COC  optical  processor,  demonstrated  and 
simulated  the  key  points  of  the  system,  begun  development  of  the  necessary  post¬ 
processor,  and  fabricated  the  necessary  components  and  support  system  to  fully 


study  this  COC  technique  for  APAR. 
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CHAPTER  3  ITERATIVE  OPTICAL  PROCESSOR  (IOP) 

3 . 1  INTRODUCTION 

In  this  chapter,  we  describe  our  recent  research  on  the  IOP  system.  Since 
this  system  has  been  outlined  in  Chapter  1  and  described  in  depth  elsewhere  [2-6], 
we  concentrate  on  our  recent  progress  on  this  system.  In  Sect.  3.2,  we  describe  the 
new  system  design  adopted  for  fabrication  and  in  Sect.  3.3  we  describe  the  system 
we  fabricated  during  the  past  year.  The  basic  and  advanced  versions  of  the  TOP 
simulator  and  their  use  in  APAR  processing  and  in  the  analysis  of  the  TOP  are  then 
described  in  Sect.  3.4.  Several  specific  IOP  operating  issues  such  as  accuracy, 
performance,  corrections,  convergence,  convergence  rate,  selection  of  the  accelera¬ 
tion  factor,  scaling  of  the  eigen-values  of  M,  etc.  are  discussed  in  Sect.  3.5. 
Initial  experimental  demonstrations  of  the  IOP  systems  use  in  APAR  applications 
are  then  presented  in  Sect.  3.6.  A  summary  and  our  conclusions  are  then  advanced 
in  Sect.  3.7. 

3.2  SYSTEM  DESCRIPTION 
3.2.1  Error  Sources 

In  our  prior  work  [2-6],  the  basic  IOP  system  was  developed  and  studied.  Our 
analysis  of  this  system  showed  three  major  error  sources  that  would  limit  the  ac¬ 
curacy  and  performance  of  the  system.  Fig.  3.1  shows  the  simplified  IOP. 

Imaging  the  LEDs  vertically  onto  the  mask  without  cross  talk  was  found  to  be 
a  severe  problem  as  was  uniformly  illuminating  each  row  of  the  mask  with  the  light 
output  from  one  LED.  Our  solution  to  these  problems  was  to  employ  fiber  optic  cou- 
;  ’  i  ,g  between  Pj  ..ndP?.We  chose  ten  LED  elements,  a  10  x  10  mask  and  thus  a  100 


-  81  - 


Figure  3.1  Schematic  diagram  of  the  iterative  optical  processor  (IOP). 

element  fiber  optic  bundle.  These  numbers  were  chosen  because  such  a  system  was 
of  adequate  size  to  allow  it  to  be  fabricated  at  reasonable  cost  and  to  be  used  in 
processing  APAR  problems  of  reasonable  complexity  and  to  allow  unforeseen  fabrica¬ 
tion  problems  to  emerge. 

The  non-linear  response  (current  input  versus  light  intensity  out)  of  the  LEDs 
tested  made  amplitude  modulation  unattractive  without  complex  circuits  to  correct 
for  the  LED  non-linearities.  An  investigation  of  laser  di  de  sources  showed  that 
they  had  much  better  linearity  over  a  very  large  dynamic  range,  but  repented  attempts 
to  obtain  a  linear  array  of  these  elements  were  unsuccessful.  We  thus  chose  to  use 
pulse  width  modulation  (PWT1)  of  the  LED  sources. 

With  these  major  changes  made  in  the  system  design  (pulse  width  modulation  and 
fiber  optic  interconnections),  we  then  considered  four  other  system  component  features 


that  would  affect  performance  of  the  final  svstcm. 


Non-uniformity  in  the  saturation  level  of  the  LEDs  existed  together  with 
the  non-uniform  response  of  the  linear  photo  diode  detector.  Non-uni  term  cou¬ 
pling  of  the  LED  sources  to  the  fiber  optic  system  and  non-uniform  outputs  from 
the  fiber  optic  array  (due  to  polishing  differences  at  the  ends  •.  f  the  elements 
in  the  fiber  optic  bundles)  were  two  other  problems.  We  chose  to  correct  for 
the  non-uniform  I.ED  saturation  levels  and  the  vertical  non-uniformitv  in  the  mask 
by  a  dynamic  RAM  in  the  electronic  feedback  system  to  adjust  the  gains  of  the  differ¬ 
ent  LEDs  on-line.  All  four  errors  and  their  residual  values  (after  RAM  correction) 
were  measured  and  can  be  corrected  for  by  fabricating  a  fixed  correction  mask  to 
be  placed  behind  the  adaptive  mask  B  at  plane  in  the  system. 

Electronic  noise  in  the  detector  and  associated  support  electronics  was  another 
problem  of  concern.  We  note  that  there  are  two  components  of  this  noise:  a  fixed 
signal  pattern  due  to  the  dark  current  profile  across  the  linear  photo  detector 
array  and  a  time  varying  component  due  to  other  noise  sources.  We  corrected  for 
the  fixed  pattern  detector  noise  by  using  a  time-multiplex  differential  output 
scheme  (see  Sect.  3.2.4)  in  which  a  single  complex  multiplication  is  divided  into 
two  parts,  each  is  run  sequentially  and  the  sequential  outputs  are  subtracted. 

The  above  system  design  and  system  operational  changes  provide  a  quite  accurate 
I  OP  system  for  APAR  use.  Residual  errors  are  still  expected,  depending  upon  the 
accuracy  to  which  the  corrections  are  fabricated  and  the  magnitude  of  the  time-vary¬ 
ing  noise  sources  in  the  system.  We  have  measured  this  residual  error  and  by  simula¬ 
tion  and  lab  verification  determined  its  effect  on  the  accuracy  of  the  final  W 
estimate  and  the  resultant  null  depth  in  the  adaptive  antenna  pattern.  Those  issues 
are  described  and  quantified  in  Sects.  3.3  and  3.6.  To  prevent  such  error  sources 
from  causing  divergence  of  the  solution  to  the  vector-matrix  equation,  we  include 

tolerance  difference  between  two  estimates  w.  and  w.  ,  When  w . -w .  ,  ,  t  ‘  , 

i  i  +  1  i  i  +  1 
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the  iterative  process  is  stopped.  This  prevents  the  system  noise  from  causing 
the  output  to  diverge  from  the  correct  answer  as  it  is  approached. 

3.3 . General  Purpose  Processor  Considerations 

This  10P  system  is  a  quite  general  purpose  optical  processor  capable  of  solv¬ 
ing  any  general  vector-matrix  equation  of  the  appropriate  form.  Since  optical 
systems  are  in  general  special  purpose,  this  architecture  is  quite  unique.  To  main¬ 
tain  the  general  purpose  features  of  the  IOP  and  to  allow  its  use  in  other  associated 
APAR  tasks  and  in  other  applications,  several  design  considerations  were  made.  These 
choices  were  also  made  to  decrease  cost  of  the  prototype  system  and  the  associated 
electronics.  We  have  used  demultiplexers  (DMUX)  at  the  LED  inputs  and  the  photo 
detector  outputs.  This  decreases  the  throughput  of  the  present  system  but  greatly 
decreases  the  associated  LED  drive  and  photo  detector  hardware  necessary.  We  have 
also  only  used  a  1  MHz  clock  frequency  for  the  LED  sources.  Although  higher  speed 
operation  is  possible  (especially  with  laser  diodes),  the  requirement  on  the  associ¬ 
ated  digital  and  analog  hardware  become  quite  cost  prohibitive  for  the  available 
funding  level.  it  is  also  possible  to  perform  addition  of  the  steering  vector  by 
appending  an  additional  row  to  the  covariance  matrix  mask  and  driving  it  with  a  sep¬ 
arate  linear  LED  array.  Instead  of  this,  we  will  perform  this  vector  addition  in 
a  ALL  (arithmetic  logic  unit)  by  digital  multiplication  in  the  electronic  feedback 
loop.  As  noted  earlier,  a  time-sequential  system  for  handling  complex  data  is  also 
used  to  decrease  the  space  bandwidth  product  (SBWP)  requirements  of  the  matrix  and 
vectors  in  the  system,  thus  enabling  more  complex  problems  to  be  addressed  on  the 
present  system.  Although  the  present  svstem  employs  a  fixed  matrix  mask  transparen  v, 
it  can  later  be  rep]a(-e(j  g  v  an  elertroded  spatial  light  modulator  (SLM)  and  thus  he 
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made  fully  adaptive. 
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Thus,  the  present  laboratory  I OP  prototype  system  represents  a  real-time 
optical  processor  for  APAR  of  sufficient  size,  speed  and  complexity  to  enable 
it  to  be  used  for  a  quite  diverse  selection  of  signal  processing  problems.  It 
also  allows  the  accuracy  and  performance  of  the  system  to  be  assessed  and  new 
iterative  algorithms  to  be  incorporated  without  major  hardware  modifications. 

3.2.3  Notation 

The  consistent  notation  to  be  used  in  describing  this  system  is  listed  in 
Table  3.1.  Underlined  lower-case  and  upper-case  letters  denote  vectors  and  matrices 
respectively.  The  letters  a,  B,  and  c  refer  to  actual  optical  and  electrical  system 
parameters  and  are  used  to  describe  the  amplitude  of  specific  functions  in  die  system. 
The  letters  x,  H,  and  v  are  used  to  refer  to  general  algebraic  equations  and  opera¬ 
tions.  The  elements  of  _x,  H,  and  v  are  bipolar,  whereas  for  explicitly  complex¬ 
valued  data  elements,  we  use  s,  M,  and  w.  Various  subscripts  .ire  used  throughout. 
Their  notation  will  be  obvious  in  the  different  cases.  The  detector  output  is  tne 
actual  vector-matrix  product,  but  to  this  we  add  another  vector  to  obtain  the  new  in¬ 
put.  i'or  notational  simplicity,  we  represent  the  detected  vector-matrix  output 
(after  iteration  i)  by  £  and  the  vector  to  be  added  by  £.  This  will  simplify  many 
vector-matrix  equations  throughout  the  text.  The  two  uses  of  £  should  not  cause 
difficulty  when  taken  in  context.  We  also  denote  versions  of  a  at  different  itera¬ 
tions  by  a.  or  a^.  A  similar  notation  is  used  for  the  algebraic  quantities  x,  HI,  and 
y . 

1.2.4  Complex-Valued  Data  Handling 

The  steering  and  adaptive  weight  vectors  £  and  w  as  well  as  the  covariance  matrix 
M,  have  complex-valued  elements.  Conversely,  the  LKD  outputs  and  the  mask  trans¬ 


mittance  in  the  IDP  are  real  and  positive  valued  quantities.  Thus,  considerable 
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attention  is  necessary  to  enable  the  10P  to  accommodate  the  complex-valued  data 
elements  necessary  l'or  the  APAR  problem. 

After  studying  various  techniques  by  which  to  operate  on  complex-valued  data 

in  a  real  and  positive  IOP  system  and  considering  the  statu1  noise  detector  pattern, 

a  sequential  bipolar-data  technique  was  chosen.  To  describe  this  technique,  we  first 

recall  that  the  elements  b  of  the  actual  optical  mask  used  at  P_  must  satisfy 

mn  1  2 

0  <  b  n  1  to  be  normalized.  We  choose  to  operate  on  bipolar  complex-valued  data 
by  performing  the  operation 
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in  the  system,  where  all  quantities  in  (3.1)  are  bipolar  and  subscripts  r  and  i  de¬ 
note  the  real  and  imaging  parts  of  the  corresponding  vector  or  matrix. 


As  noted  earlier,  we  use  two  sequential  processor  cycles  to  perform  a  complex¬ 
valued  vector-matrix  product.  In  the  first  cycle,  the  inputs  are  the  positive  parts, 

x  and  x+,  of  and  the  outputs  produced  are  the  positive  parts  of  v,  v+  and  v+.  The 
— r  — i  —  •—  -^-r  •‘-i 

same  optical  matrix  ji  is  used  for  all  vector-matrix  operations  until  M  changes.  Re¬ 
call  the  matrix  H  is  related  to  M  by 


H  =  a  -  M),  (3.2) 

where  I  is  the  identity  matrix.  We  denote  the  minimum  and  maximum  element  values  of 
H  by  h  and  h  respectively.  The  optical  matrix  B  used  is  thus  generated  from  the  matrix 
li  hv 
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or  in  algebraic  terms 

y_  =  (h-h)  [B  a+  -  B  a"]  +  h  '  x  ,  (3.8) 

—  —  —  —  —  mm 

where  the  vector  inputs  at  each  of  two  sequential  cycles  are  represented  by  a+ 
and  £  ,  thus  effectively  realizing 


The  basic  processor  is  thus  a  bipolar  vector-matrix  multiplier.  It  is  de¬ 
scribed  by  the  signal  flow  graph  of  Fig.  3.2.  The  necessary  pre-  and  post-processing, 
scaling  and  normalization  operations  are  performed  in  the  micro-processor  electronic 
feedback  system  (Sect.  3.2.5).  Another  technique  for  handling  complex  data  is 
described  in  Chapter  4  in  conjunction  with  the  wavelength  diversity  processor  (WDP) 
version  of  the  TOP  system. 


3.2.5  Micro-Processor  Feedback  System 


A  micro-processor  feedback  system  is  necessary  to  perform  the  pre-processing 
of  the  vector  x  and  the  matrix  H,  the  post-processing  to  perform  the  different  opera¬ 
tions  in  (3.7)  and  the  associated  scaling  and  biasing  needed  for  the  data.  The  micro¬ 
processor  system  also  performs  the  necessary  addition  of  the  vectors  £.»  v  or  to  the 
bipolar  vector-matrix  product.  It  also  controls  the  F.AM  correction  circuits  for  the 
photo  detector  output  and  the  LED  input.  The  LED  driver  in  the  output  demultiplexer 
and  associated  other  system  controls  are  also  included  in  the  electronic  feedback 
system  together  with  various  storage  and  readout  features  that  enable  various  vector- 
matrix  products  and  iterative  output  products  to  be  stored  in  digital  memory,  displayed 
sequentially  on  an  oscilloscope  in  single-step  or  continuous  mode,  etc.  Details  of 
this  aspect  of  the  system  are  included  in  Sect.  3.3.  A  schematic  of  this  hardware 
electronic  feedback  loop  is  included  in  Fig.  3.3. 
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Iterative  Optical  Processor:  Hardware  Block  Diagram 
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3 . 3  SYSTEM  FABRICATION 

3.3.1  Optical  Vector-Matrix  Multiplier 

A  schematic  diagram  of  the  optical  vector-matrix  multiplier  is  shown  in 
Fig.  3.4.  The  input  consists  of  a  linear  array  of  10  RCA  SC-1002  LF.Ds  which 
typically  emit  1  mw  for  a  50  ma  drive  current  at  940  nm.  The  LEDs  are  mounted 
on  0.15  inch  centers  along  a  copper  block  1.5  inches  long  and  held  in  place 
with  silver  epoxy.  A  fiber  optic  element  manufactured  to  our  specifications  is 
used  to  distribute  the  light  from  each  LF.D  across  the  matrix  mask.  At  one  end 
it  contains  an  array  of  10  apertures,  over  which  the  LEDs  are  positioned.  Each 
of  these  apertures  contains  a  bundle  of  ten  glass  fibers,  each  0.001  inch  in 
diameter,  the  bundles  are  anchored  in  brass  collars  set  in  an  A’1  block  over  the 
LEDs.  The  fibers  branch  outward  to  form  a  2-D  array  of  10  x  10  or  100  fibers 
that  are  set  in  a  stainless  steel  face-plate.  The  fibers  are  arranged  in  a 
10  x  10  rectangular  matrix  with  a  eenter-to-center  spacing  of  0.014  inch  verti¬ 
cally  and  0.0375  inch  horizontally  between  fibers.  The  fiber  ends  were  polished 
and  the  entire  unit  potted  in  black  latex  in  an  casing.  The  optical  matrix 

mask  is  sandwiched  between  the  fiber  optic  element  and  the  detector  array  as 
shown  in  the  figure.  The  detector  used  is  a  Centronics  LD-20  silicon  photo  diode 
array  which  contains  20  elements  each  measuring  4  x  0,9  mm  on  0.95  mm  (0.0374  inch) 
centers.  The  detector  elements  and  the  size  of  the  matrix  mask  are  chosen  such 
that  each  detector  element  sums  all  of  the  light  emerging  from  a  column  of  the 
matrix  mask.  This  provides  the  summation  over  the  m  rows  in  each  column  and  pro¬ 
duces  the  desired  output 


c  =  a  b 
n  m=1  m  mn 


(3.10) 


ELECTRO-OPTICAL  VECTOR-MATRIX  MULTIPLIER 


■ 


IT 
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A  photograph  of  the  optical  vector-matrix  multiplier  is  shown  in  Figure 
3.  6.  Right  to  left  are  the  LED  array  (sealed  in  white  RTV  compound)  bolted  to 
the  Tiber  optic  element,  matrix  mask  and  detector  array.  (The  components 
are  separated  for  clarity  in  the  figure.)  The  fiber  optic  element,  matrix 
mask  and  detector  are  mounted  on  mechanical  translation  stages  and  tilt 
positioners  to  facilitate  alignment  of  all  components. 

3.3.2  Electronic  Feedback  System 

The  electronic  feedDack  system  is  composed  of  the  following  four 
subsystems: 

(i)  LED  ^ulse  width  modulator  (PWM),  demultiplexer  (DMUX)  and  current 
drive, 

(ii)  time- integrating  photodetector  and  multiplexer  (MUX), 

(iii)  D/A  and  A/D  converters  and 

(iv)  microprocessor  controller. 

These  subsystems  and  their  interconnections  are  shown  in  the  block  diagram 
of  Figure  3.6. 

The  LED  board  [subsystem  (i)]  consists  of  a  clock  (reference  frequency, 

fQ),  digital  PWM,  current  drive  and  DMUX.’  Each  element  of  the  non-negative 

input  vector  a  =  {am}  is  fed  (by  the  microprocessor)  in  digital  form 

to  the  LED  board  and  is  converted  by  the  PWM  and  current  drive  to  a  current 

pulse  of  duration  a  / f  The  amplitude  of  the  pulse  is  fixed  so  that  the 
mo 

peak  power,  Pm>  radiated  by  the  m-th  LED  is  constant  and  its  energy  varies 

linearly  with  a  Amplitude  modulation  of  the  drive  current  would  result  in 
m 

nonlinear  distortion  due  to  the  exponential  (current-power)  response  character¬ 
istic  of  the  diode.  The  time-varying  output  from  the  m-th  LED  is  therefore: 

1  m(t)  =  Pm  Rect  '3.11> 

ana  the  energy  it  provides  is  Pmam/fQ.  All  M=10  LEDS  are  selected  in  turn 

by  the  DMUX  and  are  pulsed  in  this  manner.  The  time  T  required  to  enter  the 

vector  a  is  data  dependent  and  equals 

,  M=1 0 


(3.12) 


(mT Pif*oq  loioajad 
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Complex  data  is  handled  by  separating  the  input  vector  into  its  positive 
and  negative  real  and  imaginary  components  and  time-mul ti plexino  the  optical 
vector-matrix  multiplier.  Since  the  time  to  input  a  number  is  proportional  to 
its  value  and  at  least  one  of  the  two  non-negative  components  of  each  element 
of  the  array  ix  will  be  zero,  the  throughput  rate  for  bipolar  data  is 


f  m=l 
o 


;x  i,  where  f  =  1  mhz  and  -128 
m  o 


127  for  the  laboratory  system 


fabricated. 


The  total  power  incident  on  the  n-th  detector  element  is 


M=  1 0 

0n(t)  -Lit 
n  M  m=l 


mn  bmn 


(3.1?) 


where  I  (t)  is  specified  by  (3.1'}:  3 


b  -  are  the  values  of  the  matrix 
nr. 


mask;  tmn  is  the  non-uniform  spatial  transmittance  of  the  fiber  optic  element 
and  L  is  the  optical  efficiency.  Since  the  LEDS  are  pulsed  sequentially  in 
our  system  to  simplify  hardware  design  and  fabrication,  the  photodiodes  must 
time-integrate (3.13)  over  the  entire  period  T. 

The  photodetector  support  electronics  [subsystem  (ii)]  contains  a  bank 
of  10  resettable  op-amp  integrators  which  sum  each  of  the  photocurrents  over 
this  time  interval.  At  time  t=T,  the  output  voltage  of  the  n-th  integrator 


f  o„<t>  *  ytn  « 


(3. 1^) 


where  rp  is  the  responsivity  of  the  n-th  photo  diode,  0p(t)  is  given  by  (3.13), 

C  is  the  integration  capacitance  and  i  (t)  is  the  noise  current  introduced  by 

the  n-th  photodiode  and  support  electronics.  The  variations  in  fiber  optic 

transmittance  t  ,  the  peak  power  P  radiated  by  the  LEDS  and  the  detector 
mn  m 

responsivity  r  will  be  compensated  for  by  a  static  correction  mask  placed 
in  contact  with  the  matrix  mask  and  by  electronic  correction  to  be  performed 
in  the  microprocessor.  Letting  C  be  the  signal  component  of  (3.14),  we  obtain 
for  the  output 


-  ./  \  M=1( 


a  b 
m  mn 


(3.15) 


where  r  is  the  average  detector  responsivity  and  p  is  the  average  peak  power 
radiated  by  the  LEDS.  In  practice,  the  gain  (rLp/f  MC)  is  set  equal  to 
unity  by  proper  selection  of  the  integration  capacitance  C  and  adjustment 
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of  the  clock  frequency  f  .  The  analog  voltages  vn(T)  —  Cn  are  read-out 
sequentially  by  the  MUX,  converted  into  digital  form  by  the  A/D  converter 
in  subsystem  (iii)  and  fedback  into  the  microprocessor. 

The  microprocessor  controller  [subsystem  ( i v ) ] ,  which  is  depicted  in 
the  hardware  block  diagram  of  Figure  3.7,  is  responsible  for  the  schedulinn 
and  execution  of  all  operations  in  the  computation  cycle  of  the  IOP.  These 
operations  include:  LED  preprocessing  and  correction,  optical  vector- 
matrix  multiplication,  detector  post-processing,  bias  removal,  rescaling 
and  vector  addition.  The  microprocessor  contains  a  custom  designed 
arithmetic  unit  consisting  of  a  16  bit,  300  ns  TRW  multiplier,  a  16-bit 
arithmetic  ljogic  unit  (ALU)  and  a  16K  random  access  memory  (RAM)  with  a 
row-column  address  structure.  We  have  arranged  these  units  to  simplify 
software  control  and  data  paths. 

The  control  section  sequences  the  internal  data  transfers,  arith¬ 
metic  operations  and  input/output  as  well  as  the  operation  of  the  LED  and 
detector  boards  (i)  and  (ii)  in  Figure  3.6. It  consists  of  a  26K  RAM  to  store 
the  microprograms  which  operate  the  IOP,  a  Fairchild  9408  LSI  microprogram 
sequencer  to  execute  the  stored  programs  and  a  32  line  instruction  decoder 
which  activates  the  various  control  points. in  the  system.  The  instructions 
are  executed  as  either  evokes  to  activate  control  lines,  conditional 
branches  or  jumps  to  a  subroutine. 

An  extensive  IOP  program  selection  exists  for  component  and  system 
tests  and  for  system  operation.  A  console  (front  panel)  is  provided  to 
load  the  programs  into  microprogram  RAM  and  contains  all  necessary  operator 
controls  to  start,  stop  and  reset  the  microprocessor.  The  microprocessor 
controller  and  front  panel  contain  160  IC's  which  consume  about  50  watts 
of  power.  The  cycle  time  for  any  microinstruction  is  300  ns.  A  photo¬ 
graph  of  the  entire  laboratory  IOP  system  is  shown  in  Figure  3.8. 

3 . 3 . j  System  Operation 

Before  starting  the  microprocessor,  the  vector  data  /  and  xjo)  are 
loaded  into  the  data  memory  (DM)  and  the  scalar  parameters  u>,  £  and  (h-h) 
are  loaded  into  program  memory  (PM)  via  the  front  panel.  The  LED  and 

detector  correction  data,  — }  and  {|— },  are  permanently  stored  in  the  DM. 

m  n 

The  memory  maps  are  shown  in  Tables  3.3  and  3.3.  The  data  memory  (Table  3.?'  contain 
these  data  as  well  as  the  current  bipolar  iterate  x(k),  its  positive  a^ 
and  negative  a^  vector  components,  the  optically  computed  c^  and  vectors 


Column  Address  0-9 


Row  address 


i 

1 


Input  vector 


k-th  iterate  x(k) 


(1/p  } 

LED  correction  rm 


positive  optical  input  i-) 


negative  optical  input  a^ 


positive  optical  output  £, 


negative  optical  output  £2 


photodiode  correction  ^Vr  } 

n 


k-th  bipolar  output  y(k) 


>  I  OP  data 


TaM 0  3.2 


Data  Memory  Map:  DM(0  .:  63,  0  :  IP),  16  bits 


Address  0-255 


256-1023  *< 


Main  program 

Subroutines: 

(i) 

vector  accumulation 

(ii) 

LED  pre-processing 

( i  i  i ) 

Optical  V-M  multipl icatior 

(iv) 

detector  post-processing 

able  3.3 


Program  Memory  Map:  PM(0  :  1023),  ?f  bits 
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and  tne  resultant  Dipolar  vector-matrix  output  yjk)  =  H  x  ( k ) . 

The  program  structure  of  the  IOP  is  divided  into  four  major 
subroutines : 

(i)  Vector  accumulation 


(ii)  LED  pre-processing 

(iii)  Optical  vector-matrix  multiplication  and 

(iv)  Detector  post-processing 

The  subroutines  are  stored  in  program  memory  (Table  2.2)  along  with 
subroutines  for  data  loading,  initialization  and  display  and  special 
routines  for  system  tests  and  calibration.  The  main  program  (which  is 
essentially  a  series  of  calls  to  the  aforementioned  subroutines)  is  located 
in  the  top  portion  of  program  memory.  A  flowchart  of  the  computational 
cycle  of  the  IOP  is  shown  in  Figure  3.S.  Each  block  is  one  of  the  four  sub- 
rountines  mentioned  above.  After  starting  the  microprocessor,  various 
system  initializations  are  performed  [including  yjO)  =  Oj.  The  vector 
accumulation  [subroutine  (i)]  electronically  subtracts  the  input  vector 
from  the  last  bipolar  vector-matrix  product  yjk-1),  multiplies  the  diff¬ 
erence  by  the  acceleration  parameter  w  and  subtracts  the  result  from  the 
last  iterate  xjk-1)  t0  form  the  new  iterate  x(k). 


This  new  iterate  is  separated  into  the  optical  vectors  a^  and  a^  which 

are  then  electronically  multiplied  by  the  LED  correction  factors  {i->  in 

Km 

the  LED  preprocessing  routine  (ii).  Next,  two  vector-matrix  multiplications 


are  performed  optically  on  these  data  by  subroutine  (iii).  First,  the 


corrected  positive  component  a^  of  x  is  fed  to  the  LEDS;  £j  is  optically 
computed  and  stored.  The  corrected  negative  component  a^  of  x  is  then 
fed  to  the  LEDS  and  c ^  is  optically  computed  and  stored.  In  the  detector 
post-processing  routine  (iv),  these  two  vectors  £-|  and  £^  are  electronically 
subtracted  and  multiplied  by  the  detector  correction  factors  The 

M=1 0  n 

result  is  then  scaled  by  the  factor  (h-h)  and  the  bias  h[  \  x^O)] 

m=l 


is  removed  from  each  of  the  output  components  y  (k)  to  generate  the  bipolar 
output  y(k). 

This  procedure  is  repeated  until  terminated  either  manually  (by  the 
console  control)  or  by  a  preset  break-point  in  the  microprogram.  The  seq¬ 
uence  of  iterates  x_(0),  x_(l),  x(2)...can  either  be  displayed  in  real-time 
as  they  are  generated  or  stored  in  PM  for  play  back  and  off-line  analysis. 


i  ’  t.n  '/»  iterates  can  be  stored  in  the  present  system. 


l!):’ 


Start 


Stop 


FIGURE  3.9 


IOP  COMPUTATIONAL  CYCLE  TLO'/nn^i: 
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3 . A  I  OP  SIMULATOR 

Another  new  feature  of  this  IOP  program  was  to  extend  its  use  to  include 
adaptivity  in  both  angle  and  time.  This,  combined  with  the  need  to  handle 
complex-valued  data,  increases  the  required  SBWP  of  the  vectors  and  matrices. 

We  thus  used  simulation  to  study  such  advanced  cases.  We  have  written  an  IOP 
simulator  that  is  quite  interactive.  Two  levels  of  simulator,  a  basic  and  an 
advanced  version,  have  been  written.  Each  is  described  in  this  section  and  then 
several  examples  of  their  use  are  given.  Other  uses  of  the  simulator  are  in¬ 
cluded  in  Sects.  3.5  and  3.6.  All  routines  are  written  in  Fortran. 

The  APAR  .1.  FOR  routine  (Table  3. A)  calculates  the  covariance  matrix  M, 
steering  vector  s  and  the  mask  transmittance  function  1_  "  i*  from  the  given  in¬ 
put  variables.  The  operator  enters  the  number  of  antenna  elements,  the  number 
of  vector  inputs  Co  be  sampled  in  time  and  space,  the  receiver  noise  power  P^ 

as  well,  as  the  velocity,  angle,  position,  number  M,  and  strengths  P  of  the  M 

m 

noise  sources.  This  defines  the  set  of  M  noise  sources  and  the  receiving  array 
geometry.  We  also  specify  the  steering  vector  by  the  location  and  angle  of  the 
target.  From  this,  we  produce  the  2-D  antenna  receiver-processor  model  of  Fig. 

3.10.  This  new  model  has  coordinates  (n,  n')  corresponding  to  the  space  and 
time  tails  respectively,  with  N  adaptive  elements  n  in  space  and  N'  time  taps 
n'  for  each  received  element.  In  Fig.  3.10,  the  time  history  flow  of  the  re¬ 
ceived  signal  at  each  element  in  space  is  recorded  vertically  and  sampled  at  N'  time 

intervals.  A  5  x  5  matrix  or  2-D  space-time  antenna  grid  is  used  to  reduce 
romputat ions . 

We  use  an  incremental  space  increment  (antenna  eler-  •£  spacing)  of  d  =  1/2  = 

0.5  and  time  increments  T  =  0.0003  corresponding  to  a  velocity  v  =  X/AT  =  +  500. 

max  — 
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tabu:  >.4 


Ai'AK.l.  FOR  non  INF 


00 1  00 

INTECER  UMAX,  L.,  L.IIAX ,  N1 ,  N2,  LI,  L2,  NAPS!  1  20),  MAPT  (  1  20) 

(40200 

REAL  V(0  19), TH ( 0  19),P(0  19) 

00300 

REAL  PI,  T,  D,  LMDA,  N 

00400 

REAL  ARG 

00500 

COMPLEX  COV( 1.20,  1  20), STAR  <  1  20  ) 

00600 

PI =3  14159 

00  70  0 

T  -0.  001 

00800 

D=0  5 

00900 

LMU.A-1  0 

ulOOO 

C 

°**THIS  SECTION  OP  CODE  INPUlS  THE  FAf!  FIELD  ANTENNA  PATTERN*** 

01  100 

r 

t  >  1 2C0 

C 

01300 

C 

01400 

C 

01500 

WRITE! 5, 5) 

0 1  600 

5 

FORMAT ( '  ENTER  NUMBER  OF  NOISE  SOURCES  AND  ELEMENTS  IN  VECTOR  ') 

01700 

READ! 5, 10)  MMAX.LMAX 

01800 

10 

FORMAT! 21 ) 

019CO 

DO  35  M- 1 , MMAX 

02000 

WRITE (5, 30)  M 

02100 

30 

FORMAT! '  ENTER  VELOCITY,  ANGLE  AND  STRENGTH  OF  SOURCE  #  ’,12) 

02200 

READ  <5,  32  )  V  ( M ) ,  TH  <  M )  <  P()1> 

023C0 

TH(M)=TH(M)*P 1/180  0 

024C0 

32 

FORMAT (3F) 

02500 

35 

CONTINUE 

02600 

40 

WRITE (5,  45) 

02700 

45 

FORMAT! '  ENTER  VELOCITY  AND  ANGLE  OF  TARGET  ') 

028C0 

READ! 5, 50)  V(0i , TH(O) 

02900 

TH(0)=TH!0)*PI/ 180  0 

03000 

50 

FORMAT ( 2F ) 

03100 

WRITE(5. 70) 

03200 

70 

FORMAT!'  TYPE  RECEIVER  NOISE  POWER  AND  TARGET  GAIN  ') 

0 330  0 

READ ! 5, 75)  N, P!0) 

03400 

75 

FORMAT (2F) 

03500 

C 

**aWE  NOW  MAP  THE  SAMPLED  2D  ANTENNA  PATTERN  INTO  A  ID  VECTOR*** 

03600 

c 

03700 

c 

03800 

c 

03900 

c 

04000 

DO  100  L= 1  - LMAX 

04100 

WRITE (5, 85)  L 

04200 

85 

FORMAT!  '  ENTER  SPACE  TIME  COORDINATE  OF  ELFMENT  #  ',  12) 

04300 

READ  (  5,  90  )  N1,N2 

04400 

90 

FORMAT (21 ) 

045C0 

MAPS ( L ) =N1 

04600 

95 

MAP  T ( L ) =N2 

04700 

1  00 

CONTI  N'JE 

ANTENNA  ELEMENTS  n 


T  =  0.0005 


+  ’ 


'  T~ 
H- 


d  =  V2  =  0.5 


Figure  3.10  2-D  space-time  antenna  model  (n,  n’). 

9 

In  terms  of  conventional  units,  for  a  1  GHz  radar,  f  =  10  Hz,  A  =  0.3  m,  T  = 

-4 

1.5  x  10  sec  and  v  =  +500  m/sec .  These  units  and  relationships  arise  from 
the  relationship  for  the  phase  of  a  signal  ,at  element  n  and  time  tap  n'  for  a 
target  at  angle  ■■  and  velocity  v  and  a  time  delay  T  between  antenna  taps  given 
by 


'.‘n  n,  =  2"  (r.ds  in  <•  +  2vTn')/A.  (3.16) 

When  d  =  )/2  and  1  =  1,  (3.16)  becomes 

»  =  fn  sin  ••  +  n ' ( v/v  \  ,  (3.17) 

n ,  n  [_  y  max  / 

where  v  is  the  maximum  target  velocity  given  bv  v  =  A/4~'. 

This  (n,n')  mapping  is  thus  recorded  and  a  2-D  map  of  (n,n')  samples  is  pro¬ 
duced.  The  vector  s  ,  matrices  I  -  M  and  the  (n,n')  maps  are  stored  in  a  disk 
file.  To  rap  this  2-D  data  onto  the  l-D  vector-matrix  system,  we  must  map  (n,n') 
into  a  1-1)  vector  .  .  This  routine  accomplishes  this  and  stores  the  associated 
maps  on  a  disk  file.  The  routine  to  compute  s  and  M  is  given  in  Table  3.5.  In 


Fig.  3.10,  we  see  that  a  25  element  5x5  space-time  .tonne  grid  is  possible. 


-  106  - 


TABLK  3.5 


ROl'TINE  TO  COMPUTE  S  &  H 


04800 

C 

04900 

c 

05000 

C 

05100 

c 

0520 

053C0 

05400 

05500 

05600 

c 

057C0 

05800 

150 

05900 
06000 
06100 
0620 O 
06300 
06400 
065C0 
06600 
06700 

1  55 

068C0 

069C3 

070G0 

07100 

07200 

200 

07300 

07400 

202 

07500 

07600 

205 

0770  0 
07800 
07900 

300 

08000 

08100 

082C0 

08300 

084CO 

310 

08500 

08600 

09700 

08800 

320 

00900 

330 

09000 
09  1  CU 
09200 

3  50 

***COMPUTES  steering  vector  and  covariance  matrix*** 


DO  150  L= 1 . LMAX 

ARG=2  a p I  * ( MART  <  L ) * V  <  0 ) *T+MAPS ( L ) *  D*S IN  <  TH ( 0  > ) ) /LMDA 
STAR  <L)=CMPLX (SORT (P  <0) ) ,  0.  0) * 

1  CI1PLX  (COS(ARG),  SIN(ARG)  ) 

CONTINUE 
WRITE! 5, 155) 

FORMAT (  '  STEERING  VECTOR  COMPUTED  ') 

DO  200  Ll=l. LMAX 
DO  200  L2= 1 . LMAX 

IF  (  L 1  =  =L2  )  COV(Ll,  L2)=Cl!l-'LX(N,  0.  0) 

IE  <  L 1 #L2 )  COV (LI, L2  >  =  (0.  0.  0  0  > 

DO  200  M= 1 1  MMAX 

ARG=2»PI*(  ( MAP  T ( L2 ) -MART ( L 1 >  )«V(M)*T  + 

1  ! MAPS ( L2 ) -MAPS ( L 1 > )*D*S IN ( T H ( M ) ) > /LMDA 

COVILl. L2 ) =COV < LI , L2 ) +CNPLX  < P ( M ) ,  0.  0 ) * CMPLX ( COS ( ARG ) , -SIN(ARG) ) 

continue 

DO  202  LI  =  1  - LMAX 
DO  202  L2=1,LMAX 

COVILl. L2 )  =  ( -1 . 0. 0.  0 > *CQV (LI.  L2) 

IF  <L1  =  =L2>  COVILl, L2)=C0V(L1,  L2>  +  ( 1.  0.  0.  0) 

CONTINUE 
WRITE! 5. 205) 

FORMAT!'  COVARIANCE  MATRIX  COMPUTED  ') 

OPEN < UNI T=20. FILE='  STEER.  DAT') 

FORMAT ( 2F ) 

DO  310  L=l. LMAX 

WR I TE  <  20. 300  >  S  T  AR ! L ) 

CONTINUE 

OPEN ( UNI T=21, FILE='  COVAR  DAT') 

DO  320  Ll=l . LMAX 
DO  320  L2- 1 . LMAX 
WRITE(21. 300)  COVILl. L2) 

CONTINUE 

OPEN  I  UN  I  T =22.  F  I  l.E=  '  AN  TMAD  DAT  '  ) 

DO  350  L 1 = 1 ■ LMAX 

WRITE <22, 330)  MAPS! LI ). MART (LI ) 

FORMAT (21 ) 

CONTINUE 

STOP 

END 
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With  10  LEDs,  we  can  represent  10  steering  vectors  and  10  adaptive  weights.  Since 
complex-valued  data  is  necessary,  each  element  requires  four  bipolar  entries,  the 
plus  and  minus  real  and  imaginary  parts  of  each  vector  or  matrix  element.  With  the 
time-sequential  bipolar  data  technique  (Sect.  3.2),  only  two  values  are  necessary 
per  cycle  per  element.  Thus,  for  each  scenario,  we  can  use  5  or  the  25  possible 
sample!-  in  the  2-D  (n,n')  matrix  antenna  model  of  Fig.  3.10. 


From  he  APAR  .1.  FOR  routine,  we  obtain  a  2-D  to  1-D  mapping  of  (n,n')  into 
,  with  any  5  of  the  25  possible  2-D  sampled  matrix  positions  in  Fig.  3,10  possible 
for  entry  into  the  IOP.  The  third  routine  in  the  basic  simulator  is  ITFR.l.  FOR 
(Table  3.6).  This  implements  the  vector-matrix  iterative  equation 


w .  =(I-M)w.  +s 

-i  +  l  -  —  -i 


(3.18) 


in  the  form 


-i  +  l  =  -i  ”  (M  «i  -  «“)• 


(3.19) 


whore  in  (3.19),  M  w.  is  the  new  vector-matrix  product,  N  w  -  s  is  the  difference 
that  is  to  be  reduced  to  zero  and  .  is  the  acceleration  factor.  We  discuss  selection 
of  i  in  Sect.  3.5  and  its  choice  in  speeding  convergence  of  the  algorithm  and  in 
forcing  the  algorithm  to  converge.  In  general,  if  the  eigen-values  are  small,  i  is 
chosen  large,  since  tills  increases  the  eigen-values  and  makes  convergence  faster. 
Conversely,  if  the  eigen-values  -ire  large,  i.e.  1,  /  is  chosen  small  to  decrease  all 

eigen-values  to  be  less  than  i  thereby  insuring  convergence  of  the  iterative  algorithm, 
l!  one  noise  source  is  much  larger  or  smaller  than  another,  large  and  small  eigen¬ 
values  will  result  and  we  must  thus  scale  the  matrix  M  such  that  the  largest  eigen¬ 
value  fits  within  the  dynamic  range  of  M  and  such  that  all  cigen-values  are  less  than 


1 


1.  When  the  dynamic  range  of  the  eigen-values  is  large,  a  loss  in  system  noise 
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TABLE  3.6 


tTER.l.  FOR  RQL'T I NE 


00100 

C 

»TH  I S  PROGRAM  SIMULATES  THE  IDEAL  APAR  ITERATIVE 

00200 

c 

PROCESSOR  ANO  EITHER  PRIUlS  THE  WEIGHTS  AT  EACH 

ITERATION 

00300 

c 

OR  DISPLAYS  AN  ISOMETRIC  PLOT  OF  THE  FAR  FIELD 

P ATI  ERN.  ** 

00400 

c 

00500 

c 

006  CO 

INTEGER  LI, L2,  LMAX,  INCRE,  ITER, Kl, K2, KMAX 

00700 

INTECER  TEST,  DISP,  MAPS ( 1  :  20),  MAPTO:  20) 

00800 

REAL  PI,  T,  D,  LMDA 

00900 

REAL  ARG 

01000 

REAL  F ( 0 : 64, 0: 64) 

01100 

REAL  RAT  101 , RATI02 

01200 

COMPLEX  STAR ( 1 : 20 ) , COV< 1 . 20,  1 : 20), WOLD ( 1 : 20 ) , WNEW ( 1 : 20) 

0 1 30  0 

COMPLEX  E<0: 64, 0: 64) 

0 1  400 

P I =3  14159 

015C0 

T=0  001 

01600 

D=0  5 

01700 

LMDA^l  0 

01800 

c 

0 1 900 

c 

02000 

c 

02100 

WRITE (5, 1 ) 

02200 

1 

FORMAT (  '  ENTER  LENGTH  OF'  VECTOR  AND  THE  DISPLAY 

INCREMENT ' ) 

02300 

READ <  5, 2)  LMAX,  INCRE 

02400 

n 

c- 

FORMAT  <  21 ) 

02500 

0PEN(UNIT=20, FILE=  '  STEEk.DAT') 

0260  0 

3 

FORMAT <2F) 

027C0 

DO  4  L 1 = 1 , LMAX 

02800 

READ ( 20, 3 )  STAR (LI ) 

02900 

4 

CONTINUE 

03000 

OPEN ( UNI T=21 , FILE= '  COVAR.  DAT  '  ) 

03  ICO 

DO  5  L 1 = 1 , LMAX 

03200 

DO  5  L2=l i LMAX 

03300 

READ (21, 3 )  COV  <  L 1 , L2  > 

03400 

5 

CONTINUE 

03500 

WR  I  TE  (  5,  10) 

03600 

10 

FORMAT ( '  TYPE  1  TO  PRINT  W  AND  0  TO  PLOT  E  ') 

03700 

READ ( 5, 12)  DISP 

03800 

12 

FORMAT  < 1 1 ) 

03900 

ITER=0 

04000 

OPEN ( UN I T=22, F ILE=  '  ANTMAP.  DAT'  ) 

04100 

DO  14  Ll=l > LMAX 

04200 

READ(22,  13)  MAPS ( L 1 ) ,  MAPT ( L 1 ) 

04000 

13 

FORMAT ( 2 I ) 

044C0 

14 

CONTINUE 
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adaptivity  may  result  by  the  above  reduction  process  to  accommodate  the  dynamic 
range  of  the  M  mask.  This  issue  will  be  addressed  in  later  phases  of  this  work. 

The  routine  ITER. 1 .  FOR,  reads  the  disk  data  from  the  output  o'  APAR.1.  FOR, 
calculates  the  w.  and  also  K('  )  . 1 1  each  iteration  due  to  the  present  w.  estimate. 

— l  —  t 

The  latter  operation  is  performed  by  the  routine  in  Table  3.7. 

Many  sample  noise  scenarios  were  run  using  this  simulator.  In  all  eases,  the 

isometric  outputs  of  the  F(  )  adaptive  .antenna  pattern  are  shown  after  different 

numbers  of  iterative  cycles.  For  the  different  eases  to  be  described  here,  the 

output  F.(  )  antenna  pattern  that  results  from  the  applied  adaptive  weights  are 

shown  after  100  iteration  cycles.  U'e  first  consider  the  case  of  adaptivity  in 

■  only.  We  assume  one  noise  source  at  12^  and  the  signal  or  steering  vector  at 
0 

23  with  receiver  noise  power  I’  =  0.1  and  noise  power  P^  =  0.4.  The  results  are 
shown  in  Fig.  3.11  after  100  iterations.  The  pattern  shows  a  null  at  12^  (the 
angle  of  the  noise  source)  as  expected.  ,  In  Fig.  3.12,  we  show  similar  output 
results  for  the  case  of  four  noise  sources  at  +  30(l  and  +60^  each  with  power  P  =  0.1 
for  the  case  of  a  boresight  steering  vector  after  100  iterations.  Again  good  re¬ 
sults  are  observed  using  the  iterative  algorithm.  Nulls  appear  at  the  four  correct 
noise  source  angles  in  the  adaptive  pattern  shown. 

U'e  next  considered  the  case  of  three  noise  sources  .it  v  =  250,  =  -45  ; 

v  =  230,  =  +45^;  and  v  =  -250,  =  +45 For  these  eases  of  a  pattern  with 

noise  sources  at  different  angles  and  frequencies,  2-1)  adaptivity  is  required. 

After  computing  the  associated  M  and  the  resultant  adaptive  weights,  the  resultant 
adaptive  antenna  pattern  !(  )  shown  in  Fig.  3.13  was  obtained.  This  pattern  shows 
nulls  at  tiie  correct  velocity  and  angle  coordinates  corresponding  to  the  three  noise 
sources  noted  above.  It  is  thus  in  excellent  agreement  with  the  results  predicted 


be  thcorv. 
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TABLE  3.7 

R0UTJNE_  TO  COMPUTE  I-  ( ) 


04500 

C 

***THIS  BLOCK  COMPUTES  Thr-  ADAPTIVE  ANTENNA.  WEIGHTS  a# 

04600 

C 

047C0 

C 

04S00 

C 

04900 

c 

05000 

DO  15  L!=l . LMAX 

051  CO 

WOLD  1  L 1  )  =  <0.  0.  0.  0) 

05200 

1  5 

CONTINUE 

053C0 

16 

DO  50  L 1  =  1 1  LMAX 

0540  0 

WNEW < LI ) =ST  AR (LI) 

05500 

DO  50  L2= 1 , LMAX 

05600 

WNEW  (  L  1  )=UINEW<L1  >+COV(Ll  .  L2  )  *  WOLD  ( L2  ) 

057  00 

50 

CONTINUE 

05300 

IE  <MOD< ITER. INCRE)#0)  GO  TO  87 

059C0 

IE  ( D I SP= -0  )  GO  TO  62 

060C0 

WR I TE ( 5, 52)  ITER 

06100 

52 

EORMAT ( '  ITERATION  #  ',!) 

062 00 

DO  60  L 1 = 1 » LMAX 

06300 

WRITE! 5,  55)  WNEW (LI  ) 

06400 

55 

EORMAT (2F) 

06500 

60 

CONTINUE 

06600 

GO  TO  79 

C6700 

r 

***THE  FOLLOWING  SECTION  CALCULATES  THE  FAR  FIELD 

06800 

C 

ANTENNA  PATTERN  AND  PRODUCES  AN  ISOMETRIC  PLOT*** 

068C0 

c 

07000 

c 

07100 

c 

07200 

c 

073C0 

62 

KMAX=64 

07400 

DO  72  K1=0,  KMAX 

07500 

DO  72  K2=0, KMAX 

07600 

E(Ki,  K2)  =  <0  0,  0.  0) 

077C0 

DO  70  L  1  =  1 , LMAX 

07800 

RAT  1 0 1 =FLQAT  <  K 1 ) /F  LOAT (KMAX ) 

07900 

RAT  1 02=E LOAT ( K2 ) /FLOAT (KMAX ) 

08000 

ARG-2*P I * ( MAPS ( L 1 ) * D*SIN ( P I * ( RAT  101 -0  5) )  + 

08  1  00 

1  MART (LI ) * ( RAT  102-0  5))/LMDA 

0820  j 

E(K1,  K2)=E(K1, K2 ) +WNEW ( L 1 > *CMPLX < COS ( ARG '  ,  -SIN(ARG) ) 

083C0 

70 

CONTINUE 

08400 

F(K1,  K2)=REAL(E(K1, K2 ) *CON JG ( E ( K 1 .  K2 ) )  ) 

08500 

F(K1, K2)=5+AL0G10(F(K1,  K2)  ) 

08600 

72 

CONTINUE 

08700 

CALL  PL0T2D  <  64«  64,  F) 

08800 

79 

WRITE (5, 80) 

08900 

80 

FORMAT (  '  T yPE  1  TO  CONTINUE  IT'  RATION') 

09000 

READ(5,  85)  TTST 

09100 

85 

FORMAT ( I ) 

09200 

IF  ( TEST -~0 )  STOP 

093C0 

87 

I TER= I TER+ 1 

09400 

DO  90  L 1 ~ 1 , LMAX 

09500 

WOLD (LI  ) =WNEW ( L 1  ) 

096C0 

90 

CONTINUE 

09700 

GO  TO  16 

09800 

END 

suo 
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This  simulator  is  thus  a  general  routine  with  i  no  roast'd  I  leX.ibiliLv  and 
i  no  roast'd  real-time  user  interact  iuii  that  inoludos  additivo  receiver  noise,  tile 
abilitv  to  oomputo  M,  I  -  M,  the  iterative  weights  w^  at  each  iteration  and  the 
resultant  K(  ")  adaptive  far-field  noise  pattern  tliat  results  iron  the  given  w 
and  s.  An  advanced  version  of  this  simulator  includes  the  effects  ol  additive 
and  niu 1 1 i p 1 iea t  i ve  noise  error  sources  in  the  hardware  10P,  computes  the  additve 
weights  and  K(  )  witli  and  without  these  error  sources,  the  processing  gain  or 
output  SNR  of  the  antenna  pattern  and  heme  the  effect  of  TOP  error  sources  on 
the  performance  of  the  APAR  system.  We  describe  a  demonstration  of  this  advanced 
simulator  in  Sect.  .3.6  to  analyze  the  accuracy  and  performance  of  the  IOP  labora¬ 
tory  demonstration  experiments. 

3.3  C'VSTKM  OPERATION 

In  this  section,  we  consider  several  different  aspects  of  the  IOP  system  and 
their  effect  on  the  system's  design  and  its  performance.  In  Sect.  3.5.1,  we  describe 
the  setup  of  the  [OP,  measurements  of  its  accuracy  and  the  corrections  employed  in 
the  experimental  s vs  ten  assembled.  In  Sect.  3.5.3,  we  discuss  convergence  of  the 
IOP,  selection  of  the  acceleration  factor  ■  to  improve  the  rate  of  convergence  and 
how  to  determine  <  from  the  eigen-values  of  M.  Simulation  results  are  then  included 
to  verify  these  issues. 

1.5.1  System  Calibration 

I'p  ’ii  ti'i  I'ivin,:  tlie  diltcrent  svsten  components  .and  assembling  the  IOP,  we  per- 
:  >r— iti  n.inv  component  and  svstet  .  necks  and  cal  ihrat ions.  In  this  sub-section,  we 
highlight  tin  major  calibration  tests  and  the  eventual  corrections  used.  In  each 
case,  special  n i c ro-proces.se r  sub-routines  were  written  to  provide  the  individual 
required  svstem  and  component  tests  noted. 
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We  first  measured  and  corrected  for  the  LED  non-uniform  saturation  levels. 

A  test  was  made  by  pulsing  each  LED  on  in  sequence  for  256  x  4  '  1  msec  (the  PWM 

width  corresponding  to  the  largest  input  value).  With  no  mask  present,  the  ten 
fiber  optic  outputs  in  each  row  were  measured  in  sequence  with  the  detector  butted 
against  the  mask.  An  array  of  10  x  10  =  100  voltage  values  was  then  obtained. 

The  current  into  each  LED  was  then  adjusted  such  that  each  LED  produced  the  same 
10  photo  detector  voltage  output  readings  of  5  volts  corresponding  to  one  half  the 
saturation  level  of  t lie  output  from  the  detector  board.  This  ensured  that  the  system 
was  operating  in  its  linear  range  (one  half  of  saturation  voltage).  These  LED 
corrections  were  then  entered  into  a  RAM  in  the  micro-processor  electronic  feed¬ 
back  system  and  used  for  all  future  corrections. 

After  these  corrections  were  made,  the  residual  system  errors  were  measured. 

For  the  100  output  photo  detector  readings  after  RAM  correction  of  the  LEDs,  we 
measured  a  standard  deviation  in  the  data  of  0.093  with  a  mean  value  of  6.065. 

This  corresponds  to  a  non-uniformity  of  0.8a.  In  our  advanced  simulator,  we  in¬ 
cluded  this  as  a  mu  1 1 ipl ic.a t i ve  error  in  the  mask  transmittance.  This  measure¬ 
ment  and  correction  technique  corrects  for  the  differences  in  LED  saturation  levels 
and  the  differences  in  transmittance  of  the  different  fiber  optic  rows.  This 
10  x  10  =  100  element  error  matrix  with  0.82  non-uniformity  could  be  corrected  for 
by  a  fixed  mask  in  contact  with  the  B  matrix  at  P?.  However,  computer  plotters 
and  other  optical  film  recorder  systems  to  which  we  presently  hive  access  cannot 
correct  for  such  an  accuracy  of  less  than  1"  in  amplitude.  In  the  latter  phase  of 
this  project,  we  will  use  our  film  recorder  to  produce  a  correction  mask  to  the 
necessarv  accuracy.  Many  svstem  corrections  are  necessary  to  the  film  recorder 
before  this  performance  is  possible  however. 

We  then  measured  the  total  output  detector  noise.  This  measurement  was  per¬ 
formed  with  no  light  incident  on  l he  photo  detector  and  with  coaxial  cable  connected 
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to  the  detector  board  outputs.  We  measured  a  10  mV  peak-peak  variance.  This 
represents  the  temporal  noise  variance  on  a  fixed  noise  pattern.  The  latter 
noise  is  removed  by  the  2-cycle  sequential  subtraction  technique  used.  The  former 
temporal  noise  represents  an  0.16  error  source.  We  entered  this  noise  error 

source  as  additive  noise  into  our  advanced  simulator  to  determine  its  effect  on 
the  computed  weights  and  the  resultant  E(;)  adaptive  antenna  pattern.  To  do 
this,  we  produce  a  random  vector  with  a  peak-peak  noise  of  10  mV'  and  add  this  to 
the  vector-matrix  output.  This  low  measured  value  of  temporal  detector  noise  is 
so  small  that  no  additive  noise  distribution  (Gaussian,  etc.)  noise  modelling  and 
analysis  was  performed.  The  peak  detector  voltage  is  6V  and  its  LDR  is  625:1. 

We  then  measured  the  linearity  of  the  system  by  increasing  the  l.F.D  pulse  width 
from  0-500  Usee  (with  100  mA  drive  current  used)  and  measure  the  photo  detector  out¬ 
put.  The  linearity  of  the  system  was  essentially  perfect  as  expected  because  of  the 
excellent  linearity  of  the  external  integrator  used. 

A  system  light  level  budget  analysis  was  then  performed  to  determine  the  I.F.D 
power  necessary.  We  denote  the  output  detector  voltage  by 

V  =  RP'.t/c,  (3.20) 

where  K  is  the  responsivitv  of  the  detector  (0.5  A/W) ,  P  is  the  input  power  from  one 
LED,  t  =  4  x  128  =  512  500  usee  is  the  maximum  input  pulse  duration,  and  C  =  0.75 

1>K  is  the  external  capacitor  used  in  the  integrator.  For  V  =  5 v  (saturation  of  the 
photo  detector),  we  find  P  =  1.5  x  10  ^  W  of  power  from  one  LED  to  he  necessary  (one 
LED  on  is  the  case  to  be  considered  for  saturation  of  the  output  detector  for  normal¬ 
ization).  The  i.ED  used  provides  1.5  x  10  ^  W  of  intensity.  However,  a  large  (approxi¬ 
mately  a  factor  of  10)  loss  occurs  in  coupling  the  Li .11  outputs  to  the  fiber  optic 
connector.  This  is  due  to  the  0.1"  spacing  used  between  the  I  .Mi  s  and  the  l  Put  optic 


connector . 


This  spacing  was  chosen  to  improve  the  uniformity  of  the  light  into  each 
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liber  optic  bundle.  An  additional  factor  of  10  loss  occurs  when  the  LED  drive 
currents  are  adjusted  to  produce  equal  outputs.  The  quartz  window  on  the 
detector  and  other  losses  add  to  reduce  the  light  reaching  the  photo  detector. 

We  then  adjusted  the  pulse  width  such  that  a  5  volt  saturated  output  voltage  was 
obtained  from  the  photo  detector  when  T  =  512  p  sec  was  used.  The  above  values 
satisfy  these  power  requirements  and  sensitivity  values  for  the  system  and  were 
thus  used  in  the  final  design. 


These  additive  0.167.  and  mult ipl i cat ive  O-S7  residual  errors  left  after 
RAM  correction  were  then  used  in  the  advanced  simulator  and  a  mean  square  error 
in  the  adaptive  weights  of  0.17  was  computed.  This  small  error  in  the  w^  due 
to  the  residual  additive  and  multiplicitive  noise  error  sources  is  so  low  that 
it  is  expected  to  produce  no  significant  effect.  We  discuss  this  issue  further 
in  Sect.  3.6.  We  note  for  now  that  the  range  or  spread  in  the  eigen-values  of 
M  can  effect  this,  that  additive  noise  (detector)  can  effect  the  maximum  eigen¬ 
value,  and  that  multiplicative  noise  (mask)  can  cause  stability  problems  and 
osc ill  at  ions . 


3.5.2  I OP  Convergence 


In  formulating  an  I0P  solution  for  the  space  and  time  adaptive  APAR  problem, 
we  write  the  L  receiver  inputs  (complex-valued)  as 


x  = 


(3.21) 


where  the  I.  choices  can  be  any  of  the  N  N  space-time  taps  in  the  2-D  antenna  arra\ 
model.  The  covariance  matrix 


M  = 


r* 

Fax 
)  — 


(3.22) 
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is  then  computed.  We  then  estimate  the  largest  eigen-value  of  M  from 

I"  L  :  .2]  1/2 


The  TOP  then  computes 


-*■  =  ...  ...  ;  m .  . 

max  x  =  l  j  =  1  i] 


(3.23) 


-  w  + 

1  Is*  -  M  w,  \ 

-k 

l-  -  ~k) 

=  a  - 

'  M )  v.  +  1  s  * . 

—  — k  — 

(3.24) 

For  (3.24)  to  converge  to  the  solution 


w  =  M  ^  s  , 


(3.25) 


we  require  the  eigen-values  \  of  M  to  satisfy 


0  '  1  - 


(3.26) 


for  ■  =  1,  .  .  .  ,  L.  For  maximum  monotonic  convergence  at  the  fastest  possible 
rate,  we  choose 


and  we  approximate  the  number  of  iterations  k  necessary  by 


k  l/log(l  -  P  •)  . 

max  ^  \  r  / 


where  P  is  the  receiver  noise, 
r 


(3.27) 


(3.28) 


To  prove  (3.27)  and  (3.28),  we  consider  a  two  element  antenna  with  receiver  noise 
power  P  and  a  single  wide-hand  jammer  at  -j  with  power  Pj.  For  this  case 
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M  = 


P  +  P, 
r  1 


exp  (-jv) 


Pj  exp  (+.]  , )  ' 


P  +  P, 
r  1 


where  :  =  ( 2  '  d  /  '• )  sin  ■  .  The  eigen-values  *  of  H  are  the  solution  of 


det[M  -  IX]  =  0 


for  the  case  in  (3.29),  the  solution  of  (3.30)  yields 


3  -%  +  Pl  +  - 


2  2 

P  -  2P  P, 
r  r  1 


=  P  +  P 
2  r 


1  -  /(Pr  +  Pj)2  -  Pr2  -  2V 


If  P  P.  , 

r  1 


If  pr  pr 


>1  -  2Pj  +  Pr 
a2  Pr. 


1  =  A2  =  Pr  +  Pr 


In  the  usual  case,  P  P  ,  and  (3.26)  and  (3.31)  yield 


or  approximately 


(2P1  +  Pr)  ' 


0  2,P(  •  2, 


or  for  i-l  (no  scaling) 


P1 


(3.29) 

(3.30) 

(3.31a) 

(3.31b) 

(3.32a) 

(3.32b) 

(3.33) 

(3.34) 

(3.35) 


(3.36) 
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Continuing  the  proof  of  (3.27)  and  (3.28),  we  next  determine  A  for  this 

max 

case.  Since  M  is  Hermetian,  A  is  given  bv  (3.23), 

-  max 


A  =  2P,  +  P  =  A 
max  1  r  1 


(3.37) 


for  •  P^.  For  the  case  of  one  jammer  of  power  P  and  1.  adaptive  weights  each 

with  noise  P  , 
r 


A  =  i.P,  +  P  .  (3.38) 

max  1  r 


We  next  consider  the  k-th  iterative  estimate  of  the  weights  w^  and  compare  it 
to  the  final  steady  state  solution  w  to  determine  the  choice  of  <  and  its  relation 
to  A.  For  a  boresight  b  =  oj  target  direction,  s*  =  [1,  1  ]T  and  is*  =  [  t,  'i]T  and 
we  find 


Pl(l__-  e  Jv)  +  Pr 

(,  .  ,-v)  h  _  ,,„P  -  1 

(l  -  e"jV\  ( 1  - 

,  \k  -  1 
'*21 

2(P1P1  *  ?r2) 

2t 

“  2 

1 

— * 

+ 

T) 

U~e+i'\  h- 

'  ‘  k  ~  ^  i 

L  2(PlPr+Pr2) 

■>} 

1 

?  '» 

*"  2 

J 
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Substituting  (3.32),  we  find  that  as  k  approches  infinity,  the  terms  in  k  approach 
zero  and  in  steady  state 


>l(l  -  e'j,)/2P1»r~] 

Jlf1  "  e+j')/2PlPr_ 


(3.40) 


Comparing  (3.39)  and  (3.40),  we  see  that  for  a  rapid  convergence  we  can  set  the  accel¬ 
eration  factor  in  (3.24) 


i  =  ]/• 

max 


1/A 


(3.41  ) 


(3.27).  Then,  the  terms  in  (1 


in  (3.39)  equal  zero.  The  remaining 


J21 


transient  terms  in  (3.39)  become 


Q  -  exp (+j  ,)j  Q  -  >2./>^k  !/2Pr 


(3.42) 


and  are  thus  decreased  in  magnitude  (after  k  iterations)  by  the  factor 


or 


[M1  -  b)'(a  -  -2)lk 

'J1  -  x2y{i  -  v’ijyk  ■ 


(3.43) 


As  .•  approaches  /„*  the  choice  in  (3.41)  most  improves  the  rate  of  convergence  of 
w,  to  the  steady  state  solution  in  (3.40). 


Having  derived  (3.27),  we  now  consider  the  proof  of  (3.28).  We  first  find 

the  performance  measure  of  our  adaptive  array  as  the  SNR,  the  ratio  of  the  power 

in  the  direction  of  the  target  of  strength  to  the  sum  of  the  noise  power  in 

the  antenna  due  to  the  M  noise  sources  of  strengths  P  at  angles  r'  .  Thus,  at  itera- 

m  m 

tion  k  of  the  I0P, 

.  2 

SNR(h)  w  t - p°  i|k(f))T2  •  (3.44) 

m  m  1  k  V  raj  1 


In  (3.44),  we  denote  the  antenna  pattern  at  iteration  k  bv  F.,  (:).  This  is  easilv 

k 

found  from  the  weights  w^  at  iteration  k  bv  a  simple  Fourier  transform  ns 


V-)  = 


I. 

•  =  1 


wR(.) 


exp 


!""  j  2"  '•  d  sin  f'J  / '  , 


(3.45) 


In  our  performance  measure  of  the  use  of  the  10P  system,  specifically  the  iteration 
technique,  we  define  the  process  gain  PC  of  the  system  as  the  improvement  obtained 
in  the  antenna  SNR  as  a  function  of  the  number  of  iterations.  We  thus  define  l’C.(k) 
as  a  ratio  of  the  SNR  at  iteration  k  to  the  initial  SNR  at  iteration  zero,  i.e. 


PC  (It) 


SNR(k) 
SNR (0)  ' 


(3.46) 
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We  now  consider  the  proof  of  (3.28)  by  considering  the  effect  of  the  choice 


of  •  on  Ptl.  We  denote  the  number  of  iterations  necessary  to  achieve  the  given 


PC  without  scaling  by  k  and  the  number  of  iterations  necessary  to  achieve  the 
same  PC  with  scaling  by  k^..  Then,  from  (3.43),  we  can  relate  k^.  to  k^,  by 


(3.47) 


where  and  '•  9  are  the  largest  and  smallest  eigen-values  of  M.  '  ^  and  > 

.ire  less  than  1;  from  (3.32),  >  =  P^_  for  '  ^  •>  >0;  and  from  (3.27),  =  l/">. 

This  shows  that 


(3.48) 


Thus,  choosing  i=  l/>  will  cause  a  shift  in  the  PC  versus  k  graph.  This  also 

verifies  (3.28)  and  shows  the  proper  choice  of  .■  is  most  important  when  P  .  =  1 

or  •  .In  other  cases,  choosing  <  =  1/*,  speeds  up  convergence  bv  de- 

min  max  1 

creasing  the  number  of  iterations  necessary  by  the  factor 


i<bT[T 


(3.49) 


We  next  investigated  the  above  issues  using  the  simulator  (Sect.  3.4).  We  first 

considered  the  effect  of  scaling  on  I’F(k)  for  the  case  of  one  jammer  at  10^.  We 

assumed  receiver  noise  1’  =  0.001,  signal  power  =  0.001  and  different  noise  inter- 

runee  powers  I*  =  0.1,  0.01  and  0.001.  We  calculated  PC(k)  with  and  without  scaling 

versus  iteration  number  k.  For  these  sealing  cases,  we  used  •  =  l/'^  and  for  the 

other  cases  .  =  1 .  In  these  cases,  ■  was  less  than  1  and  scaling  bv  .  onlv  affected 

max 

the  convergence  rate  rather  than  whether  the  1  OP  would  converge  to  the  steady  state 
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solution.  In  all  oases,  the  use  of  sealing  improved  the  number  of  iterations 
necessary  for  the  system  to  converge  to  steady  state.  In  Fig.  3.14,  we  show 
one  example  of  the  outputs  obtained.  In  all  cases,  PC  versus  k  with  scaling 
produced  a  shifted  version  of  the  PC  versus  k  curve  with  no  scaling  as  predicted 
by  (3.47)  and  (3.48).  In  Fig.  3.14,  ;  =  0.05  and  =  0.001  anti  we  '  ind  k^.  =  0.03 

k.  in  (3.47).  Fig.  3.14  shows  that  in  this  case,  the  use  of  scaling  improves  the 
number  of  iterations  necessary  for  convergence  to  steady  state  by  a  factor  of  20 
or  1/0.05  in  perfect  agreement  with  theory.  In  the  case  of  Pj  =  0.001,  an  even 
larger  speed  up  factor  of  125  was  obtained.  The  null  depths  obtained  in  the  three 
cases  were  48,  20  and  11  dB  for  the  cases  of  decreasing  jammer  noise  P  =  0.1, 

0.01  and  0.001. 

In  general,  we  expect  P  P,,.  Also,  in  manv  cases,  we  must  choose  .  to 

m  0 

ensure  that  the  IOP  converges  to  steadv  state.  This  case  arises  when  '•  1. 

max 

From  (3.39),  for  the  case  of  five  adaptive  elements,  as  in  our  antenna  model,  ■ 

max 

will  be  greater  than  1  if  Pj  •  |l  -  P^j/5.  To  address  this  case  of  .  •  1  there 

is  a  relaxation  factor  rather  than  an  acceleration  parameter  as  in  the  case  when 
1),  we  chose  P  =  0.01  =  Pf)  and  P  =  0.1,  1.0  and  10.0  ••  P  .  The  resultant 
PC  versus  k  graphs  are  shown  in  Fig.  3.15.  Those  data  show  that  the  k  necessarv  to 
reach  steady  state  increases  as  P  increases.  Likewise,  the  null  depth  obtained  also 
increases  with  k.  These  results  are  as  expected.  In  this  case,  the  |l  -  .  P^J  eigen¬ 
values  of  the  system  increase  due  to  the  scaling,  and  the  ^1  -  P^j  terms  in  M  de¬ 
crease  less  rapidly.  Other  analogous  cases  were  simulated  and  showed  the  same  trends 

that  increasing  P  is  similar  to  decreasing  P  since  both  increase  k  and  P('.(k). 
m  r 

We  a  I -to  investigated  the  PC.  or  null  depth  possible  for  different  separations 
between  the  target  and  signal  for  the  case  P^  =  0.01  and  P  =  1.0  with  P^  =  0.01. 
Scaling  was  used  and  the  results  are  shown  in  Fig.  3.1b  for  'lie  case  of  k  “  1000 
and  different  angular  separations  between  the  target  (beam  steering  direction) 


Figure  3.14  K fleet  <>f  se.ilinR  on  process  pain  IT,  ns  a  function  of  the  number  of  iterations 
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and  the  jammer.  These  data  show  that  the  null  depth  or  PG  obtainable  decreases 
as  ('  decreases  as  expected  and  that  the  system  provides  adequate  nulls  and  SNRs 
for  fi  ■  1^.  Finally,  we  investigated  the  performance  of  the  IOP  system  with 
multiple  jammers  (M  =  1-7)  at  5^*  increments  centered  at  30*^  with  P^  =  1.0  >>  P^  = 

P^  =  0.01.  We  found  PG  in  steady  state  to  decrease  as  M  increases  as  expected. 

The  initial  E  (()  patterns  were  calculated  and  plotted  for  the  seven  jammer  case. 

The  E(’.‘)  plots  showed  clear  nulls  at  the  noise  sources  as  desired. 

These  simulations  verify  the  use  of  the  10P  and  the  simulator  and  the  impor¬ 
tance  of  proper  selection  of  <  and  its  effect  on  the  performance  of  the  APAR  pro¬ 
cessor  ihrough  SNR,  PG  and  number  of  iterations  k  necessary. 

3.6  EXPERIMENTS 

We  now  consider  one  detailed  IOP  lab  demonstration  experiment.  In  Sect.  3.6.1 
we  describe  the  noise  source  scenario  and  the  associated  21)— 1 n  antenna  model  used. 

The  case  of  an  APAR  processor  with  adaptivity  in  space  and  time  operating  on  complex¬ 
valued  data  is  chosen.  In  Sect.  3.6.2,  expressions  for  s,  M,  and  w  and  the  associated 
optical  masks  are  described.  In  Sect.  3.6.3,  experimental  IOP  laboratory  data  on  the 
optically  computed  weights  is  presented  and  compared  to  the  exact  values  obtained 
by  simulation.  The  performance  of  the  IOP  system  for  APAR  is  best  obtained  by  computing 
the  resultant  E(0)  adaptive  antenna  pattern  produced  by  applying  the  optically  computed 
weights  to  the  antenna.  The  SNR  obtained  and  the  depth  of  the  location  of  the  nulls 
produced  by  the  adaptive  system  is  the  performance  criteria  used.  To  assess  the  accu¬ 
racy  of  the  optical  IOP  system,  the  theoretical  SNR  and  null  depths  for  a  perfect  pro¬ 
cessor  and  for  a  processor  with  the  IOP  errors  were  digitally  calculated.  The  additive 
and  mul t ipl ic it ive  noise  error  source  values  measured  for  the  IOP  (Sect.  3.5)  were  in¬ 
cluded  in  the  advance J  simulator  model.  The  K ( • ' ) ,  SNR  and  null  depths  as  well  as  PG 
expected  from  the  IOP  with  its  given  accuracy  are  computed  and  compared  to  the  exper- 
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3.6.1  Noise  Scenario 

The  APAR  noise  scenario  chosen  was  one  noise  source  (jammer)  with  power 

Pj  =  1.0  located  at  ;•  =  0°  with  a  velocity  =  0  relative  to  the  antenna. 

The  receiver  noise  power  was  =  1.0  and  the  signal  power  =  0.1  and  thus 

P  P^.  We  assume  that  the  signal  is  at  =  45°  with  v^  =  250  relative  to 

the  antenna.  A  5  x  5  antenna  model  with  N  =  5  taps  (five  antenna  elements) 

and  N'  =  5  time  taps  for  the  signal  received  at  each  antenna  element  was  used. 

We  assumed  a  two  element  antenna  with  spacing  d  =  a/2  =  0.5  or  >  =  1  and  time 

delavs  between  taps  of  it  =  T  =  0.0005  =  l/4v  corresponding  to  a  maximum 

max 

target  velocity  of  500  as  described  in  Sect.  3.4.  The  four  adaptive  elements 
chosen  from  this  2-D  antenna  model  were  two  adjacent  antenna  elements  with  two 
adjacent  time  taps  from  each.  The  associated  2-D  and  1-D  mapping  is  shown  in 
Fig.  3.17. 

3 . b . 2  M  And  Mask  Determination 

The  covariance  matrix 

M  =  z  (t)  J  (t)  (3.50) 

with  elements 

m.  .  =  ."  z.  (t)  z.  (t)  dt  (  3.51) 

-tj  t  j 

was  computed  in  our  simulator  assuming  a  narrow-band  jammer,  uncorrelated  receiver 
noise  and  signal  power  P^  much  less  than  the  noise  jammer  power  P^ .  The  matrix  1' 


obtained  (complex-valued  elements)  is 


The  largest  eigen-value  of  M  was  then  determined  as 

=  7.42  =  • 


-  .fs 


* 


3 


(3.53) 

The  acceleration  factor  used  in  the  IOP  was  chosen  as 

(3.54) 

The  complex-valued  M  in  (3.52)  was  then  partitioned  into  its  real  and  imaginary 


. ...  ...  m .  .  m . 

i  =  l  ]  =  1  l.]  lj 


=  1/A  =  0.13. 

max 


bipolar  part.  and  M.  and  these  were  arranged  as 
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(3.55) 


fhc  optical  mask  B  was  obtained  from  H  =  1  -  M  bv  scaling  H  by  dividing  its  elements 
bv  (li  -  h)  =  2  and  then  biasing  this  up  by  addition  of  h/(h  -  h)  =  0  to  each  element 
where  h  and  h  are  the  maximum  and  minimum  values  of  the  elements  of  H.  The  optical 


mask  is  thus  described  hv 
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The  steering  vector  for  =  0°  and  =  0  is 


-o.h4  +  in. *2 
0.99  +  jO.27 
0.62  +  j0.84 
0.27  -  i 0 . 99 


It  is  easily  decomposed  into  its  real  and  imaginary  parts  as 


(3.57) 


s*  s  v 
— r  — l 


T=  [-0.84  0.99  0.62  0.27  0.62  0.27  0.84  -0.99"1  1 .  (3.58) 

t.  _ 


This  vector  was  electronically  entered  into  the  mi cro-prucessor  electronic  feed¬ 
back  system  of  the  IOP  with  the  necessary  scale,  bias  and  acceleration  factors. 

3.6.3  Experimental  Data 

The  above  vector  and  mask  elements  were  entered  into  the  I OF  system.  The 
initial  starting  value  used  for  w  (the  initial  LED  outpu:s  at  the  first  iterative 
evcie)  wds 


— in i L ,  k  =  0  =  101 


(3.59) 


The  IOP  computes 


x  ( k )  =  x  (  k  -  I  '  +  .[v  — M  x(k  -  1)  ] 


(3.60) 


lor  each  iteration  k  where  all  parameters  are  bipolar.  With  eight  elements  for  v 
and  x  and  an  8  x  8  mask  for  H,  the  eight  x  outputs  at  iteration  k  are 


x  (  k )  =  !  x  t  ( k )  , 


,  >:8(k)'T, 


(3.61 ) 


where  x-x,  and  x,-x„  correspond  to  x  and  x.  respectively  because  of  the  data  format 
I  i  >  8  — r  -l 

used.  ITie  complex-valued  output  weights  w(k)  at  iteration  k. 


w  ( k )  =  W||(k),Wj.,(k),w.1j(k),w,1(k)'‘,  (  i .  6  2  ) 
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to  be  applied  to  the  four  time  and  space  taps  in  Fig.  3.17  are  computed  from  the 

x  -x„  outputs  in  (3.61)  as 
1  o 


l,(k) 

=  Xj (k) 

+ 

j*5(k) 

(3.63a) 

^2  (It) 

=  x9(k) 

ix6(k) 

(3.63b) 

2i(k) 

=  x  ( k  ) 

+ 

,i  x  7  (  k ) 

(3.63c) 

.,(k) 

=  x,(kl 

•  ‘  •  o v  •  • '  * 

f  3 . 63d ) 

3.6.4  Analysis 


In  an  actual  system, 


'12 


(t),  z21(t),  z??(t) 


nr 


22 


.J 


antenna  output 


the  weights  would  be  applied  to  the  z(t) 
antenna  outputs  and  the  results  added  to 


[■ll(t). 
form  the 


v 

out 


(t,k) 


2  1 


i=l  j  =  ] 


z  . .(t)w. . ( t )  . 
t)  i.) 


(3.64) 


In  Fig.  3.19,  we  show  the  experimentally  obtained  outputs  corresponding  to 
the  x(k)  data  obtained  from  the  IOP  after  iterations  0,  5,  50.  These  scope  output 
are  the  x(t)  =  Jj<^(k),  .  .  .  ,  x^(k)j  1  outputs  from  (3.60)  at  k  =  0,  3,  and  50. 

The  experimental  values  and  the  exact  values  predicted  by  theory  are  listed  in 


Table  3.8.  From  the  steady  state  (k  =  50)  values  given,  we  find  an  rms  error 


(3.65) 


of  0.26  rms  or  approximately  0.1  of  a  division  on  the  scope  displays  in  Fig.  3.19. 


11  If," 
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TABLE  3.8 

ANTENNA  WEIGHTS  IN  STEADY  STATE  (k  =  50) 


1 

f 

Parameter 

|  Experimental 

Theoretical  (Exact) 

i 

i 

X1 

| 

:  -0.9 

T 

-1.0 

I 

x2 

| 

0.76 

x3 

1  0.35 

0.,  1 

i 

i 

x4 

1 

0.1 

.  - - 

0.05 

[ 

x5 

j  °'3 

0.46 

x6 

|  0.2 

0.13 

; 

x7 

0.75 

0.65 

— 

i _ 

X8 

i  -’-I 

i 

~i — . - - — - - 

-1.1 

We  then  computed  the  adaptive  antenna  pattern  v)  produced  by  the  theoretical 
(exact)  data  and  the  experimental  data,  as  our  comparison  measure  we  use  the  antenna 
SNR 


SNR 


P-oliifco.vol  j; 


(3.66) 


The  experimental  cat  a  gave  on  SNR  =  14. 79  dB,  whreas  the  exact  theoretical  values  pave 
an  SNR  =  14.95  dB.  These  are  in  excellent  agreement.  We  then  used  the  advanced  sim¬ 
ulator  with  the  additive  and  multiplicative  system  errors  included  and  found  an  SNR  = 
15  +  0.75  dB.  Thus,  the  experimental  data  is  well  within  the  range  predicted  by  the 
measured  system  errors  and  is  in  excellent  agreement  with  the  theoretical  values,  with 
an  error  between  experiment  and  exact  theory  of  onlv  0.16  dB. 
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3.7  SUMMARY  AND  CONCLUSION 

The  I  OP  concept  for  APAR  processing  has  proven  to  be  most  excel  lent  The 
highlights  of  our  recent  work  on  this  approach  to  APAR  processing  are  summarized 
below. 

(1)  Design  and  fabrication  of  a  new  10?  system  with  fiber  optic 
interconnections,  pulse  width  modulation  and  a  microprocessor 
electronic  feedback  system  (Sect.  3.3). 

(2)  Use  of  a  new  technique  to  allow  the  IOP  to  operate  on  complex¬ 
valued  data  (Sect.  3.2.4). 

(3)  Development  and  demonstration  of  a  new  10P  simulator,  antenna 
model,  and  2-D  to  1-D  mapping  with  sp.'ce  and  time  adaptivity 
and  including  IOP  error  sources  (Sect.  3.4). 

(4)  Demonstration  and  measurement  of  the  excellent  accuracy  of  the 
IOP  system  to  be  less  than  0.83'  and  with  less  than  17  error  in 
the  antenna  SNR  obtained  (Sect.  3.5.1;  . 

(5)  Development  and  demonstration  of  a  new  technique  to  ensure 
convergence  to  the  steady  state  solution  and  fast  convergence 
of  the  IOP  system  (Sect.  3.3.2). 

(6)  Real-time  demonstration  and  analysis  of  the  laboratory  IOP  system 
for  antenna  adaptivitv  in  space  and  time  (Sect.  3.6). 

The  weights  and  antenna  SNRs  obtained  experimentally  were  within  1  of  the 
theoretical  limit.  This  excellent  system  performance  in  the  IOP  lahoratorv  system 
fabricated  represents  a  ma  jor  new  optical  data  processing  architecture  t  hat  appears  to 


he  quite  attractive  and  useful  for  APAR  processing  and  other  applications. 


CHAPTER  4 


WAVELENGTH  DIVERSITY  PROCESSOR 


4  .  I  INTRODUCTION 

In  tliis  chapter ,  we  consider  the  use  oi  color  to  provide  an  added  dimension 
to  optical  processors  with  attention  to  their  applicability  to  various  aspects 
ot  the  APAR  problem.  In  Sect.  4.2,  we  provide  a  genera’  overview  of  the  wave¬ 
length  diversity  processor.  Specific  details  on  the  design  of  the  system  we 
assembled  are  provided  in  Sect.  4.3.  The  original  reason  16]  tor 
considering  the  use  of  a  wavelength  diversity  processor  was  to  allow  complex- 
value.!  data  to  be  efficiently  processed  in  our  non-coherent  processor  [5] 

In  Sect.  4. 1 ,  we  review  and  summarize  the  various  possible  techniques  for  complex- 
valued  data  processing  in  non-coherent  systems.  The  specific  use  of  the  wave- 
’ ongth  diversity  processor  in  processing  complex-valued  data  is  then  addressed 
in  Sect.  4.3.  A  new  technique  is  introduced  in  Sect.  4.5  that  producer-  the 
complex  output  directly,  without  the  post  processing  required  in  the  system  we 
previously  described  |5|. 

In  subsequent  sub-se  tions,  we  consider  other  advanced  aspects  of  this 
wavelength  diversity  porcessor  for  various  APAR  operations.  included  are  the 
system's  use  as  a  multi-channel  1-D  vec tor-mat r ix  convolver  (Sect.  4.6),  a 
matrix-matrix  mu  ti,'ier  (Sect.  4.71,  a  matrix  inverter  (Sect.  4.8)  and  a 
covariance  matrix  computer  (Sect.  4.9).  Experimental  verification  and 
demonstration  are  included  for  the  applications  in  Sects.  4.5  and  4.7. 

GENERAL  PROCESSOR  DESCRIPTION 

The  general  schematic  for  the  wavelength  diversity  APAR  processor  is 
shown  in  Pig.  4.1.  The  input  at  P  can  contain  I  linear  laser  diode  (ID)  or  light 
emitting  diode  (LED)  arrays,  lor  a  color-separated  wideband  1  iglit  source  and  2-D 
spatial  light  modulator  (SUM)  .vstemj. 
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•  ach  emiitir  ,  at  different  wavelengths  A  =  to  ^  .  Advanced  I,D  technology 
[23,  24]  and  advances  in  optical  communications  [25]  appear  to  allow  fabrication 
ot  such  an  input  LI)  array.  In  the  system  shown,  the  output  light  from  each  I.D 
array  is  spread  horizontally  (to  uniformly  illuminate  each  row  m  of  the  mask  at 
1’,,)  and  imaged  vertical  lv. 

n 

We  denote  the  output  from  P.  bv  a  and  the  transmittance  of  the  mask 

A  "  m  ■ 

at  1’  bv  b  •  Leaving  P„ ,  we  find  the  product  b  a  .  Note  that 
B  mn  B  '  mn  m. 

the  light  from  each  linear  LD  array  is  combined  in  P  ,  thus  super  impcs  i  r.  nil  1. 
input  colors  ■  and  hence  providing  essentially  white  light  illumination  of  11  , 
horizontally.  Behind  P  ,  we  place  a  grating  that  separates  the  different  input 
wavelengths  '•  horizontally  in  Pf  with  all  \  light  going  to  output  PQ ,  all  ^  , 
light  to  P  ,  etc.  at  P  .  Plane  PR  is  imaged  onto  Pc>  The  linear  detector  arrays 
at  P(.  have  sufficient  height  so  that  the  rows  of  P„  are  imaged  within  each  detector 
element  at  P^.  This  provides  the  summation  over  m  of  the  P?  outputs. 

The  final  ?r  output  is  thus 

M 

c  ,  =  :  b  a  (4.1) 

1  n  nil  rn  • 

m 

where  PQ  in  Fig.  2.1  is  c^,  P^  is  c^  ,  etc.  We  will  use  a,  b  and  c  to  denote 
optical  values  and  x,  H,  and  y  to  denote  data  values  in  this  chapter. 

The  P  output  plane  thus  contains  L  linear  1-D  detector  arrays  each  with 
N  elements,  or  equivalently  one  L,  •  N  element  arrav.  This  1 -I)  output  detector 
topology  is  preferable  to  use  of  2-1)  arrays,  both  from  fabrication  considerations 
ami  From  the  ease  with  which  appropriate  horizontal  spacings  between  the  P^, 
etc.  sections  of  the  P  array  can  be  achieved  (compared  to  the  ease  in  obtaining 
the  necessary  vertical  separations  between  rows  of  a  2-D  arrav  or  between  verti¬ 
cal  I  v  stacked  1  —I)  arrays). 


A 
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The  general  format  for  P  ,  P^  and  the  data  incident  on  P(.  is  shown  in 
Fig.  4.2.  As  shown  in  Fig.  4.2,  P  has  L  elements  horizontally  and  M  elements 
vertically.  We  denote  the  width  of  each  element  (1.1)  width)  by  p  ,  their  center- 
to-center  horizontal  spacing  by  p  and  tlieir  vertical  spacings  by  v.  The 

cl 

horizontal  size  of  the  P,  mask  is  denoted  bv  ' p,  and  the  vertical  (as  well  as 

I)  '  b 

horizontal)  size  of  an  element  of  P  is  chosen  as  "v  (this  allows  1:1  imaging 

D 

vertically  of  P  onto  P  )•  The  pattern  incident  on  P  is  shown  in  the  bottom  of 
Ad  C, 

Fig.  4.2.  It  contains  1,  2-D  patterns  each  of  M  X  X  elements.  Each  of  these 

M  X  X  patterns  incident  on  P  is  the  product  of  one  column  of  A  at  one  1  and  the 

P  mask.  The  center-to-eenter  separation  between  each  M  X  X  array  is  'p  .  The 
P)  o 

P  detector  with  its  long  height  integrates  or  sums  the  M  outputs  in  each 
column  at  P  .  In  Sect.  4.3,  we  will  relate  the  'p  and  ip  parameters  in  Fig.  4.2 
to  the  system's  design  factors. 


In  Fig.  4.3,  we  show  the  functional  block  diagram  of  the  wavelength  diversity 
processor  of  Fig.  4.1.  The  system  thus  forms  at  its  output  the  products  of  I. 
input  vectors  a  and  the  matrix  B. 


! 
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Fig.  -*.3  Functional  block  diagram  of  the  wavelength  diversity  APAR  processor. 
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Ho t ore  desc r ibi ng  the  assembled  system  and  its  parameterization  and  design, 

wo  discuss  an  advanced  version  of  the  system  in  which  three  matrix  masks 

and  B.j  (with  center-to-cent or  spacing  '.P  )  are  placed  at  P  (Fig.  4.4).  Consider 

tiie  case  of  three  linear  id)  input  arrays  a^ ,  a^  and  a9  each  emitting  at  a  different 

■  (Fig.  4.4  top).  Mach  l.D’s  output  is  spread  horizontally  to  illuminate  all  of 

P  and  each  linear  l.D  arrav  is  imaged  vertically  onto  P  .  The  P  output  is  again 
B  B  B 

passed  through  a  grating  that  separates  the  colors.  The  spacings  between  the  P 

B 

masks  is  \P  ,  the  same  as  the  separation  between  each  P  output  at  P  .  With 
c  *  C 

proper  choice  of  spacings  and  grating  frequencies,  the  pattern  incident  on  P^ 
is  as  shown  in  Fig.  4.4  (bottom). 


The  P  output  pattern  (Fig.  4.4,  bottom)  obtained  for  this  case  is  best 

•A  .A 

seen  by  considering  each  a  input  separately.  The  a^  input  at  '■  multiplies  all 

A 

three  B  masks  at  and  produces  B  a  at  P  .  Because  the  three  masks  at  P 

B  —  (J  C  B 

A 

are  spaced  hv  '.p  and  because  P  is  imaged  horizontally  onto  P„,  the  B  a.  output 

C  D  O  —  U 

at  P(.  i insists  of  three  terms:  B^a^,  B^a^  and  B^a^  separated  horizontally  at 

P  .  hv  r  .  We  now  consider  the  a,  input  at  P,  at  The  P  ,  output  resulting 

U  '  v‘  1  A  1  c 

^  A  ^ 

from  this  input  is  is  a(  and  again  contains  three  terms:  Jig3!*  — iai  anc*  — 

each  separated  b’/  \p  .  By  proper  choice  of  and  the  frequency  of  the 

c  1  U 

grating,  tin*  '  d  i  1‘  ference  causes  the  B  pattern  at  P(  to  be  displaced  from 

the  B  a  pattern  at  P.  hv  p  which  is  also  the  width  of  an  individual  B  a  pattern 

(with  1:1  imaging  assumed  from  P  to  P  )  .  To  achieve  a  ,'.p^  shift  at  P  due  to 
■ ,  we  must  sat  i  sfy 


fodio 


where  the  terms  involved  arc  discussed  in  Sect.  4.3.  The  P(,  patterns  are  displaced 


due  tc 
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Fig.  4 
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.4  General  patterns  for  the  wavelength  diversity  APAR  processor. 
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Fig.  4,5  Schematic  of  the  experimental  wavelength  diversity  APAR  processor. 
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('.in:  C  inn  i  lit;  with  our  discussion  of  the  anatomy  of  the  P  pattern  in  Fig.  4.4 
(bottom),  the  I’  .  output  If  a.,  due  to  the  a,,  input  at  • is  displaced  by  p  _  from 
the  h  a  output.  in  Fig.  4.4  (bottom),  we  show  the  separate  B  a.  outputs  at 
I’  separated  vertically  to  clarify  what  the  contents  of  P^.  are.  In  actuality, 
the  tiiree  P  outputs  shown  are  added  vertically  on  the  output  detector  to  obtain 
the  desired  output. 


4  .  3  FXPIlk  I  M.KNTAI.  WAV F.l.FNCTH  Ml'LT  iPl.ilXED  PKOCFSSOR 

1  he-  wavelength  diversity  APAR  processor  that  we  assembled  is  shown 
schema t ical ly  in  Fig.  4.5.  Since  linear  LD  source  arrays  were  not  available 
t>'  us,  we  produced  the  equivalent  P  pattern  using  an  arc  lamp  source  and  mask. 
Referring  to  Fig.  4.5,  white  light  from  a  1000  watt  Xenon  arc  lamp 
was  sepurat.u  i.  a  dit  reret.t  colors  by  reflection  from  a 

blazed  grating  (1  at  f()  =  1200  line/mm  spatial  frequency.  Lens  images  the 
slit  vertically  onto  P  , 


I 


( -* .  2) 
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whereas  L2  images  Che  slit  horizontally  across 


dl  +  d2  +  d3 


c2 


(4.3) 


Thus  incident  on  P^,  we  find  separate  wavelengths  or  colors  of  light  A^  =  >Q, 

etc.  at  different  horizontal  locations  across  P  .  At  P^  we  place  a  mask  whose 

transmittance  in  x  =  i  and  y  =  m  corresponds  to  the  desired  a  ,  pattern. 

1  m- 


Lens  L^  performs  1:1  imaging  of  PA  vertically  onto  Pg, 


1  +  1  =  1. 

d5  d6  fc3 


(4.4) 


Lens  L?  also  images  horizontally  onto  PR 


l 

d 


3 


+ 


1  _ 

d4  +  d5  +  d6 


i 


f 


c  2 


(4.5) 


This  effectively  recombines  the  different  A^ 

from  G,  at  P  so  that  P  is  illuminated  with 
IB  B 

,•  from  all  L  sources  is  spread  horizontally 


colors  of  light  incident  on  PA 

essentially  white  light  (i.e.  all 

across  P  as  in  Fig.  4.1). 

B 


If  a  grating  G^  of  spatial  frequency  f^  is  illuminated  with  light  of  wave¬ 
length  >  ,  then  the  first-order  diffracted  spot  in  a  plane  ?A  a  distance  d  from 
G  lies  a  distance  ^  from  the  on-axis  location,  where 


p 


f  A d  . 


(4  ,6a) 


In  the  system  in  Fig.  4.5,  the  spatial  frequency  at  is  f^  and  the  spatial 
frequency  at  l,  is  f^,  vhere 

fo  =  *U _ 

dl  +  d2  dl  +  d2  +  d3 


(4.7) 


-  14  5  - 

since  t‘()  =  f^  if  d^  =  U  and  the  effective  frequency  f^  should  be  used  in  (4.6a) 
with  d  =  d,.  Solving  (4.7)  for  f*  at  the  plane  of  L  and  substituting  this 

•-4 

expression  for  f  in  (4. ha)  with  d  =  d,  ,  we  can  relate  input  •*■ ,  horizontal  dis¬ 
ci  a 

placement  p  in  and  lens  spacings  by 

1’  =  '  fQ  d4(dj  +  d2)/(d1  +  d2  +  d  ).  (4.6b) 

Rewriting  (4.6b)  in  terms  of  the  wavelength  separation  AX  =  X  +  -  A  f  between 

adjacent  vectors  and  the  spatial  separation  ,'P  (see  Fig.  4.2)  between  the  input 

vectors  in  I’  ,  we  obtain 
A 

=  ('■•■)  f()  d4(d|  +  d,)/(d  +  d2  +  dj) .  (4.6c) 

The  1’  masks  used  were  5  T  element  binary  arrays.  They  were  drawn  on 

Amber li th  A3A  masking  film  and  then  photoredueed  20X  onto  Kodak.  1A  high  resolution 

plates.  For  the  final  l\  inputs,  'p  =  2.5mm  and  6p  =0.15  mm  and  TV  =  1  mm. 

A  a  a 

All  cylindrical  lenses  used  had  focal  lengths  f  =  100  mm.  Satisfying  all  focusing 
and  imaging,  conditions,  all  d.  values  in  <4.2)-(4.5)  are  found  and  AX  =  48  nm  in 
(■’..<>>  results.  This  value  appears  to  be  typical  of  the  A>  difference  to  be 
expected  between  1.0  arrays  of  different  colored  light  [26]  With  6p  -  0.15  mm, 

ti 

the  spectral  width  of  each  ■  is  S'  5  0.64  nm.  The  value  lies  between  the  typical 
spectral  line  width  of  an  LTD  (30-40  nm)  and  a  I.D  (0.01  nm)  | 24 ] . 


We  now  return  to  the  details  of  the  rest  of  the  system  in  Fig.  4.5.  Lenses 
ill  I..  image  1’^  with  1:1  magnification  onto  P(., 


d8  *s4  +  ‘s5 


(4.8a) 


d7  d9  +  d10  *  s4  '  fs5’ 


(4.8b) 
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where  f  =  381  mm  is  the  focal  length  of  all  spherical  lenses  used.  This  1:1 
s 

imaging  from  P  to  P  ,  maintains  matched  P  and  P  ,  plane  sizes,  with  the  separations 
B  C  B  C 

between  _B  a  outputs  equal  to  the  separation  p  between  the  mask  patterns. 

It  is  thus  necessary  to  ensure  that  the  separations  at  P  between  the  different 

wavelengths  (of  separation  .'.A)  also  equals  Ap_.  The  grating  G?  performs  the 

separation  of  the  colors  at  P^.  These  colors  were  initially  spatially  separated 

horizontally  at  P  by  Ap  as  in  (4.6c)  and  later  recombined  at  P  by  condition 
A  a  B 

(4.5).  The  separation  Ap  at  P^  between  the  different  is  related  to  their 
wavelength  separation  ^  \  by 


A,> 


,Afodio 


(4.9) 


Since  each  vector  in  f  at  PA  is  not  just  one  A  but  contains  a  range  6a  of 

wavelengths  (due  to  the  finite  slit  width  Sp  used  (Fig.  4.2)  as  well  as  due  to 

practical  6A  ranges  for  LDs ) ,  the  Sp  spatial  spread  due  to  a  6A  is  given  by 

u 

(4.6c)  with  p  =  Sp  and  A  =  6>,  or  Sp  -  f’  (6A)d.  The  P„  output  will  also  have 

cl  cl  U  v 

a  spatial  spread  or  resolution  Sp  due  to  Sp  given  by  (4.6a)  with  d  =  d^,  and 

p  =  Sp  .  The  horizontal  image  resolution  in  P  due  to  6p  is  thus 

C  C  3 

6p  =  Sp  (d  /  d  )  (d,  +  d.,  +  d  )  /  (d  +  d,)  (4.10) 

c  a  1 U  4  1  /  i  1  z 

Experimentally  ip^  <_  0.1  mm  was  obtained.  With  hi)  sources,  typical  spectral  line 

widths  are  0.01  nm;  therefore  we  expect  about  a  10X  improvement  in  horizontal  image 

rose 1 ut ion . 

When  the  actual  system  of  Fig.  4.5  is  used  with  multiple  linear  LD  input 

arrays  at  P,,  lens  I,  is  placed  f  _  to  the  right  of  P  .  This  collects  all 

A  2  c  2  A 

light  from  the  LDs  and  images  each  LD  column  at  P  horizontally  across  P  ,  thereby 

A  B 

illuminating  P  with  white  light. 
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4  . 4  REVIEW  OF  COMPLEX  AND  BIPOLAR  DATA  HANDLING 

In  this  section,  we  review  the  major  techniques  by  which  complex-valued 
data  and  bipolar  data  can  be  handled,  represented,  proces  .cd  and  operated  upon 
in  a  noil-coherent  vector-matrix  ;  rocessor  of  the  type  used. 


Fig.  4.6  Schematic  of  basic  vector-matrix  non-coherent  APAR  processor. 


Consider  first  the  svstem  of  Fig.  4.6  whose  input  at  PA  is  a  1-D  vector  a 

a  n 

With  a  mask  b  at  P  ,  the  P,  output  is 
m .  n  b  o 


n  m=  1 


a  , 


(4.11) 


where  M  elements  at  i’  ,  N  X  M  at  P  at  N  at  P  ,  are  assumed.  The  P  output  in 

A  B  C  C 

(4.11)  is  the  vector  matrix  product 


=  K  a 


(4.12) 


where  (  )  denotes  a  vector  and  (  )  a  matrix.  Since  such  non-coherent  processors 


can  operate  only  on  positive  valued  data,  new  techniques  are  necessary  when  the 
data  values  are  bipolar  or  complex  (the  latter  is  the  case  in  APAR) .  We  first 
consider  two  methods  to  process  bipolar  data  and  then  two  techniques  for  complex- 
valued  data.  In  Sect.  4.5,  we  describe  a  new  and  preferable  method  for  complex¬ 
valued  data  processing  using  the  wavelength  diversity  processor.  In  Chapter  3, 
we  consider  a  new  method  when  multiple  wavelength  sources  are  not  available. 

This  latter  technique  is  the  one  used  in  our  IOP  system,  the  one  discussed  in 
Sect.  4.5  is  the  (",.e  used  in  the  wavelength-diversity  processor  (the  subject  of 
the  present  chapter) . 

in  the  first  bipolar-valued  data  processing  method  considered,  we  bias  x 

'  m 

and  K,,  data  as  (1/2  +  x  )  where  -1/2  4  1/2,  so  that  the  recorded  data 

ui ,  1 1  n  n 

a  and  b  satisfy  0  .<  a  ,  b_  _  <  1 .  The  output  from  detector  n  at  P  will 

then  be  the  sum  of  the  light  leaving  the  n-th  column  of  P  or 

B 

C  =  :  (1/2  +  1l  n)  (1/2  +  X  ).  (4.13) 

n  m  m ’n  m 

The  correct  y  output  (4.11)  can  be  obtained  from  c  by 
n  n 

v  =  <•  -  :■  (a  / 2  +  b  12)  +M/4.  (4.14) 

n  n  m  m  m ,  n 

The  2  a^/2  and  M/4  terms  in  (4.14)  can  be  optically  computed  by  addition  of  one 

extra  column  to  B  with  uniform  trans ni ttance  of  one-half.  At  the  correct  detector 

output  in  P  ,  we  obtain  1  a  / 2  +  M/4.  which  can  thus  be  subtracted  as  a 

C  mm 

constant  from  all  observed  c  outputs  in  an  electronic  post  processor.  Computing 

b  / 2  could  also  be  performed  optically  in  one  iteration  with  all  a  =  1/2. 

m  ti: ,  1 1  m 

The  c  output  in  P  is  then  21  /2  +  M/4.  This  technique  reduces 

the  useable  dynamic  range  of  the  I. FT)  input,  the  B  mask  and  the  output  detector,  but 
it  requires  electronic  post  processing  plus  an  extra  iteration  of  the  optical 
system  (or  electronic  computation  of  1  b  for  each  new  B) .  This  technique  only 
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i 

i 

\ 

i 

I 

t 


increases  the  space  bandwidth  (SBW)  of  li  to  (S  +  1)  (M  +  1)  from  NM  and  is  thus 
at  t rat t i ve . 


In  the  second  method  of  bipolar  data  processing  we  consider,  each  component 


is  divided  into  its  positive  (  )  and  negative  (  )  parts.  The  output  from  P 


■*+  T  _  A  + 

is  thus  [x  x  ]  where  the  top  half  of  the  LED  source  array  is  fed  with  x  and 
mm 


the  lower  half  with  x  .  Note  that  this  scheme  now  requires  a  2M  element  input 
m 


I.KD  arrav.  The  B  mask  at  l1  is  similarly  decomposed  into  four  sections  (increasing 


its  SBW  bv  a  factor  of  4  from  MN'  to  2M  *  2N)  as  in  (4.15).  The  c  output  at  P 

n  C 


is  thus 


—  “ 

—  - 

—  — 

f,  + 

_ 

H+  H~ 
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X 
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A  'N  element  linear  output  detector  is  thus  needed.  The  >'  outputs  appear  on 

_ 

the  first  \  detector  elements  and  the  v  outputs  on  the  last  N  detector  elements. 

n 

The  increased  SBW  of  l>  -P(-  is  •>  disadvantage  of  this  system.  However,  it  requires 
no  electronic  post  processing  and  does  not  reduce  system  dynamic  range. 


Specific  quant  i  I  icat  ion  of  t  lie  dvnamic  range  of  the  a  .  b  and  c 
11  '  m  m ,  n  n 

elemental  values  and  the  P  -P  .  component  specifications  (SBW  and  dynamic  range) 

A  l  ■ 

for  exact  scenarios  is  needed  to  enable  a  decision  between  these  two  bipolar 
data  processing  schemes  to  he  made.  We  now  consider  two  techniques  for  complex- 
valued  data  processing. 


\n  attractive  method  for  performing  vector-matrix  mu l t i p 1 icat ions  on 
,  onp 1 cx-va i  tied  data  was  described  bv  Bond  |  27  \  for  (XI)  processing  and  Goodman,  et 
,i  l  .  !  ,‘S  I  for  optical  processing.  It  was  then  applied  to  a  vector-matrix  elect  ro- 

oprical  processor  ro  perform 
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the  1)FT  by  Goodman,  et  al .  [29]  and  for  AFAR  processing  by  Psaltis,  et  al  . 

1 5 1 .  In  this  technique,  each  element  is  decomposed  into  its  residues  along 
unit  vectors  spaced  by  120  =  2’i/3  in  the  complex  plane  as 

x  =  Xy  exp  (  jr  0)  +  x  ^  exp(i2’i/3)  +  x.,  exp(j4ji/3)  (4. 16) 

A 

For  simplicity,  we  denote  the  three  components  of  x  along  the  three  directions 

-A 

by  Xq,  x-,  and  x ' ,  with  similar  notation  used  for  the  three  components  of  the  mntri 
^  and  the  output  vectors  "v. 


we  describe 

Fig .  4.6 

as 

r-~ 

j  -  6 

t 

Pq 

:  *i 

= 

Uj  .Uq  h2 

L;i 

2  H  *Q 

(4  .17) 


'2) 


In  (4.17),  all  components  are  real  and  positive.  The  P^  input  is  2M  long,  P 
is  3N  •  3M  and  P  is  3N.  The  P^  mask  must  be  properly  arranged  as  shown  in  (4.17) 
The  Increased  SBW  of  P^-P^.  by  factors  of  3  or  9  is  the  major  disadvantage  of 
such  a  svstem.  This  system  arrangment  was  the  one  we  used  in  our  earlier  phase  1 
work  [5 1  to  realize  and  demonstrate  electro-optical  processing  of  complex¬ 
valued  APAR  data. 


The  second  complex-valued  data  processing  technique  previously  described 
]6]  and  demonstrated  by  us  used  color-multiplexing  and  the  complex  decomposition 
in  (4.16).  In  this  processor  (see  Fig.  4.7)  three  linear  LD  or  LED  input  arrays 
are  used  at  with  different  wavelength  outputs  ».  -\0.  The  inputs  to  the  three 
1.1)  input  arrays  are  and  the  R  mask  at  P  was  arranged  as 


I!  = 


-hi  ;~*2  ^ 


(4.18) 
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The  B  mask  was  illuminated  with  all  input  wavelengths  as  before  and  the  three 
separate  as  were  then  separated  vertically  (not  horizontally)  at  by  a  grating 
between  P.}  and  P.^ .  The  output  P  patterns  on  the  three  linear  output  phot  odet  ect  or 
arrays  are  thus 


co  = 

C  On 

lb0; 

h 

O 

OJ 

031 

C1  = 

C  In 

fBo- 

*1 

Vai 

(4.19) 

-1. 

C2  = 

C2n 

lBo 

*1 

B2la2 

and  the  system  so  described  by 

ft 


The  aesired  v ^ .  v“,  outputs  are: 


yo 

=  H  -  x- 
-0  0 

+ 

l\x- 
-2  1 

+ 

a 

ti-X„ 
-1  2 

V  - 

1 

=  1!  -  x  - 
-1  0 

+ 

H~x- 
-0  1 

+ 

A 

H~X" 
-2  2 

V  - 
'  2 

=  H-x- 
-2  0 

+ 

H-x- 

-1  1 

+ 

H-x. 
-0  2 

'  *  ■» 

the  observed  c^, 

out  purs 

at  P, 

are : 

co  = 

L,VV. 

-iao  • 

— 2ao 

_» 

X 

1 

‘i = 

;  -0H1  ' 

-r'l : 

— 23 1 

_\ 

V 

_  > 

c  2  = 

1  B  a  ' 
0  2  • 

-la2 

-2a2 

=  B  [a0  ,  a2] 


(4.20) 


(4  .21) 


(4  .22) 


.a  realize  (4.21)  from  (4.22)  we  combine,  the  detector  array  outputs  as 


shown  in  Figure  4.B; 


— 

B  a 

+ 

B, 

+  B  a 

•  0 

0  0 

2 

—2 

_> 

_k 

y  = 

15  a 

+ 

15  a 

-r  15,,  a 

1 

V  1 

1  0 

-> 

_> 

_ 

V  „  = 

B  a 

+ 

B ,  a , 

+  B  ;i 

) 

0  2 

1  ! 

-  2  i 

The  system  of  Fig.  ..7  reduces  the  1-0  input  und  output  SBW  requirements 
by  a  factor  of  three  (by  use  or  three  1-0  input  l.!1  arravs  .it  different  ■  and 
three  1-1)  output  detector  arravs  stacked  vertically)  and  to-  SBW  ot  the  I*, 
mask  B  to  (  •  N  •  In  Sect.  -.r',  a  superior  co  1  ot  -mu  1 1  i  p  1 <. xed  processor 

without  electronic  post  processino  is  dese  r  i  bed . 


4  . 5  comit.i.:x-vai.it:i)  wavkixxcth-p i vers  i tv  pkqckssor 

In  this  section,  we  describe  new  and  preferable  bipolar  and  complex- 
valued  APAR  processors  and  techniques  using  the  wavelength-diversity  processor 
We  also  provide  experimental  confirmation  and  demonstration  of  the  system. 


The 


P  , 
A 


K'e  first  consider  the  new  bipolar-valued  data  processor  to  realize  (4.15) 

basic  system  of  Fig.  4.1  is  used  with  the  data  formatted  as  shown  in  Fig. 

->  x  _  + 

For  this  case,  we  decompose  x  into  x  and  x  and  H  into  11  and  II  .  At 
we  place  three  linear  1.1)  arrays  each  emitting  at  a  different 


and 


masks  side-bv-side,  B 


0 


We  describe  the  three  linear  l.D  outmits  at  P  ^  hv  h^,  ; 
,  a,  =  x  ,  and  a.,  =  x  .  At  1’  ,  we  place  two  matrix 

1  l  n 

ii+  and  Bj  =  H  .  The  output  at  l’(.  will  thus  have  four 


terms  spatially  separated  horizontally: 
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co 

’Vo 

ji 

C  , 

=  B,  a„ 

+  B  a 

1 

-1  0 

Yi 

_v 

c  _ 

=  B,a, 

+  B  a 

2 

-1  1 

-0  0 

C3 

-  v2- 

The  central  and  c^  outputs  at  P  are  the  desired  bipolar  components  of 

Hxr 

--  ^  — -+  +-•- 

y  =  c^=H_x  +Hx 

-*  — k 

v  =  c2  =  H  X  +  H  x  (4  .25) 

From  (4.25),  we  see  that  the  use  of  one  additional  linear  LD  input  array 

— »  —  j 

provides  the  correct  y  and  y  bipolar  output  components  on  the  two  parts  of  a 
linear  output  photodetector  array  in  P  .  Use  of  the  prior  color-multiplexing 
scheme  with  bipolar  data  would  require  electronic  post-processing  to  sum  and 
combine  the  proper  photodetector  outputs  as  well  as  two  linear  output  detector 
arrays  stacked  vertically. 


We  next  consider  a  second  method  of  processing  bipolar-valued  data  using 
wavelength  diversity.  This  approach  combines  the  methods  of  (4.14)  and  (4.15). 
This  technique  uses  two  linear  LD  input  arrays  (rather  than  three  as  above)  and 

-V  -*-f 

one  mask  (not  two  as  above).  The  LD  outputs  at  X ^  and  are  a^  *  x  and 

a  =  x  and  the  components  of  the  B  mask  are  b  =  (b  +  1/2).  The  resultant 
1  —  m  ,n  m  ,n 


outputs  on  the  two  parts 

of  a  linear 

photodetector 

array  at  P  are: 

y 

n 

=  b 

m  m,n 

aA  -E  ;1A 

0 ,  m  m  0 ,  m 

(2 

ji  _ 

V 

*  n 

=  ;:  b 
m  m  ,n 

a  a 

1  *  m  m  1 ,  m 

12. 

(4  .26) 

We  can  generate  the  last  two  terms  in  (4.26)  by  use  of  an  additional  column  of 
B  with  transmittance  1/2  as  noted  earlier.  The  outputs  from  the  two  associated 
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photodetectors  are  the  last  terms  in  (4. 26)  and  thus  can  be  electronically 
subtracted  from  the  observed  c  outputs  to  produce  the  desired  v  and  y 
outputs  as  in  (-*.26).  Thus  in  order  to  reduce  the  .5151.'  requirements  of  P^  and 

Pu.  we  trade  oft  dynamic  range  in  P^  (by  adding  a  bias)  and  re noire  some  electronic 
post  processing. 

We  now  consider  a  wavelength  diversity  processor  to  realize  the  complex¬ 
valued  vector-matrix  operations  required  in  APAR.  This  system  does  not  require 
the  electronic  post-processing  of  our  prior  system  [A  .  1)  nor  does  it  need  2-D 
or  stacked  1-D  output  detector  arrays. 


We  consider  forming 


v  =  H> 


(4  .27) 


where  all  elements  can  be  complex.  At  P  ,  we  place  five  (not  3)  linear  LD 


arrays  emitting  at  > 


0 


..nd  with  inputs  -  a ^ ,  where 


ao  =  xo’  a L 


Xi’  11 2  x2’  ;i3  X0’  a4  X1 


(4  .28) 


and  xq-x9  are  described  by  (4.16).  At  PR,  three  B  matrices  B H.,  are  placed 
side  bv  side  where 


V  *-1 


=  HJ' 


(4  .29) 


and  ll^-!!7  are  the  projections  of  (I  as  in  (4.16).  The  resultant  P(,  pattern 


contains  terms  in  five  wavelengths. 
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_x  ^  A  A 

At  Aq,  we  find  J3a0  containing  the  three  terms  BqSq,  JL^aQ  and  — 2ao  seParated 

by  Ap^,  which  was  chosen  equal  to  the  Ap^  separation  of  Hq,  and  at  Pg*  The 

\  ^ 

Ba  output  at  A,  is  shifted  right  by  Ap  from  the  Ba  output  with  the  Ba,  output 
—  11  c  u  t 

i  A  i 

shifted  right  by  Ap^  from  the  Ba^  output.  The  output  equals  Ba^  and  is  dis¬ 
placed  from  the  Ba^  output  by  Ap^  (all  A^+^-  ^  values  are  assumed  to  equal  AA 

x  a 

and  AA  is  assumed  to  correspond  to  a  Ap^) .  We  show  the  Ba^  “  Ba^  output  on  the 

x 

top  line  in  the  P^  pattern  in  Fig.  4.9,  displaced  by  Ap^  from  the  pattern 

i  X 

and  by  3Apc  from  the  BaQ  pattern  in  XQ  light.  The  Ba^  pattern  is  similarly 
shown  beside  the  Ba^  pattern  since  Ba^  =  Ba^ .  It  is  displaced  by  3Apc  from  the 
A^  pattern  for  reasons  noted  earlier. 


Fig.  4.9  Formats  for  the  novel  complex-valued  data  processor  for  APAR. 


157 


As  before,  the  actual  patterns  at  P^  are  superimposed  vertically 

and  are  shown  separated  vertically  only  for  clarity.  The  detector  forms 
the  sum  vertically  in  m  of  the  P  outputs  shown.  The  outputs  of  the  three 
central  detector  arravs  are 


Zk 
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=  v- 

=  H-x- 

+ 

H  -  X  - 

+ 

H  -  x  - 
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•  2 

-2  0 

-1  1 
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-  H-x- 

+ 
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=  II  X 

+ 

II -x  • 

+ 

H  •  x ' 

4 

•  1 

-  1  0 

-0  1 

—2  2 

We  note  that  each  of  the  output  components  (>’q,  v]  ,  and  yy)  in  (4.30)  contain  three 
•'  . .  Prom  (4.30)  we  see  that  the  outputs  from  the  I.  ■  N  =  3N  linear  output 
detector  elements  are  the  0,  1  and  3  projections  of  the  complex  y  vector-matrix 
product  in  (4.27).  Thus,  by  the  use  of  five  input  LD  arrays  (rather  than  3), 
no  output  P  electronic  processing  (summing  and  reorde. ing  of  data)  is  needed  to 
obtain  the  components  of  the  complex  output  (as  in  Fig.  A. 8). 


To  demonstrate  the  use  of  the  data  formats  in  Fig.  4.9  in  the  system  of 
Fig.  4.1,  a  specific  complex  vector  x  and  matrix  li  were  chosen: 


i- 

— 

1  +  exp  (j2n/3i 

0  +  1 

1  +  exp  (j4n/3) 

6  +  2 

1 

= 

6 

exp  ( j 2 n / 3 ) 

i 

exp  (  j 4 ti / 3 ) 

o 

U  -  J 

0+2  1  +  2 


(4.32) 


0  +  2  1  +  2 


Flio  short-hand  0,  1,  2  notation  is  used  for  simplicity  of  description.  We  write 

A 

the  three  0,  1,  2  components  of  x  and  H  explicitly  as 


xc  =  1 

u 


0  ,  x“  =  0 


(4.33) 


10  10  0 


'0  1  0  1  ol 


0  0  111 


0  10  0  1 


10  10  0 


0  1110 


10  0  10 


0  10  0  1 


1110  0 


(4.34) 


0  0  10  1 


01010. 


10  0  10 


0  0  10  1 


110  0  1 


10  0  11 


The  format  of  Fig.  4.9  was  thus  used  in  the  laboratory  version  in  Fig.  4.5 
of  the  system  of  Fig.  4.1  with  the  vector-matrix  data  of  (4.33)  and  (4.34).  The 
pattern  is  shown  in  Fig.  4.10a  for  the  vector  in  (4.33).  Each  of  the  five 
vector  inputs  shown  are  actually  in  different  colors  or  X  but  are  printed  in 


black  and  white.  The  P  pattern  containing  H;,  H-  and  H;  in  (4.34)  is  shown  in 

Fig.  4.10b.  The  output  pattern  is  shown  in  parts  in  Fig.  4.10c.  The 

output  in  colors  and  A  ^  is  shown  first.  Each  is  H  x^ .  Because  =  3A' 

corresponds  to  a  3p^  separation,  the  and  A  ^  pattern  contains  two 
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la,  I  '  i  I  ' 
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COMPLEX  VECTOR  OUTPUTS 


(b) 
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i  **  i  ^  i 

MJf/Vl 


tix6 


(C) 


iiii 

ii  ii  au  n  ■ 

ii  if  i  ii  ii  c 

i  i  i  i 

i  ■  i  ■ 


HXj 


Hxj 


ii  in 


composite 


Spoco  Integrating 
Dotector  Array 


(d) 


Fig.  4.10  Example  of  complex-valued  vector-matrix  multiplication. 
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H  Xq  patterns  side-by-side.  Since  x^  contains  three  I's,  only  the  first  three 
rows  of  H  appear  as  shown. 

Tlie  output  pattern  in  X.  and  A.  is  shown  next.  Since  a  =  a  =  x - , 

14  14  1 

-K  A 

this  pattern  is  two  11  x~  patterns  side  by  side.  Since  xj  contains  a  1  as  the 
first  and  fourth  elements,  only  the  first  and  fourth  rows  of  H  appear  in  the  H  X' 
output  as  shown.  The  11  x^  output  is  shown  next. 

As  noted  before,  these  patterns  are  superimposed  vertically  at  P  .  This 
composite  P  output  pattern  is  shown  in  Fig.  4.10d.  The  central  15  columns 
of  the  Fig.  4.10d  pattern  are  incident  upon  a  15  element  photodetector  array 
shown  in  Fig.  4.10e.  The  detector  integrates  each  column  in 

Fig.  ^ . lOd .  The  first  5  detector  outputs  in  these  15  elements  are  the  five  elements 
of  y2,  the  central  5  detector  outputs  are  the  elements  of  y^  and  the  right  5  detector 
outputs  are  the  elements  of  the  "y^  component  of  the  complex  y  =  H  x  output. 

When  the  system  of  Fig.  4.1  is  used  as  an  iterative  processor,  a  5  element 
complex  vector  with  components  0,  1,  2  is  added  to  the  respective  15  detector 
outputs  and  fed  back  to  PA  as  the  6,  1,  2  components. 

To  show  that  the  experimentally  obtained  output  is  correct  for  the  specific 
example  in  (4.33)  and  (4.34),  we  computed  y  =  H  x  and  found 


y  -  iy0,  yx,  y2.  y3.  y4J 

=  [40  +  3 i  +  22,  20  +  31  +  52,  40  +  41  +  32,  20  +  11  +  52,  40  +  31  +  22]  (4.35) 

=  [20  +  i,  i  +  32,  6  +  i,  i  +  32,  20  +  i]. 

The  first  representation  in  (4.35)  lists  the  five  elements  of  y.  The  computation 

... 

of  y  =  li  x  yields  the  second  expression  in  (4.35).  Since  0  +  1  +  2  =  0,  we 


subtract  these  common  factors  and  obtain  the  last  expression  in  (4.35). 
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The  detector  outputs  from  left  to  right  are 


y  =  y6.  yj  -  y20"-y24  yoo‘"yo4  yio" *yi4‘ 


From  (4.35),  the  expected  P  detector  output  is 


Y  =  Y,  Vq  >'j  =  [03030201021111  1]. 


(4  .36) 


(A  .37) 


The  optical  output  data  (obtained  by  summing  the  column  outputs  in  Fig.  4.10d) 


25352  |  42424  |  33433  » *y *  "y'.'y.. 

i  i  2  U  1 


(4.36a) 


Combining  the  first  elements  of  y^,  y^  and  y '  and  then  the  second  etc  elements. 


we  obtain 


243  523  344  523  243, 


The  symmetry  of  these  outputs  was  set  by  the  choice  of  matrix  used. 

Subtracting  the  offset  or  lowest  value  from  each  pair  of  three  values  (e.g.  for 
the  first  set  232,  the  lowest  value  is  2.  Subtracting  2  from  each  element 


of  this  set  yields  012)  yields 


021  301  011  301  021. 


(4.36b) 


We  now  recall  that  the  three  elements  in  each  of  the  above  groupings  corre¬ 
spond  to  the  201  component  values  respectively,  we  write  the  above  output  as 


2  0  +  1  3  2  +  1  0  +  1  3  2  +  1  2  0  +  1. 


(4.37) 


This  detector  output  pattern  is  shown  in  Fig.  4.10f  alter  electronic  post  pro¬ 
cessing.  This  pattern  in  (4.37)  agrees  with  the  predicted  result  in  (4.35), 
thus  verifying  the  system's  operation  as  a  vector-matrix  multiplier  of  complex¬ 
valued  data  forAPAR.  The  subtraction  of  the  smallest  component  of  each  y 
element  output  from  all  five  of  those  output  components  is  easily  achieved 
bv  a  comparato r/ subt ract or .  This  is  necessarv  because  the  components 


0,  1,  2  are  not  linearly  independent.  A  simple  geometrical  construction  to 


allow  computation  of  the  components  of  the  complex  data  is  shown  in  Fig.  4.11. 
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4.6  MULTICHANNEL  1-D  CONVOLVER 


An  alternate  and  more  general  formulation  of  the  wavelength  diversity 

processor  of  Fig.  4.1  as  a  multi-channel  1-D  or  vector-matrix  convolver  is 

possible.  We  can  advance  this  formulation  with  reference  to  Fig.  4.4.  We 

denote  the  L  linear  vectors  at  P.  by  the  index  V. ,  the  R  matrices  at  P„  by 

A  B  1 

the  index  r,  and  the  L  linear  N  element  outputs  at  P  by  the  index  •  also. 
Summing  the  vector  outputs  at  9,-2  in  P  ,  we  obtain 


Vo  +  Vi  +  v> 


L0  2 


(4.38) 


or  in  more  general  notation  with  f  and  r  as  the  a  and  B  subscripts, 


c  =  L  B  a 
£  r  ~r  i-r 


For  ■  =  2  and  r  =  0,  1,  2,  (4.39)  becomes 


"2  =r=0  — r^2-r  =  V2  +  Vl  +  V<T 


(4.39) 


(4.40) 


By  inspection  of  (4.39)  and  Fig.  4.4,  we  see  that  the  form  of  the  output  at  a 
given  N  element  detecior  in  P  is  that  of  a  type  of  convolver,  specifically  a 
mu  1 1 i-ehanne L  1-D  vector-matrix  convolution.  This  general  operational  descrip¬ 
tion  in  (4.39)  of  the  system  improves  the  power  and  flexibility  of  the  wave- 
length-dive  rs  i  tv  processor  witli  wavelength  *  as  the  convolution  shift  variable. 
We  thus  rewrite  (4.39)  as 


c(>.)  =  /  b(r)  a  (  /  -  r )  d  r 

c ( \ )  =  /  n(r)  a(X-r)dr,  (4.41) 

from  which  the  convolution  with  as  the  shift  variable  is  apparent.  The  P 


output  is  shown  in  Fig.  4.12  for  this  general  case. 


/ 
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4.7  MATRIX-MATRIX  MULTIPLICATION 


In  this  section,  we  consider  a  new  application  of  the  wavelength  diversity 
processor,  to  perform  a  matrix-matrix  multiplication.  We  describe  two  new 
electro-optical  systems  that  achieve  this  operation.  Consider  formation  of 


C  =  B  A. 


(4.42) 


where 
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Substitution  shows 
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to  he  t!,<  desired  result. 


Consider  the  simplified  version  of  Pig.  4.1  shown  in  Fig.  4.13.  The  input 
at  !’  is  N  vectors  a  each  of  length  M  corresponding  to  N  spatially  separated 

1  Ml  i  1.1  ) 

wavelengths  ■  .  L.u-h  linear  Id)  vector  is  imaged  verticalv  and  spread  horizon- 

n 

tally  across  1’,,.  where  .1  mask  B(m,n)  is  placed.  1*2  2s  Imaged  onto  I’ 3,  while 
each  'n  is  separated  vertically 
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at  P^  by  a  grating  as  before.  The  P  pattern  is  thus  described  by  N  vectors 

c  (n)  eacli  of  length  N  and  each  formed  in  a  different  1  light  where  p  »  '  f  z, 
pn  n  s 

l  is  the  input  wavelength,  f  is  tiie  spatial  frequency  of  the  grating 
between  P>  and  P.^  and  z  is  the  distance  from  l\;  to  P  )  . 

To  realize  (4.44)  using  Fig.  4.13,  we  place  A  at  as  shown  in  Fig.  4.14. 

T 

When  the  R  mask  is  placed  at  1*2,  the  P.^  output  containing  eight  elements 
.■suits.  The  detector  array  at  P^  contains  four  elements  as  shown.  Over- 

layiu.  the  detector  and  outpuL  patterns,  we  see  that  each  detector  integrates 
the  sum  of  the  .  orresponding  two  output  P^  components.  The  four  detector  out¬ 
puts  are  thus  the  components  of  the  desired  C  output  matrix  in  transposed  form. 

We  now  consider  a  matrix-matrix  multiplier  using  a  1-D  detector  and  1 -I) 
input  vectors  of  longer  length.  To  visualize  this  system,  recall  the  descrip¬ 
tion  of  Fig.  4.1  in  (4.1), 

M 

c  =  ::  b  a  (4,i) 

■n  n,  mn  m. 

This  is  the  product  of  the  two  matrices  A  and  B, 

C  =  BA,  (4.45) 

where  A  is  an  I,  X  M  matrix  and  B  is  N  X  M. 

To  demonstrate  this  use  of  the  system  of  Fig.  4.1,  we  consider  the  matrices 

- m  — — —  ►  n 
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Substitution  into  (4.45)  yields 


3  12  2  1 

12  2  13 
2  2  13  1 

2  13  12 

13  12  2 


(4  .47) 


To  perform  the  matrix-matrix  product  in  (4.46)  and  (4.47)  on  the  system  of 

Fig.  4.1,  the  experimental  system  of  Fig.  4.5  was  used  with  five  different  input 

wavelengths  to  distinguish  the  five  rows  of  A.  The  input  V  pattern  is  shown  in 

A 

the  top  of  Fig.  4.15,  in  actuality  this  image  was  in  color  with  each  column  being 
a  different  color.  This  2-D  input  pattern  with  color  horizontal  with  index  > 
and  with  five  elements  per  column  with  index  m  describes  the  2-D  matrix  A(>,m). 
Incident  on  Pn,  the  five  rows  are  illuminated  with  different  combinations  of 

D 

input  wavelengths  as  shown  in  the  center  of  Fig.  4.15.  The  P  mask  B(n,m)  in 

D 

(4.46)  is  shown  to  the  center  right  in  Fig.  4.15.  The  output  plane  pattern 
at  P  is  shown  in  the  bottom  of  rig.  4.15. 


The  detector  at  P  sums  the  five  entries  in  each  column  in  P  .  The  25 
detector  outputs  (from  inspection  of  the  bottom  figure  in  Fig.  4.15)  are  (from 
1 ef t-to-r ight ) 
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3122112213221312131213122. 


(4.48) 


We  see  that  these  detector  outputs  are  the  elements  of  the  matrix  C  in  (**.47) 
thus  proving  the  system's  use  as  a  matrix-matrix  multiplier. 
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4.8  MATRIX- INVERSION' 


We  now  consider  yet  another  new  application  of  the  wavelength-diversity 

processor,  its  use  in  inverting  a  matrix.  This  applieaiton  is  of  direct  use 

in  AFAR  processing.  Consider  the  matrix-matrix  multiplier  of  Fig.  9.13  as 

described  in  Sect.  4.7  with  the  feedback  shown  applied.  We  describe  this  10P 

matrix  processor  with  P  input  W.  ,,  PD  mask  f H ] ,  output  W..  After  addition 

A  1  “l  D  1 

of  a  matrix  F  to  the  output,  the  IOP  system  is  described  by 

W  =  H  W  +  F,  (4  .49) 

—i+l  -  -i  — 

where  all  quantities  are  matrices  and  the  subscripts  refer  to  successive  output 

iterations.  In  steady  state  W.  =  W.  ,  ,  =  W  and  (9.49)  becomes 

-i  “i+l  — 

W  =  H  W  +  F 

=  F(J-H)'1  (9.50) 

If  we  let  H  =  I-M  and  F  =  J_,  then  the  output  is 

W  =  M'1.  (9.51) 

The  system  has  thus  performed  the  inversion  ot  the  matrix  M  at  P^.  The 
system  and  algorithm  used  are  shown  in  the  functional  diagram  in  Fig.  9.16. 

Note  that  the  output  matrix  to  be  added  is  the  identity  matrix  and  thus  this 
potentially  complex  operation  is  greatly  simplified. 

W-£L  W  +  I  =  W 
M  W  =  I 

fo  =  m 

Fig.  4.16  Functional  di  'gram  of  the  matrix-inversion  processor. 
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4 .9  COVARIANCE  MATRIX  COMPUTATION 

A  difficult  processing  step  in  the  AFAR  system  involves  the  calculation  of 
the  covariance  matrix,  M(?,n)  of  the  noise  field,  namely: 


M(P,n)  =  Xy(t)x*(t)  =  1  /  x(,  ( t )  x*  ( t  )df 


T  T 


(4.52) 


where  x.,(t)  and  xn(t)  are  the  noise  voltages  received  by  the  ?th  and  nth  antenna 
elements.  Rewriting  (4.52)  as  a  summation  over  t  we  obtain: 


M(J,n)  "  I  x(t,f)x*(t,n). 
t 


(4.53) 


We  note  that  (4.53)  is  identical  in  form  to  (4.1)  and  therefore  conclude  that  in 

principle,  the  matrix-matrix  multiplier  can  be  used  to  compute  the  covariance  matrix 

by  writing  the  time  history  of  each  of  the  f,  n  =  0,  1,  2  •  •  •  N-l  received 

signals  along  the  m  spatial  coordinate  in  P^  and  where  x^(t)  is  written  along 

the  ■ th  column  in  P,  and  x*(t)  is  written  along  the  nth  column  in  P„ .  Since  the 
An  b 

received  signals  are  complex  valued,  the  SBWP  required  in  PA  and  to  represent 


the  data  must  be  3X  as  large  and  the  resultant  output  P^,  SBWP  must  be  9X  larger 

than  that  required  to  perform  a  matrix-matrix  multiplication  of  positive  valued 

data.  The  P  ,  P  and  P  data  formats  would  be: 

A  o  (, 


where  N  is  the  number  of  antenna  elements  and  2N  is  the  minimum  number  of  time 
samples  required  to  estimate  the  covariance  matrix  M  as  a  time  average. 

Since  the  above  calculation  is  essentially  a  matrix-matrix  multiplication, 
an  experimental  demonstration  was  not  conducted.  Further  experimental  work  is 
recommended  in  this  area  in  order  to  assess  the  accuracy  of  the  overall  APAR 
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4  .  10  SUMMARY 

In  this  chapter,  we  have  described  a  new  wavelength-diversity  processor. 
The  use  of  this  system  as  a  vector-matrix  multiplier  and  iterative  processor 
with  complex-valued  and  bipolar  valued  data  elements  was  presented.  Three  new 
optical  systems  using  three  new  algorithms  to  represent  and  process  complex 
and  bipolar  data  in  non-coherent  electro-optic  vector-matrix  processors  were 
dose  r ibed . 


The  system's  use  as  a  matrix-matrix  multiplier,  a  matrix  inverter,  a  1-D 
multichannel  convolver  and  for  computation  of  the  covariance  matrix  M  were  also 
described.  In  each  case,  new  system  data  arrangements  were  used  to  realize  new 
and  more  powerful  electro-optical  processors  and  new  adaptive  algorithms  for 
various  APAR  processing  applications. 

A  laboratory  wavelength-diversity  electro-optical  processor  was  assembled. 

It  was  used  to  experimentally  demonstrate  and  verify  its  use  in  complex-valued 
vector-matrix  multiplication  and  as  a  matrix-matrix  multiplier.  In  all  instances, 
experiments  and  theory  agreed.  The  system  offers  the  promise  of  increased 
processor  capacity  (due  to  the  added  degree  of  freedom  that  wavelength  provides) 
and  is  of  use  in  developing  new  algorithms  (as  the  systems  described,  fabricated 

and  demonstrated  have  shown). 

The  highlights  of  the  WPP  approach  to  APAR  and  the  new  research  performed 
during  the  past  year  include: 

(1)  WDP  concept  and  demonstration  of  it. 

(2)  Use  of  the  WDP  concept  in  a  new  complex-valued  data  representation 


( Sec.'  t  .  4.4). 
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(3)  Modification  of  the  original  WDP  concept  to  decrease  the  electronic 
post-processing  necessary  and  demonstration  of  this  concept  (Sect. 
4.5). 

(4)  Laboratory  WDP  system  formulation  and  demonstration  using  an  art- 
la,  np  and  grating  (Sect.  4.3). 

(5)  Multi-channel  convolver  description  and  demonstration  of  the  WDP 
system  (Sect.  4.6). 

(6)  Matrix-matrix  multiplication  description  and  demonstration  of  the 
WDP  system  (Sect.  4.7). 

(7)  Description  of  the  use  of  the  WDP  system  for  matrix  inversion  (Sect. 
4.8)  . 

(8)  Description  of  the  use  of  the  WDP  for  covariance  matrix  computation 
(Sect .  4.9). 


CHAPTER  5 


-  i  7  j  - 

SUMMARY  AND  CONCLUSION 

We  summarize  the  highlights  of  our  programs  in  the  area  of  coherent  optical 
correlator  (COC) ,  iterative  optical  processor  (IOP)  and  wavelength  diversity 
processor  (WUP)  and  thei"  use  in  APAR  processing. 

Our  new  COC  work  on  APAR 'processing  has  resulted  in  many  new  algorithms  and 
system  architectures  as  well  as  in  the  fabrication  of  many  components  for  the  COC 
processor.  In  the  past  one  year  we  have: 

(1)  Developed  a  new  COC  concept  using  1-D  AO  cells  rather  than  2-D  SLMs 
(because  they  are  more  easily  fabricated  and  more  readily  available) 
and  TI  rather  than  SI  correlators  (because  their  longer  integration 
times  allow  better  noise  statistical  estimates)  (Sect.  2.2). 

(2)  Fabricated  and  performed  initial  testing  of  two  AO  cells  (Sect.  2.3). 

(3)  Developed  a  new  simulator  to  handle  multiple  phased  array  signals 

with  adaptivity  in  space  and  frequency,  to  compute  weights  and  calculate 
corrected  antenna  patterns,  plus  conventional  Fourier  transform  and 
correlation  routines  (Sect.  2.5). 

(4)  Devised  a  new  adjunct  antenna  concept  that  allows  full  space  and  fre¬ 
quency  APAR  data  to  be  obtained  using  only  a  two-channel  processor 


(Sect  2.1). 
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(5)  Developed  a  new  hybrid  Cime  and  space  integrating  AO  correlator 
architecture  that  combines  the  best  features  of  the  T1  and  the 

SI  systems,  and  provides  a  2-D  display  of  the  angle  and  frequency 
location  of  the  far-field  antenna  noise  pattern  from  the  two- 
element  adjunct  antenna  data  (Sect.  2.6). 

(6)  Designed  and  fabricated  a  hardware  electronic  support  system  and 
designed  and  nearly  completed  fabrication  of  a  multi-purpose 
computer-driven  electronic  support  system  for  complex  noise 
source  scenarios  (Sect.  2.7). 

(7)  Developed  a  self-correcting  post-processing  system,  a  new  fast 
iterative  modified  projection  method  technique  and  an  advanced 
projection  concept  to  compute  the  adaptive  weights  for  narrow- 
band  and  wide-band  jammers  that  differ  in  both  angle  and  fre¬ 
quency  (Sect.  2.8). 

(8)  Performed  an  experimental  verification  of  a  TI  correlator  and 
used  it  to  demonstrate  residue  arithmetic  computations  in  a 

new  time  position  coding  scheme.  Residue  arithmetic  is  an  advanced 
technique  whereby  numerical  computations  can  be  performed  in 
parallel  with  no  carries  and  with  reduced  dynamic  range  require¬ 
ments.  This  is  necessary  to  realize  an  accurate  optical  computer 
(Sect.  2.9) . 

(9)  Demonstrated  the  use  of  the  TSI  system  for  computation  of  the  2-D 
space  and  frequency  output  antenna  pattern  from  an  adaptive  array 
using  an  adjunct  antenna  (Sect.  2.10). 

Thus,  in  conclusion,  we  have  devised  a  new  COC  optical  processor,  demonstrated  and 
simulated  the  key  points  of  the  system,  begun  development  of  the  necessary  post¬ 
processor,  and  fabricated  the  necessary  components  and  support  svstem  to  fully 
study  this  COC  technique  for  APAK. 
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The  IOP  concept  for  APAR  processing  has  proven  to  be  most  excellent.  The 
highlights  of  our  recent  work  on  this  approach  to  APAR  processing  are  summarized 
below. 

(1)  Design  and  fabrication  of  a  new  IOP  system  with  fiber  optic 
interconnections,  pulse  width  modulation  and  a  microprocessor 
electronic  feedback  system  (Sect.  3.3). 

(2)  Use  of  a  new  technique  to  allow  the  IOP  to  operate  on  complex¬ 
valued  data  (Sect.  3.2.4). 

(3)  Development  and  demonstration  of  a  new  IOP  simulator,  antenna 
model,  and  2-D  to  1-D  mapping  with  space  and  time  adaptivity 
and  including  IOP  error  sources  (Sect.  3.4). 

(4)  Demonstration  and  measurement  of  the  excellent  accuracy  of  the 
IOP  system  to  be  less  than  0.8%  and  with  less  than  1%  error  in 
the  antenna  SNR  obtained  (Sect.  3.5.1). 

(5)  Development  and  demonstration  of  a  new  technique  to  ensure 
convergence  to  the  steady  state  solution  and  fast  convergence 
of  the  IOP  system  (Sect.  3.5.2). 

(6)  Real-time  demonstration  and  analysis  of  the  laboratory  IOP  system 
for  antenna  adaptivity  in  space  and  time  (Sect.  3.6). 

The  weights  and  antenna  SNRs  obtained  experimentally  were  within  1%  of  the 
theoretical  limit.  This  excellent  system  performance  in  the  IOP  laboratory  system 
fabricated  represents  a  major  new  optical  data  processing  architecture  that  appears  to 
be  quite  attractive  and  useful  for  APAR  processing  and  other  applications. 
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The  highlights  of  the  WPP  approach  to  APAR  ana  the  new  research  perrorrc 
during  the  past  year  include: 

(.  1 )  WDP  concept  and  demonstration  of  it. 

(2)  Use  of  the  WDP  concept  in  a  new  cor'p le:-:-va i ued  data  representation 
(Sect.  4.4)  . 

(3)  Modification  of  the  original  WDP  concept  to  decrease  the  electronic 
post-processing  necessary  and  demonstration  of  this  concept  (Sect. 
4.5). 

(4)  Laboratory  WDP  system  formulation  and  demonstration  using  an  arc 
lamp  and  grating  (Sect.  4.3). 

(5)  Multi-channel  convolver  description  and  demonstration  of  the  WDP 
system  (Sect.  4.6) . 

(6)  Matrix-matrix  multiplication  description  and  demonstration  of  the 
WDP  system  (Sect.  4.7). 

(7)  Description  of  the  use  of  the  WDP  system  for  matrix  inversion  (Sect. 
4.8). 

(8)  Description  of  the  use  of  the  WDP  for  covariance  matrix  computation 
(Sect.  4.9) . 
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